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Abstract 


In this study, we define some tridigional matrices depending on two real 
parameters. By using the determinant of these matrices, the elements of 
(s, t)-Pell, (s,t)-Pell Lucas and (s,t)-modified Pell sequences with even or 
odd indices are generated. Then we construct the inverse matrices of these 
tridigional matrices. We also investigate eigenvalues of these matrices. 


1 Introduction and Preliminaries 


Special integer sequences are encountered in different branches of science, art, 
nature, the structure of our body. So it is a popular topic in applied mathematics. 
Two of the special integer sequences are the Pell and Pell Lucas sequences. 
By changing the initial conditions but preserving the recurrence relation the 
Pell Lucas sequence is obtained. The authors investigated some sum formulas 
for Pell numbers in [2]. Gulec and Taskara generalized the Pell and Pell 
Lucas numbers by using two parameters in [3]. The authors generalized the 
modified Pell sequence similarly in [4]. By the determinant of the tridiagonal 
matrix, the values of the Fibonacci and Lucas numbers are demonstrated in 


[8]. Feng gave Fibonacci identities via determinant of the special tridiagonal 
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matrix in [9]. Seibert and Trojovsk studied the factorization of the Fibonacci 
and Lucas numbers using determinants of tridiagonal matrices in [10]. In [11], 
it is demonstrated that the permanents of some special tridiagonal matrices are 
equal to Fibonacci numbers. Falcon displayed some equalities of k-Fibonacci 
numbers with the determinant of some special matrices in [12]. Catarino obtained 
the n-th elements of k-Pell, k-Pell-Lucas, and modified k-Pell sequences by the 
determinants of some tridiagonal matrices in [13]. In [14], properties of hyperbolic 
generalized Pell numbers are studied. The recurrence relations for Pell, Pell 
Lucas and modified Pell sequences are P, = 2Pn—1 + Pn-2, Po = 0, Pi = 1, 
Qn = 2Qn-1 + Qn-2, Qo = 2, Qi = 2 and qn = 2Gn-1 + G-2, M=1, M=1 
for n > 2, respectively in [1]. There are some generalizations of these integer 
sequences. For example, the generalizations for Pell, Pell Lucas, and modified 
Pell sequences called (s,t)-Pell, (s,t)-Pell Lucas, (s,t)-modified Pell sequences 


are defined by the aid of the following recurrence relations respectively 


(n(s,t) = 28n—1(5, t) Tr tn—2(s, t), 99 (s, t) = 0, g1(s, t) =1, 
Rn(s,t) = 2%WsRn-1(s,t) + tRn—2(s,t), Rol(s,t) = 2, Ri(s,t) = 2s, 
Nn(s,t) = 2sX,-1(s,t) + tXp_a(s,t), No(s,t) =1, Nil(s,t) =s, 





for n > 2 in [8, 4]. 
Some elements of (s,t)-Pell and (s,t)-Pell Lucas sequences are given in the 


following tables: 


n (s,t) — Pell numbers 

1 1 

2 2s 

3 As? +t 

4 85? + 4st 

5 16s* + 12st + t? 

6 328° + 32s%t + 6st? 

7 648° + 80s*t + 24st? + £3 
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n (s,t) — Pell Lucas numbers 

1 2s 

2 Ag? 5 OF 

3 88° + 6st 

4 16s* + 16st + 2t? 

5 328° + 40st + 10st? 

6 645° + 96s4t + 36st? + 23 

7 12887 + 224s°t + 112st? + 14st? 


The elements of (s,¢)-modified Pell sequences are obtained by taking half of 
the (s,t)-Pell Lucas numbers. 

Special integer sequences are obtained with special numerical choices for 
(s,t)-Pell and (s,t)-Pell Lucas, (s,t)-modified Pell sequences. For example, if 
s =t = 1, then we get well-known Pell, Pell Lucas, modified Pell sequences. If 
s=1,t=hk, then we get k-Pell, k-Pell Lucas, and modified k-Pell sequences. 


Let us consider a tridiagonal matrix as 





ee | 
c de 

| c de | 

Hire | m | 

| c de 

cd] 
Then 

det A, = ddet An—1 — ce det An_2. (1) 


The inverse of a matrix A can be obtained by the formula A~! = toon" 
where (cof(A))" is the transpose of the cofactor matrix A or adjugate matrix of 


A [4]. Let T be a nonsingular tridiagonal matrix as 
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at by 
Cy ag bo 
T= | C2 ; | 
; bn—1 
| Cn-1 an | 


Usmani [6] gave a formula for the inverse of this matrix T~! = (t;;) as 


—1)'19 7 C5.-C4—10;-19541 if 2 >J 


ft (-1)') £0j...b;-1-10)41 fi <j 
7) 
( os 


where 


e 9; verify the recurrence relation 0; = a;0;-1 — b;-1q_10;_-2 for i = 2,....n, 
with the initial conditions 0) = 1 and 6; = a,. Observe that 6, = det (T). 





ed; verify the recurrence relation ¢; = a;¢;41—bicg;¢j42 for t=n-—1,...1, 


with the initial conditions ¢@n41 = 1 and dy = ay. 


In [7], if the tridiagonal matrix is given in the following form 





a 
a 0b 
T= Cc a os . 
b 
c a 
then the eigenvalues of this matrix are 
rT 

Ay = a + 2Vbecos(———), PEA 2 es Oi (3) 
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1.1 Some properties of tridiagonal matrices A,,, by (s,t)-Pell 


sequence 


Theorem 1. Assume that Apn is ann x n tridiagonal matrix defined as 


is t 
—1 2s t 


| —l 2s 


Apn = | . 
be Hy 6 
—1 2s | 
Then the determinant of Ap» is 


det(Apn) = @n41- 


Proof. The proof is made by induction applied on n. For n = 1, we have 
det(Api) = 92 = 2s. Assume that det(Apn—1) = @n, det(Apn) = 9n+41 for 
n > 2. Then 


det(Apni1) = 2sdet(Apn) — (—t) det(Apn_1) 
= 28On+1 + tn = On+2- 


(s, t)-Pell sequence is also obtained by using the following tridiagonal matrix 


with complex entries. Assume that A, is an n x n matrix defined as 


ie it 
2 Qs at 


a |, 
- | 
| | 
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Then it is easily seen that the determinant of A, is also (n + 1)th element of the 
(s, t)-Pell sequence 
det(An) = Qn+1- 


For the inverse of Ap», by using (2), it is obtained that 
aj = 25, bb =t, G =-—1, 
9, = Pi41, $7 = Pn—j+2- 
Therefore the inverse of Ap, is displayed by 
4 i Dee aa jitpeg eee 
Apales = tats 
99} On—i+1 5 if i>j 


aoe 
The elements of the cofactor matrix are given as 


cof (Ap 2) §4{9n—j+1, chip J 

a (-1)!YtJojOnin, if i>g 
It is well-known that |cof(Apn)| = |adj(Apn)| = [Age = Te By using 
cofactor matrix, we get some properties of (s,t)-Pell sequence. For n = 2, we get 


991 (92 9191 


|cof(Ap,2)| = 
—tpi91 £291 








9% + tot = 3. For n = 3, we get 


(21 (93 (21992 P1991 
|cof(Ap3)| =| —te1p2 262 g2—1 
Poigr —tyep1 9391 


|cof(Ap3)| = 9193(—53—1 + t—3~7) 
+t9102(01 032 + tes) + tet (1 e3 — 705) 
= 19593(3 + ter) + te? 03(03 + te1) 


= (93 + te)’ eI "2 
03 — toy 
= g{(93+ iia) 
2 (83 = aay 
ae 
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For n = 4, we get 


§1 994 971993 (192 1971 
—tp193 9293 9292 9201 
Peig. —te2p2. 9392 ~3—1 

oie, Pee —tese1 ~p4pi 


|cof(Ap,a)] = 


93 = (t@2 + ©3)(to1w3 + 21)”. 

Eigenvalues of the matrices Ay, construct the spectra of the Ap». By using 
the property (3), the sequence of the spectra of Ap, for n = 1,2,3,5 is 

iS hs ApS {28} 

n=2=>, = {28+ ivt, 2s —iV/t} 

n=3=> dA, = {28 + iV 2t, 2s — i/2t, 28} 

n=5=> rA» = {28 + i22V4, 2s + iV3t, 2s, 2s — 124, 2s — iV/3t}. 

If s =t = 1, the sequence of the spectra of the matrices Ap, for n = 2,3, 4,5,6 
is computed by using the Matlab Program as 

So ={2+i, 2-i}, 

$3 = {2+ /2i, 2, 2- V2}, 

oie 2 + 1.6180339887498967, 2 + 0.618033988749895i, . 

2 — 0.6180339887498957, 2 — 1.618033988749896i 


rae 2 + 1.732050807568879i, 2 — 1.732050807568879i | 
ay ee ee 
2 + 1.8019377358048387, 2 — 1.801937735804838i, 
Se= < 2+0.445041867912629i, 2 — 0.4450418679126294, 


2 + 1.2469796037174687, 2 — 1.2469796037174687 
Evidently, the product of eigenvalues is the determinant of the matrix and the 


sum of eigenvalues is the trace of the matrix. Therefore 


Sod = tr(Apn) = 2sn 
i=l 


and 
n 


[> = det (Ap,n) = Pn41 = [es + V2i cos(— 


j=l 


™) 
+1 





)). 
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1.2 Some properties of tridiagonal matrices E,,,, by even (s,t)-Pell 
sequence 


Assume that Ep »(p) is an n x n tridiagonal matrix defined as 


28 0 
t 4s? + 2t t 
| t 4s*4+2t °°. | 
Fyn = | ; ; | ‘ 
t 
t 4s? + 2t | 
Then the determinant of Ep. is computed by (1) as 
det Lyn = 2n- 


For the inverse of Ep», the values in (2) are computed as 


ay = 2s, 

Gy = Ag? OF, > 2 

bj = t, i>] 

b = 0, 

Cp SE; 

6, = ga, t21 

oy = 5 P2n-j42) peel 


Therefore the inverse of Ep, is displayed as 





0, if i=1 
(Epn)(is) = (—1)*9t9* o_9692(n—54:1) eaae if i<j 
(—1)*4t*9 993-29 2(n-i41) tapas i a oeg 





If all entries of the matrix are real and nonnegative, then the matrix is called 


positive. All eigenvalues are real if the matrix positive and tridiagonal [4]. 
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Therefore all eigenvalues of Ey», are real if s,t > 0. If we choose s = ¢t = 1, 
then the sequence of the spectra of the matrix E,, for n = 2,3,4,5,6 is given in 
the following result with the help of the Matlab program 

So = 12.6) 

S3 = 12,5, 7} 

S4 = {2, 4.5857864376269042, 6, 7.414213562373092} 

Ss “| 2, 4.3819660112501022, 5.381966011250107, 

6.618033988749892, 7.618033988749891 

Se = {2, 4.267949192431121, 5, 6, 7, 7.732050807568871} . 

Evidently 

>> Ai = tr(Epn) = (n — 1)(48? + 2t) + 28 and [] Aj = det(Eyn) = 2n- 

If we take care of the spectra, one of the eigenvalues is 2 = 2s ifs =t=1 
for all positive integer n. And the minimum eigenvalue of spectra converges to 


2 = 2s, the maximum eigenvalue of spectra converges to 4s? + 4t. 


1.3. The properties of tridiagonal matrices O,,,, by odd (s,t)-Pell 


sequence 


Assume that Op» is an n x n tridiagonal matrix defined as 


As? +t t 
t Ag? Ot t 
On = t As? + 2t 
t 
P As ape 


Then the determinant of O,,, is given by (1) as 


det Opn = 92n+1- 
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For the inverse of Op», the values are computed by (2) as 





a, = 4s*+t, 

a; = 49°42t, i>2 

bb = Gg=t, i>1 
G6= Pisi, t21 

oj = 5 Pini) ae 





(O52) 49) = (— 1)" 925-1 @2(n—41) Tapa’ if i<j 
(—1)*9t°9 995-1 92(n-i41) Zapp? if hes 





Matrices Op, are symmetric so the eigenvalues are real. If s = t = 1, the sequence 
of the spectra of the matrices Op» for n = 2,3,4,5,6 is given in the following: 
So = {4.381966011250105, 6.618033988749895} 
S3 = {4.198062264195162, 5.554958132087372, 7.246979603717467} 
S4 = {4.120614758428182, 5, 6.347296355333861, 7.532088886237956 } 
ie 4.081014052771005, 4.690278532109430, 5.715370323453431, ! 
6.830830026003771, 7.682507065662362 


8 4.058116365147897, 4.502978503657798, 5.290790225914926, 
6= ; 


6.241073360510647, 7.136129493462313, 7.770912051306419 
Then 


> Xi = tr(Opm) = (n — 1)(45? + 2t) + (48? + t) = n(45? + 2t) -t 
[]Ai = det(Opn) = @2n41- 


If we take care of the spectra, minimum eigenvalue converges to 4s7. The 


maximum eigenvalue of spectra converges to 4s? + 4t. 


Theorem 2. /f X; is an eigenvalue of the matrit Opn, then A; + 8s + 4 is an 
eigenvalue of Opn(s +1). 


http://www.earthlinepublishers.com 


On Generating Tridiagonal Matrices of Generalized (s,t)-Pell, ... 239 





Proof. A; is an eigenvalue of Op», then 


|Opn — AiL| 
4(s +1)? +t-(A +8844) t 
t 


t 
t A(s+1)?+2t-— (4; + 85 +4) 


|Opn(s +1) — (A; + 88 + 4)J]. 


1.4 Some properties of tridiagonal matrices Ag, by (s,t)-Pell 
Lucas sequence 


Assume that Ag, is an n x n tridiagonal matrix defined as 


Then the determinant of Ag, 


det Agn = Rn. 
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(s, t)-Pell Lucas sequence is also obtained by using the following symmetric matrix 


with complex entries. Assume that A, is an n x n tridiagonal matrix defined as 


2s 2it 
a 2s at 
1 2s 


Then the determinant of A, is given by (1) as 
det An = Rn. 


The sequence of the spectra of the matrices Ag, for n = 2,3,4,5,6 is given in 
the following: 
So = {2 + 1.414213562373095i, 2 — 1.4142135623730957} 
S3 = {2, 2+ 1.7320508075688777, 2 + 1.732050807568877: } 
J 2+1.8477590650225737, 2 — 1.847759065022573% 
| 2 +.0.765366864730179i, 2 — 0.765366864730179% 
Ss = 2 + 1.9021130325903067, 2 — 1.902113032590306:i 
,2, 2+1.1755705045849467, 2 — 1.1755705045849462 
2 + 1.9318516525781377, 2 — 1.931851652578137%, 
Se = ¢ 2+ 1.4142135623730962, 2 — 1.4142135623730967, 
2 + 0.517638090205041i, 2 — 0.5176380902050412 
Evidently, > A; = tr(Agn) = 28n and [J Ay = det(Agn) = Rn. 
For the inverse of Ag», by using (2), it is obtained that 


4 


a = 25, +>0 

by = 2t, b =F, ai>1 
qq = -l, i>0 

% = 1, &=Ri, i> 
oj = Pn-j+2, J > 0. 
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Therefore the inverse of Ag, is the following matrix 


(-1)' "HR, n—jaa (a) ge, if i<j 
ep), if i=1,j>1 
(AG JGR = a if i=j=1 
epee if j=l,i>l 
Wy—1n—i41 Fe if i>j 


1.5 Some properties of tridiagonal matrices Eg, by even (s,t)-Pell 


Lucas sequence 


Assume that Eg is an Xx n tridiagonal matrix defined as 


As? + 2t oF 
t Ag Oe. f 


_ : 
t ed 


Then the determinant of Eg» is given by (1) 
det Egwn = Ron. 


For the inverse of Eg.n, the values are computed as, 


a = As*+2t, i>0 
b, = 2t,b:=t, i>1 
Gj = t, 1 > 1 
Ao oo t; 6; < Rai, 
§72(n—j+2) : 
eS ENE, Sa 
; 26 ’ i ee 


where 








= hae SR Rei-292(n—j-+1) if acy 
(oi Tait? Raj—202(n—i41) ee 
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The sequence of the spectra of the matrices Eg, for n = 2,3,4,5,6, are given 
in the following 


S = {4.585786437626905, 7.414213562373095} 
S3 = {4.267949192431122, 5.999999999999997, 7.732050807568876} 


_ J 4,152240934977428, 5.234633135269817, 
‘ 7.847759065022574, 6.765366864730179 


ae 4.097886967409691, 4.824429495415054, 6, 
2 7.175570504584945, 7.902113032590311 


4.068148347421865, 4.585786437626904, 
Se = 4 5.482361909794957, 6.517638090205042, 
7.414213562373094, 7.931851652578145 


Then, )> \; = tr(Egn) = n(4s? + 2t) and [] A; = det(Egn) = Ron. 
If we take care of the spectra, minimum eigenvalue converges to 4s?. The 


maximum eigenvalue of spectra converges to 4s? + 4t. 


1.6 Some properties of tridiagonal matrices Og, by odd (s,t)-Pell 
Lucas sequence 


Assume that Og, is an n x n tridiagonal matrix defined as 


then the determinant of Og jn 


det Og,n = Ron-1. 
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For the inverse of Og n, the values are computed by (2) 


ay = 2s, 

a; = As? +21, 7>2 

cy = —2s, bh =t 

Oj = eS ae 

Oo = hee] Ras: eS 
o = PR jen 


If s =¢t = 1, then the sequence of the spectra of the matrices Og, for n = 
2,3,4,5,6 is given in the following 
S2 = {2.585786437626905, 5.414213562373095} 
S3 = {2.657076917222828, 4.529316580128842, 6.813606502648331} 
_ J 2.665585781661023, 4.258036215697276, 
ia 5.741963784302725, 7.334414218338978 
Ss = 2.666546286933180, 4.149626499380563, 5.132659714556408 
6.474094390850837, 7.577073108279015 
2.666653286578264, 4.096987803756087, 
Se = ¢ 4.780438656766319, 5.834296887585696, 
6.913255181694661, 7.708368183618974 
The sequence of maximum eigenvalue is increasing and converges to 4s? + 4t = 
8. So we can say Jim, max(A(O,(k)) = 4s? + 4t. 
Evidently, 3> A; = tr(Ogn) = (n — 1)(48? + 2t) +28 and [] \; = det(Ogn) = 
Mogeas 


1.7 Some properties of tridiagonal matrices A,,, by (s, t)-Modified 
Pell sequence 


Assume that Ag» is an n x n tridiagonal matrix defined as 
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then the determinant of Ay, is 
det Ajx:— Nn: 


(s,t)-modified Pell sequence is also obtained by using the following symmetric 
matrix with complex entries. Assume that A, is an n x n tridiagonal matrix 


defined as 


Then the determinant of A, is given by (1) as 
det Ay, = Nn. 


The sequence of the spectra of the matrices Ag, for n = 2,3,4,5,6 is given in 
the following: 
So = {1.5 + 0.8660254037844387, 1.5 — 0.8660254037844387} 
_ J 1.480159709001947, 1.784920145499028 + 1.307141278682045i, 
= 1.784920145499028 — 1.3071412786820457 
1.895123382259650 + 1.5524918200618802, 
1.895123382259650 — 1.5524918200618802, 
1.604876617740350 + 0.5068439018059832, 
1.604876617740350 — 0.506843901805983: 
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1.941818622615235 + 1.691279149514195z2, 
1.941818622615235 — 1.6912791495141952 
os 1.583716458826265, 
1.766323147971634 + 0.8855567602321662, 
1.766323147971634 — 0.8855567602321662 


1.964532008416013 + 1.7753035706951052, 
1.964532008416013 — 1.775303570695105z2, 
1.678391681959849 + 0.3590790226703332, 
1.678391681959849 — 0.3590790226703332, 
1.857076309624141 + 1.1595156563348332, 
1.857076309624141 — 1.1595156563348332 


S6= 


Evidently, 5° 4 = tr(Agn) = 2sn — s and [] A; = det(Agn) = Xn. 


For the inverse of Ay», by using (2), it is obtained that 


a, = 8, a=2s, 1>1 
b = t, i>0O 

q = -l, i>0 

O = 1, &=%, a>1 
03 = Pn—j+2, J > 0. 


Therefore, the inverse of Aj, is the following matrix 


(—1) 8 1N, 1 nj) if ey 
ah eS ib oh, 1 
(Agni) = x if i=jal 
rae 1 yao 
ee rea eet if +> J 
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1.8 Some properties of tridiagonal matrices E,, by even 
(s, t)-Modified Pell sequence 


Assume that Ey, is an n x n tridiagonal matrix defined as 


287 +t t 
t 4s?42t ¢ 
Eqn = t 
t 
t 4s? + 2¢ 
Then the determinant of Ey, is given by (1) 
det Eqn = Non. 
For the inverse of Eqn, the values are computed by (2) 
a) = 2874+t, aj =4s7+2t, i>1 
b = t,7i>0 
Gj = b 1 > 1 
Ao = 1, 6; = Rai, 
§92(n—j+2) : 
i IS ee > 1. 
’; 2s 2 J a 


Therefore, we get 





E,n)74.3 = Gas Tate? Nei—-22(n—j41) LE eS. 
( ) 2sNon 2j-2§2(n—i+1) Wwe>y 








The sequence of the spectra of the matrices Ey for n = 2,3,4,5,6, are given 
in the following 
Si = {2.697224362268005, 6.302775637731995} 
S3 = {2.669941260432018, 5.201639675723405, 7.128419063844577} 
6 2.667028328506892, 4.700596159221179, 
6.156443643217495, 7.475931869054436 
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2.666706814990114, 4.449448536406551, 

S5 = 4 5.511997050564387, 6.720107682566076, >, 
7.651739915472869 

2.666671127024368, 4.309895427899775, 
Se = § 5.098757535396342, 6.105336258757300, 

7.067037793312944, 7.752301857609272 
Clearly, >> A; = tr(Eqgn) = n(4s? + 2t) — 28? —¢ and [] A; = det(Egn) = Nan. 
If we take care of the spectra, minimum eigenvalue converges to 2s”. The 


maximum eigenvalue of spectra converges to 4s? + 4t. 


1.9 Some properties of tridiagonal matrices O,, by odd 
(s, t)-Modified Pell sequence 


Assume that Og, is an n x n tridiagonal matrix defined as 


ae | 
t-s43 ADE t 


Os = t Age Lig . 


Then the determinant of Ogn 
det OG = Non-1- 


For the inverse of Ogn, the values are computed as 


ay = 8s 

a; = 4s%+2t,i>2 

cq = —s, bb =t 

bb = qg=t, i>l 

Gor = dy BHR, ee 
Gi W294?) ae 


2s 
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If s =¢t = 1, then the sequence of the spectra of the matrices O,,, for n = 
2,3,4,5,6 is given in the following 
Sq = {1.208712152522080, 5.791287847477920} 
S3 = {1.218716204021644, 4.862194798386097, 6.919088997592260} 
1.219199198376559, 4.504733902179790 
5.898455994015420, 7.377610905428227 
1.219222421149124, 4.332385534359584, 
S5 = 4 5.293995934620866, 6.555626057939926, 
7.598770051930496 
1.219223537248886, 4.235965700260784, 
Se = 4 4.930440201174632, 5.932935651316058, 
6.960679940971054, 7.720754969028587 
The sequence of maximum eigenvalue is increasing and converges to 4s*+4t = 
8. So we can say Jim max(A(On(k)) = 4s? + 4¢. 
Evidently, > \; = tr(Ogn) = (n — 1)(48? + 2t) + 8 and J] \; = det(Ogn) = 
Non-1- 


S4= 
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Abstract 


Using the Al-Oboudi type operator, we present and investigate two special 
families of bi-univalent functions in D, an open unit disc, based on ¢(s) = 
ws s > 0, a modified sigmoid activation function (MSAF) and Horadam 
polynomials. We estimate the initial coefficients bounds for functions of 
the type gg(z) = z+ s @(s)d;z) in these families. Continuing the study 
on the initial esdticas of these families, we obtain the functional of 
Fekete-Szego for each of the two families. Furthermore, we present few 
interesting observations of the results investigated. 


1 Preliminaries 


Let the set of complex numbers be denoted by C and the set of normalized regular 


functions in D = {z € C: |z| < 1} that have the power series of the form 
ca . 
g(z) = 2+ doz” + dgz* +... =2+ So dj2, (1.1) 
j=2 


be indicated by A and the set of all functions of A that are univalent in D is 


symbolized by S. The famous Koebe theorem (see ) ensures that any function 
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" satisfying z = g~'(g(z)), w = g(g7*(w)), lw] < ro(g) 


and ro(g) > 1/4, z,w € D, where 


g € S has an inverse g~ 








g + (w) = fw) =w — dow" + (242 — d3)w® — (5d8 — Sded3 + da)w*+.... (1.2) 


A function g of A is said to be bi-univalent (or bi-schlicht) in D if g and its 
' are both univalent (or schlicht) in. The set of bi-univalent functions 
having the form is indicated by }>. Historically investigations of the family 
)> begun five decades ago by Lewin and Brannan et al. [9]. After few years, 


Tan [40] found the initial coefficient bounds of bi-univalent functions. Later, 


inverse g~ 


Brannan and Taha presented and investigated certain subsets of > similar 
to convex and starlike functions of order o (0 < o < 1) in 9. Some interesting 
results concerning initial bounds for certain special sets of }> have been appreared 


in [11], [18] and [32]. 


Let the set of real numbers be R = (—oo, 00) and the set positive integers be 
Ne=Ng\{0}-=41, 2, 3, ...1. 

Recently, Hérzum and Koger [21] (see also [20]) examined the Horadam 
polynomials H;(x) (or H;(x,a,b;p,q) ). It is given by the recurrence relation 


Ag) = pallyeal@) + qijo(@), Bye) =a, Hox) = be, (1.3) 


where j € N\{1,2}, x © R, p, g, a and b are real constants. It is easy to see from 
(1.3) that H3(x) = pbx? + qa. The generating function of the sequence H;(z), 
j EN, is as below (see [21]): 


_ at (b— ap)xz 


1.4 
1 — prz — qz?’ Le) 


where z € C is independent of the argument x € R, that is R(z) # z. 


Few particular cases of H;(x, a,b; p,q) are: 
i) Hela, dy el 1) = Fe), ti) (a1, 252,=1) = Ua); 
it) te Ne 2 A) ee), 1) Bate 2,41) = ae), 
0) e222, =O) amd, ot) Aa (a, 1.22, 1) = Pe): 
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They are named as Fibonacci polynomials, second type Chebyshev polynomials, 
first type Chebyshev polynomials, Lucas polynomials, Pell-Lucas polynomials and 


Pell polynomials, respectively. 


In the literature, the estimates on |d2|, |d3| and the famous inequality 
of Fekete-Szeg6 were determined for bi-univalent functions related to certain 
polynomials like Fibonacci polynomials, (p,q)-Lucas polynomials, second kind 
Chebyshev polynomials and Horadam polynomials. We also note that the above 
polynomials and other special polynomials are potentially important in statistical, 
physical, mathematical and engineering sciences. Additional information about 
these polynomials can be found in [7], [8], [16], [17], and [42]. More details 


about the famous Fekete-Szeg6 problem connected with Haradam polynomials are 


available with the works of [I], [2], [3], [26], [31], [88] and [4]]. 

The recent research trend is the study of bi-univalent functions linked with 
any one of the above mentioned polynomials using well-known operators, which 
can be seen in the research papers [4], [13], [25], [28], [34], [36], and [89]. 
Generally interest was shown to estimate the initial Taylor-Maclaurin coefficients 
and the celebrated inequality of Fekete-Szego for the special families of 5° that 
are being introduced using known operators. 

In this work, we present two special sets of }> using Al-Oboudi type operator 
which was precisely defined in the paper [19]. We determine the initial coefficient 
bounds and also obtain the relevant connection to the celebrated Fekete-Szeg6 


functional for functions in the defined families. 


Let Ay denote the set of regular functions of the form 
Pe. es 
go(z) = z+ > b(s)d;z7, 
j=2 


where ¢(s) = eo s > 0, is a MSAF. Note that ¢(0) = 1 and hence A; := A 
(see [14]). 


Definition 1.1. For gg € Ag, k € No, 8 => 0, an Al-Oboudi type operator Ds : 
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Ag > Ag, is defined by 
D9(z) = ge(z), Dage(z) = (1 — B)g6(z) + 8294 (2), ---, Dgge(z) = Da(D§*9(z)), 
zeED. 


Remark 1.1. If gg(z) = 2+ > ¢(s)djzI € Ag, z € D, then 
j=2 


D5g9o(2) =z+ ya + (j — 1)8)*6(s)djz*, z € D. 
j=2 


When ¢(s) = 1, we get the Al-Oboudi operator [5], which reduces to the 
Salagean operator [29], if G=1. 

For regular functions g and f in 9, g is said to subordinate to f, if there is a 
Schwarz function ~ in D, such that ~(0) = 0, |¢(z)| < 1 and g(z) = f(¥(z)), z € 
®. This subordination is indicated as g < f or g(z) ~ f(z). Specifically, when 
f €S in 9, then g(z) ~ f(z) is equivalent to g(0) = f(0) and g(D) C f(D). 

Inspired by the articles [6], and the trends on functions € 5° , we present 
two special families of > by using Al-Oboudi type operator, which is as in 
Definition 1.1 and Horadam polynomials H;(x) as in the relation having 
the generating function (1.4). 


Throughout this paper, fg(w) = Gy (w) is an extension of g~! to D given by 


(1.2), p, g, a and b are as in (1.3) and G is as in (1.4), unless and otherwise 


mentioned. 
Definition 1.2. A function g in > having the power series (1.1) is said to be in 
the family SGs-(z,7,u,k, 6,0(s)), OS y < 1,420, KE No, B = 0 and ¢(s) 
the MSAF, if 
z(Dgge(z))! + wz?(Dhgg(z))” 
yDhge(z) + (1—y)z 





a+G(xz,z),z€D 


and 
w(Di fg(w))’ + wa? (Di fg(w))” 


Di fs(w) + (1 — 9) 





<~1l-—a+G(z,w),w €D. 
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We note that i) w = 0, ii) y = O and iii) y = 1 lead the family 
SGs-(z,7,H,k, 8, @(s)) to the below mentioned subfamilies: 

1. SKy(a, 7, k, 6, ¢(s)) = SGy(a, 7, 0,k, 6, 0(s)) is the set of functions g 
€ > satisfying 


2(Digo(2))' 
7D590(2) +(1-y)z 


w(D5 few)! 


WDhfyu) + aya 1 ee) 





~< 1—a+G(a,z), and 





where z, w € D. 
2. SLyx(x, bw, k, B,0(s)) = SGy-(z, 0, uk, 8, O(s)) is the family of functions 
g € >> satisfying 


(D§9o(2))' + uz(Dbgg(z))” < 1- a+ G(a,z) 


and 
(D5 fs(w))’ + pw(D§ fo(w))” <1—atG(a,w), 
where z,w € D. 
3. SMs(x, 1,k, 8, O(s)) = SGs-(z,1,H,k, 8, O(s)) is the family of functions 
g € > satisfying 
(Se) : (Ses 
Di9o(2) ‘ 


ee) s (“a 
Di few) J” 
where z, w € D. 

Letting k = 0 and ¢(s) = 1 in the Definition 1.2, we obtain the family 
SNx (2, 7; H) = SGy-(z,7, u,0, 8,1) of functions g € )) satisfying 


2g'(z) + w279"(2z) wf" (w) + po? fw) 
yg(z) + (1 — yz Vf(w) + (1 y)w 
where z, w € D, f(w) = g-!(w) is as given by (1.2), a is as in (1.3) and G is as in 
. 


=< a+G(x,w), 





<1-—a+G(a,z) and 
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Definition 1.3. A function g € )> having the power series is said to be in 
the family SB5-(z, 7,7, k, 8, ¢(s)), OS y<1,721,k © No, 8 = 0 and ¢(s) the 
MSAF, if 

2{(D§9¢(z))')" 


<~1l-—a+GQ(z,z),z€D 
yDh9e(z) + (1 —y)z 





and 
w[(D5fo(w))']” 
Dh fole) + 0 =) 
Note that the certain choices of 7 lead the family SBs>(z,7,7,k, 8,(s)) to 
the following two subclasses: 
1. SP3(z,7,k, 8, 6(s)) = SBs>(z,0,7,k, 8,0(s)) is the set of functions g 
€ > satisfying 





<~1l-—a+G(z,w),we®D. 


[(D§90(2))'” < 1-a+G(az,z),z€D and (Di fe(w))']” < 1-a+G(r,w),w € D, 


2. SNy-(z,7,k, 8, O(s)) = SB (az, 1,7,k, 6, 0(s)) is the class of functions g 
€ > satisfying 
2{(D§go(z))'" w[(D5 fo) 

D590(2) D5 fo(w) 

SMs-(x,7,k, B,0(s)) is the family of Al-Oboudi type r-bi-pseudo-starlike 
functions associated with Horadam polynomials involving the MSAF. 

On taking k = 0 and ¢(s) = 1 in Definition 1.3, we get the family 
SQx(2,7,T) = SBy (2,7, 7,0, 8,1) of functions g(z) in )) satisfying 

w(f"(w))" 
yf (w) + (L- yw 
where z, w € D, f(w) =g7!(w) is as given by (1.2), a is as in (1.3) and G is as in 
(1). 








< 1-a+G(az,z), z €D and < 1-a+G(r,w),w €D, 


<1l—a+G(a,z) and <1l=a+ G(a,w), 








Remark 1.2. We note that i) SB5>(z,7,1,k, 8,6(s)) = SKs(z, 7, k, 8, 0(s)), 
ii) SMy(a, 1, k, B, o(s)) = SKy(z, 1, k, B, $(s)) = SMy (za, 0, k, 8, 6(s)) and 
iii) SP5(a, 1, k, 8, o(s)) = SKy(z, 0, k, 8, O(s)) = SLy-(z, 0, k, 8, 0(s)). 
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Remark 1.3. i) For u = 7 = 0, the class SNy(az,0,0) = H»(x) was studied by 
Alamoush [2] and ii) For » = 0 and y = 1, the family SNy(z,1,0) = S$(x) was 
investigated by Srivastava et al. [32]. 


Remark 1.4. i) For y = 0, the family SQy(x,0,7) = S$(a#,7) was investigated 
by Abirami et al. and ii) For 6 = 1, the family SMy:(z,7,k, 1,¢(s)) = 
My(a, T, k, &(s)) was considered in [25]. 





Remark 1.5. In a special situation, if we choose a 1lb=p=2,¢ —1 
and x — t, the generating function (1.4) reduces to the second type Chebyshev 
polynomials U;(t), which is explicitly given by 





(t= sw) 





3. i) _ sin(g + 1)u 


(t) = (7 +1) oF | -7,74+2; 
U;(t) (Gj+1) 2 1 (5.5425; 9 simap 


in terms of the Gauss hypergeometric function .2F;. In this particular 
situation, the bi-univalent function families SGyx(x,7,u,k, 8,¢(s)) and 
SB (x, 7, 7,k, B,0(s)) would become the families SGy>(t, y, u,k, B,0(s)) and 
SBs(t, 7, 7,k, B,0(s)), respectively. The families SB5(t,1,7,0, 8, ¢(s)) = 
AOs(t,T, O(s)) and SBys(t, 1, 7, 0, 6,1) = AYS-(t,7) were studied earlier in 
and [7], respectively. 


In Section 2, we derive the estimates for |dz|, |d3| and the inequality of Fekete- 
Szego for functions of the form € SGy\(z,7,u,k, B,o(s)) and we also 
present some observations of our result. In Section 3, we derive the estimates 
for |d2|, |d3| and the Fekete-Szegé inequality for functions of the form (1.1) 
€ SBy-(2,7,7,k, 8, ¢(s)). Few interesting consequences of the result are also 


considered. 


2 Estimates for Function Family SGs(z,y, u,k, 8, 0(s)) 


In the following theorem, we determine the initial coefficients bounds and the 


inequality of Szeg6 for functions in SGy>(z, 7, u,k, 8, (s))- 
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Theorem 2.1. LetO0<y<1,u~>0,k © No, 8 > 0 and G(s) be the MSAF. If the 
function g € SGs(x,7,u,k, 8, 0(s)), then 















































\ds| < |ba|./|ba| 
~ (1+ B)¥b(8)/](y? — Qu + 3)y + 3(Qu + 1) (62)? — (H+ 1) — 7)? (pba? + es 
2.1 
1 (bx)? | ba: | 
\d3| < + (2.2) 
(1+ 28)*d(s) [2(ut+1)—7)? ((@u+1)-7) 
and ford ER 
|d3 — 5d3| 
|bar| . (1+28)*65 
(1428) (8) BQu+1)=7) i |1— tages | S 4 
< bal? [1 —a+28"6 
[bal +8)? o(5) al (1428)"5 |x 7 
(14+28)* @(s)|(q7-(244+3)7+3(2u-+1)) (bx)? —(2(u+1)—y)?(pba2+qa)| ? (1+8)?*d(s) | — “? 
(2.3) 
where 
1 2 2 pbx? + ga 
= 2u+3 + 3(2u+1 2 me 
Se peay 7 Oe + 87 + 32u +1) - C+D) 9)? (PAE 
(2.4) 


Proof. Let g © SGy(z,7,u,k, B,0(s)). Then, for two regular functions I, N 
with IN(0) = 0, |Nt(z)| < 1, 9t(0) = 0 and |MN(w)| < 1, z,w € D and on account of 
Definition [1.2] we can write 


z(Dhgg(z))! + w2*(D¥ge(z))” 
yDhg9g(z) + (1-2 





=1—a+G(x,M(z)) 


and 
w(D5fs(w))! + pw?(Dé fo(w))” 
DR foe) + = 7) 





Or, equivalently 


z(Dhg¢(z))' + wz?(Dgge(z))” 
yDhgo(z) + (1 —y)2 





=1-—a+4+ Ay(x) + Ho(x)m(z) + H3(x)(m(z))? +... 


http://www. earthlinepublishers.com 


Coefficient Bounds for Al-Oboudi Type Bi-univalent Functions ... 259 





and 

w(D5 fo(w))’ + uw? (D5 fo(w))” 
Dif) += 7) 

From (2.5) and (2.6), in view of (1.3), we find 


z(Dgge(z))! + wz?(Digg(z))” 
yDhgo(z) + (1— yz 





= 1-a+H,(x) + Ho(x)n(w) + H3(x)(n(w))* 4+... 


(2.6) 





= 14 Ho(x)myz+([Ho(x)m2+H3(x)m?]z7+... (2.7) 


and 


w(Di fs(w))! + me? (D5 fo(w))" 


yD fg(w) + (1 - yw = 1+ Ho(x)myw + [Ho(x)ny + H3(x)nq]w? + .... 





(2.8) 

It is well known that if |9(z)| = |miz + mgz” + m3z2+...| < 1,2 € D and 
|Nt(w)| = |nyw + Now? + ngw® + | <1,w€9%, then 

|m;| < Land |nj| < 1(¢ € N). (2.9) 


We easily get the following by equating the corresponding coefficients in (2.7) and 


(2.8): 





(1+ 8)*6(s) (2(u + 1) — y)dz = Ho(x)my (2.10) 
(1+28)*6(s)(3(2u+1)—y)d3—(1+8)?*6?(s)(2(ut+1)—7)7 dj = Ho(x)m2+H3(x)m} 
(2.11) 

— (1+ 8)*4(s) (2(u+ 1) — 7) dp = H(z) (2.12) 


—(1 + 28)*6(s)(3(2u + 1) — yds + (1 + B)?*G?(s)(7? — 2(w + 2)y + 6(2pu + 1)) 43 
= Ho(ax)ng + H3(a)n?. 


(2:13) 
From and (2.12), we easily obtain 
my, = —ny (2.14) 
and also 
2(1 + 8)°*$?(s)(2(u + 1) — 7)*d3 = (my + nf) (Ho(2))?. (2.15) 
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If we add ( and (2.13), then we obtain 


2(1+ B)**6"(s)(7? — (2Qu+3)7+3(2n+ 1))d3 = Ho(x) (my +g) + H3(x)(mj +07). 
(2.16) 
Substituting the value of m? +n? from (2.15) in (2.16), we get 


(H(a))? (mg + ng) 





& = 
2 2(1 + B)?*o?(s) [(y? — (24 + 3)7 + 3(2u + 1))(ho(a))? — (2(u + 1) - a 
which yields on using (2.9). 
After subtracting (2.13) from (2.11) and then using (2.14), we obtain 
2k _ 








(1 +3 2(1 + 28)*$(s)(3(2u + 1) — ¥) 
Then in view of (2.15), becomes 
(H(a))? (mj + nj) H(x)(mz — ng) 
2(1 + 26)¥6(s)(2(u+1)— 7)? © 2(1 + 28)*d(s)(3(2u + 1) — 7)’ 
which yields on using (2.9). 
From and (2.18), for 6 € R, we get 


lds — 63] = |Ha(2)| (76.2) | 





dz = 


1 
2(1 + 28)*b(s)(3(2u + 1) - a) ia 








r} 





1 
(nr ”) — 30+ 9p) o(s)(3Qn + 1) — 5) = 
where 


T (6, x) 


_ (1+28)* 
2(1 + B)?*¢?(s) [(7? — (Qu +. 3)7 + 3(2u-4 
In view of (1.3), we conclude that 





Casa 8) — 5) (Hp 
iG . — (2(u + 1) — 7)? Hs(x)] 














|H2(x)| 3 1 

ds — 63) < J ROCA) §°SIFOM S cE OBE 
IHANIT6.2)1 5 INGE scoarectsoys 

which gets (2.3) with J as in (2.4). This evidently ends the proof of Theorem 
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By setting i) uw =0, ii) y = 0, iii) y = 1 and iv) k = 0, ¢(s) = 1 in Theorem 


2.11 we have the following four corollaries, respectively. 


Corollary 2.1. If the function g € SKyx(a, 7, k, 8, ¢(s)), then 


|ba|./|bx| 






































|d2| < ’ 
(1+ B)* G(s) /|(7? — 3y + 3)(bx)? — (2 — 7)*(pba? + ga)| 
\ds| < 1 | bee? ES et] 
*~ 1 +28)F6(s) [@—7? * 3-4 
and for some 6 ER, 
[bar ji - (14+26)*6 
5 (1+28)* $(s)(3—7) ; (1+8)?*6(s) | = 
|ds ~~ dd5| < |ba|3 p28) 5_ 
(1+8)?* 6(s) ; |1- (1+22)*6 a: 
(1+28)*9(s)|(9?-37+3) (bx)? —(2—7)? (pba? +qa)| ero) 
pba? + ga 
where J, = a5 q7 —-3y+3-(2 vw 5a ) 








Remark 2.1. For y = 6 = 1, Corollary [2.1] reduce to Corollary 2.1 of Magesh et 
al. and Corollary [2.1] further coincide with Corollary 2.1 of Abirami et al. [1], 
when k = 0 and ¢(s) = 1. Corollary [2.1] coincide with Theorem 2.2 of Alamoush 
[3], when y = k = 0 and ¢(s) = 1 and also we obtain Corollary 1 and Corollary 3 


of for k= 0,4 = ole) =I. 
Corollary 2.2. If the function g € SLs(x, y, k, 6, 0(s)), then 


|ba|./|bx| 


(1+ B)*b(s) V13(2n + 1) (bx)? — 4(u + 1)? (pba? + qa) 


1 D2 
aaa | = £ |ba| | 











|do| < 
































(1+ 28)¥d(s) [4(u +1)? 3(2u+1) 
and ford ER, 
|ba| (1+2)*65 

3. | 3@eD +28) ACs) : | abe] S 2 

lds — dd3| < ibn) [1—-G+200"5_ 
(1+8)?* 6(s)  |l- 7 (14+26)F5_ | s J 
(1+28)¥ 6(s)[32u+1) (bx)? —4(u-+1)? (pba?-+Ga)| 1+f)** G(s) | — °?? 
2 

_ 4(u+1)2 ( pbx* + ga 

where Jo = 1 BOD) ( 52? )} 
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Remark 2.2. For 4. = k = 0 and ¢(s) = 1 Corollary [2.2] coincide with Theorem 
2.2 of Bl. 


Corollary 2.3. If the function g © SMs-(x, u,k, B,0(s)), then 





























Id) < |ba|./|bax| 
(1+ B)F4(s)4/ [](4u + 1) (bx)? — (24+ 1)?(pbx? + ga)]] 
1 b2 2? |bx| 
|d3| < ki 2 
(1+ 26)Fo(s) | (2u4+1)? © 2(34 +1) 
and ford ER, 
br . |, G428)6 
4». | ARTA STS) ; | caaEses| S 4 
|d3 — dd5| < \bo|3 [1 @-+2ayks 
(1+8)?* 6(s) : r (14+28)*6 a: 
(FBO Gut 1) (Ox)?—Ont iba aay) (+P @s)| = > 











2 a 
where Js = agaq (40 Pl) = (Qu 1)? (Pirtt* ) 
Remark 2.3. Corollary coincide with Theorem 2.1 of Magesh et al. [26], 
when k = 0 and ¢(s) = 1. Also we obtain Corollary 2.1 of [25] from Corollary 
[2.3] when pp = 0 and 6 = 1. 

Corollary 2.4. If the function g(z) € SNsx(x, y, 1), then 


























\dy) < |ba:|./|bax| 
~ V\(q? = (Qu + 8)y + 3(2u + 1)) (bx)? — (2(u + 1) — 7)? (pba? + qa)| 
2 
b 
Ids) < (bz) a |ba| 
(Qu+1)—y)? (824+ 1)—-7) 
and ford ER 
|ba| F 
|d3 — 5d2| < 3(Qu+1)—7 j|L— 9d] S Ja 
et \ba-|® |1—3| lls > J 
[72 (23) 7+ 3H) (be)? (AL) —y)2 (pba Gayl? ali 
where 
_ 1 2 2 ( pba? + qa 
Ih a |? — Cu + 8)7 + 82u4+ 1) — 4 1) 9) (AES 





Remark 2.4. By choosing appropriate values for parameters y and yp in Corollary 
we obtain Theorem 2.2, Theorem 2.1 and Corollaries 1, 2 of [3], and [31], 


respectively, as it can be seen from earlier remarks. 
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3. Estimates for the Function Family 
DBS (x, VT; k, B, (s)) 


In the next theorem, we find the first two Taylor-Maclaurin coefficients and the 


inequality of Fekete-Szeg6 for functions in SBy-(2,7,7,k, 6, 0(s)). 


Theorem 3.1. LetO0<y<1,7>1,k € No, 8 > 0 and G(s) the MSAF. If the 
function g € SBs(z, 7, T,k, 8, (s)), then 


\dy| < |ba:|,/|bax| 


(1+ B)FG(s) IV? + (7 — 7) (27 + 1) (bx)? — Or — 7)?(pba? + ga)]’ 




















(3.1) 
1 bax)? ba 
|ds| < ; | Oz) 5 4 [ba (3.2) 
(1+ 28)*9(s) [27—-)? (87-7) 
and ford ER, 
|b(x)| Ace (14+26)"6 [<0 
1428)" d(s)(37— , 1+B)?" ¢(s)! = 
Ids — 5d2| < (14+:28)¥ 6(s) (37-7 se ituaasi gai (1+8)?* G(s) 
(1+)? 4(s) | (14+28)*6 |>a 





(1+28)* 6(s)|(9?+(7-7) (27 +1)) (bx)? — (27)? (pbx? +qa)| (1+8)?* 6(s) 


(3.3) 


where 








1 pbx? + qa 
Qe 24 Qr+1 2 a 

Bag | + = ver +) - era» (PES 
Proof. Let g © SBsx(x,7, T,k, 8, (s)). Then, for some regular functions IN and 
MN such that 9t(0) = 0, |M(z)| = |myz + mez? + mgz3 4 ...] < 1,9(0) = 0 and 
|St(w)| = |nyw + now? + ngw3 +... <1, z,w € D and on account of Definition[1.3} 


we can write 





2{(Dige(2))'" 
yDh9g(z) + (1— yz 





=1—a+G(z,M(z)),zED 
and ; 
w[(Defo(w))']” 
yD) f(w) + (1 — y)w 
Following the procedure similar to the proof of Theorem [2.1] one gets 





=1-—a+G(z,N(w)),w ED. 


(1+ 8)*(2r — 7)$(s)d2 = Ho(x)my (3.4) 
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(1+8)°*¢?(s)(?—-277y+2r (7-1) )d3+(14+-28)*6(s)(37-7)d3 = Ho(x)m2+H3(x)mj 
(3.5) 

— (1+ )*(2r — y)6(s)dz = Ho(x)m (3.6) 

(1 + B)?*6?(s)(7? — 2(7 + ly + 27(7 + 2))d3 — (1 + 28)"6(s)(37 — yds (3.7) 


= Ho(x)ng + H3(ax)n3 


The results (3.1)-(3.3) now follow from (3.4)-(3.7) by adopting the procedure as 
in Theorem 














By setting i) y = 0, ii) y = 1 and iii) k = 0, ¢(s) = 1 in Theorem |3.1| we 


have the following three corollaries. 


Corollary 3.1. /f the function g € SP5x(a,7,k, 8, (s)), then 

















Ids) < |ba|./|bax| 
(1+ B)F6(s) Ir @r + Dba)? — 472 (pba? + ga) 
1 (bx)? |ba| 
|d3| < 7 A 
(1+ 28)*@(s) | 4r oF 
and ford ER, 
[b(sx)| . [zy — G+28)*6 
> | RAKES) ; [1 atime] S & 
lds — dd3| < 1 0+28)*5 |i 3 
(1+8)2* 6(s) ; I (1+26)*5 | Q 
(1+26)*¢(s)|r(27+1)(bx)?—47? (pba?+qa)| ’ C+6)?* 6(s) | — °°? 











where Qy = } \(2r + 1) —4r (erst) . 





Remark 3.1. Corollary [3.1] coincides with Theorem 2.1 of [3], when k = 0 and 
T=]0\s) =. 


Corollary 3.2. If the function g © SNyx(a,7,k, 8, O(s)), then 











\da| < |ba| y/|ba 
(1 + B)¥o(s) /|(r (27 — 1) (bx)? — (27 — 1)? (pba? + qa)| 
il (bx)? _ [ba 
3! = Gy 2Bya(s) [Or — I? * Gr-1) 
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and ford ER, 
b(x)| . — _(1+28)*6 
; BF 6(s)(r=1) ; | ~ aeaeaes| S % 
|d3 — 6d5| < 1 4426)5_| pos 
(148)? 6(s) ; 1 (14+26)*6 + QO 
(1+26)* ¢(s)|(7(27—-1)) (bx)? —(27—-1)? (pba? +qa)|? (1+6)?* 6(s) | — °"?? 








where Qz = ea (rer 1)) — (27 — 1)? (Bers ) 





Remark 3.2. Corollary [3.2] reduces to Theorem [2.1] of [25], when 6 = 1 and also 
the results of Corollary [3.2] coincide with Theorem 2.1 of Abirami et al. [I], when 
k =0 and ¢(s) =1. 


Corollary 3.3. If the function g © SQsx(x,7,7), then 


Ids) < |ba:|./|bax| 


~ VA (9? + (r — 9) (27 + 1))(bx)? — (27 — 7)? (pba? + ga)| 


(ba)? [ba 








o] 




















d3| < 
lS Gap Gra) 
and ford ER, 
|b(x)| -|1—6|< Q 
lds — dd] < 9 7 aly 
|1—6|]b2| -|1-6/ > 
I(y?2+(7-7) (27 +1)) (bx)? (27-7)? (pba? +qa)|? oe 
where 
1 pbx? + qa 
2= geal eer 43) - Gra} (Pe) 





Remark 3.3. Corollary [3.3] reduces to Theorem [2.1] of [1], when y = 1. 


4 Conclusion 


Two special families of holomorphic and bi-univalent (or bi-schlicht) functions 
are introduced by using Al-Oboudi type operator involving a modified sigmoid 


activation function associated with Horadam polynomials. Bounds of the first 
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two coefficients |d2|, |d3| and the celebrated Fekete-Szeg6 functional have been 
fixed for each of the two families. Through corollaries of our main results, we 
have highlighted many interesting new consequences. 

The special families examined in this research paper using AlOboudi 
type operator could inspire further research related to other aspects such as 
families using q-derivative operator [22], [85], meromorphic bi-univalent function 
families associated with Al-Oboudi differential operator and families using 
integro-differential operators [27]. 
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Abstract 


In this paper, we introduce and clarify a new presentation between the 
n-exact sequence and the n-injective module and n-projective module. Also, 
we obtain some new results about them. 


1 Introduction 


Category theory formalizes mathematical structures and their concepts in terms 
of a labeled directed graph called a category, whose nodes are called objects, and 
their edges called arrows (or morphisms). This category has two basic properties: 
the ability to compose the arrows associatively and the existence of an identity 
arrow for each object. The language of category theory has been employed to 
formalize concepts of other high-level abstractions such as sets, rings, and groups. 
Several terms were utilized in category theory, including the 4morphisma that is 
used differently from their usage in the rest of mathematics. In category theory, 
morphisms obey specific conditions of theory. Samuel Eilenberg and Saunders Mac 


Lane introduced the concepts of categories, functors, and natural transformations 
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in 1942-45 in their study of algebraic topology, to understand the processes that 
preserve the mathematical structure. Category theory has practical applications 
in programming language theory, for example, the usage of monads in functional 
programming. It may also be used as an axiomatic foundation for mathematics, as 
an alternative to set theory and other proposed foundations. In mathematics, an 
abelian category is a category in which morphisms and objects can be added 
and in which kernels and cokernels exist and have desirable properties. The 
motivating prototype example of an abelian category is the category of abelian 
groups, Ab. The theory originated to unify several cohomology theories by 
Alexander Grothendieck and independently in the slightly earlier work of David 
Buchsbaum. Abelian categories are very stable categories. For example, they 
are regular and satisfy the snake lemma. The class of Abelian categories is 
closed under several categorical constructions, for instance, the category of chain 
complexes of an Abelian category, or the category of functors from a small category 
to an Abelian category are Abelian as well. These stability properties make 
them inevitable in homological algebra and beyond. This theory has significant 
applications in algebraic geometry, cohomology, and pure category theory. The 
Abelian categories are named after Niels Henrik Abel. An exact sequence is a 
concept in mathematics, especially in group theory, ring, and module theory, 
homological algebra, as well as in differential geometry. An exact sequence is a 
sequence, either finite or infinite, of objects and morphisms between them such 
that the image of one morphism equals the kernel of the next. Homological algebra 
is the branch of mathematics that studies homology in a general algebraic setting. 
It is a relatively young discipline, whose origins can be traced to investigations in 
combinatorial topology (a precursor to algebraic topology) and abstract algebra 
(theory of modules and syzygies) at the end of the 19th century, chiefly by Henri 
Poincaré and David Hilbert. The development of homological algebra has closely 
intertwined with the emergence of category theory. By and large, homological 
algebra is the study of homological functors and the intricate algebraic structures 
that they entail. One quite useful and ubiquitous concept in mathematics is 
that of chain complexes, which can be studied both through their homology 


and cohomology. Homological algebra affords the means to extract information 
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contained in these complexes and present it in the form of homological invariants 
of rings, modules, topological spaces, and other atangiblea mathematical objects. 
A powerful tool for doing this is provided by spectral sequences. From its very 
origins, homological algebra has played an enormous role in algebraic topology. Its 
sphere of influence has gradually expanded and presently includes commutative 
algebra, algebraic geometry, algebraic number theory, representation theory, 
mathematical physics, operator algebras, complex analysis, and the theory of 
partial differential equations. K-theory is an independent discipline that draws 
upon methods of homological algebra, as does the noncommutative geometry of 
Alain Connes. This paper is organized as follows. 

In this paper, we show to prove the important theorems of n-injective modules 
and n-projective modules. Finally, we recall the definition of n-projective module, 


and we give an open problem about some theorems of n-projective modules. 


2 Preliminaries 


All rings R in this paper are assumed to have an identity element 1 (or unit) 
(where rl = r = Ir for all r € R). We do not insist that 1 4 0; however, should 
1 =0, then R is the zero ring having only one element. 

In this section, we recall some of the fundamental concepts and definitions, 


which are necessary for this paper. For details, we refer to [4,6,7,9,10,11]. 


Definition 2.1. An R-module M is injective provided that for every 
R-monomorphism g : A —> B between R-modules, any R-homomorphism 
f : A—  M can be extended to an R-homomorphism h : B —> M such that 


hg = f; 1.e., the following diagram commutes 


g 


0 ——> A —*> B 


les 
M 


7 


Definition 2.2. An R-module P is projective provided that for every 
R-epimorphism g : A —> B between R-modules and R-homomorphism f : 
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P —+» B, there exists an R-homomorphism f : P —> B, there exists an 
R-homomorphism h : P —+ A such that gh = f; ie., the following diagram 


commutes 


P 
ra 
A a B= 0 


Definition 2.3. A left R-module F is a free left R-module if F' is isomorphic to 
a direct sum of copies of R; that is, there is a (possibly infinite) index set B with 
F = ®sepRp, where Ry = (b) = R for all b€ B. We call B a basis of F. 


Definition 2.4. Let M be an R-module. An element m € M is divisible provided 
that for any r € R that is not a righ zero-divisor, there exists an x € M such that 
m = ra. We also say that M is a divisible module module provided that every 


element of M is divisible. Note that a divisible group is a divisible Z-module. 


Definition 2.5. Let C be an additive category and f : A — B a morphism in 
C. A weak cokernel of f is a morphism g : B —+ C such that for all C’ € C the 


sequence of abelian groups 
COCNZA CBO) 564.0) 
Definition 2.6. A category C is abelian if 


1. C has a zero object. 

2. For every pair of objects there is a product and a sum. 
3. C Every map has a kernel and cokernel. 

4. C Every monomorphism is a kernel of a map. 


5. C Every epimorphism is a cokernel of a map. 
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Definition 2.7. A category C is additive if 


1. Hom(A, B) is an (additive) abelian group for every A, B € obj(C) 


2. the distributive laws hold: given morphisms 


f 


Gat fey Bee yy 





and 


ty g AR ey 





where X and Y € obj(C), then 


b(f +9) = bf + bg 


and a + [underloverb 
(f+g)a= fa+ga 


3. C has a zero object. 
4. C has finite product and finite coproduct. 


Definition 2.8. An abelian group D is said to be divisible if given any y € D 
and 0 #n € Z, there exists x € D such that nx = y. 


Example 2.9. 


1. Note that Q is a divisible Z-module since for every q € Q, where q = § for 
integers a,b © Z with b £0, and for every 0 4 z € Z, there exists x € Q 


such that « = & so that q = za. 


2. Note that Z is not a divisible Z-module since there is no x € Z with 3 = 2a. 
Definition 2.10. Let Mjjcz be a family of R-modules, and let fi;<z be a family 
of R-homomorphisms such that M;_1 dts Mi for every i € Z. Then the sequence 


eS) ea hee ee ee (2.1) 
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is said to be exact provided that Im(fi — 1) = Ker(f;) for every i € Z. Note that 





(=A BA eG 0 (2.2) 
is exact if and only if f is an R-monomorphism, g is an R-epimorphism, and 
Im(f) = Ker(g) 
This type of sequence is called short exact. 


Definition 2.11. Let C be an additive category and f : A —> B a morphism in 
C. A weak cokernel of f is a morphism g : B —+ C such that for all C’ € C the 


sequence of abelian groups 
cic, c’) 45 e(B,c") 4 c(A, Cc") 


is exact. Equivalently, g is a weak cokernel of f if fg = 0 and for each morphism 
h: B — C’ such that fh = 0 there exists a (not necessarily unique) morphism 
p:C —+C’ such that h = gp. These properties are subsumed in the following 


commutative diagram: 


Clearly, a weak cokernel g of f is a cokernel of f if and only if g is an epimorphism. 


The concept of weak kernel is defined dually. 


Definition 2.12. A morphism f : A —> B in C is called ¥-monic if 


) 


c(B, x) “4 e(a, x) —0 


is exact for any object X € XY. A morphism f : A —> X in C is called a 
left V-approximation of A if f is ¥-monic and X € XY. The subcategory ¥ is 
said to be a covariantly finite subcategory of C if any object A of C has a left 
X-approximation. We can defined -epic morphism, right ¥- approximation and 
contravariantly finite subcategory dually. The subcategory ¥ is called functorially 


finite if it is both contravariantly finite and covariantly finite. 
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Definition 2.13. Let C be an additive category and d° : X° —> X! a morphism 


in C. An n-coker of d° is a sequence 





1 2, nm 
(discal a a ee Ea 
such that, for all Y € C the induced sequence of abelian groups 
n+l dr >, qn-1 di 1 qo 0 
OS C(X"™ VY) Se C(X", VY) SS... C(x, Y) C(x”, Y) 


is exact. Equivalently, the sequence (d',...,d”) is an n-coker of d° if, for all 
1<k<n-—1 the morphism d’ is a weak cokernel of d*~!, and d” is moreover a 


cokernel of d”~!. In this case, we say the sequence 


gOS Sa Ti A ye a aan 





is right n-exact. 


Remark 2.14. When we say n-cokernel we always means that n is a positive 
integer. We note that the notion of 1-cokernel is the same as cokernel. we can 


define n- kernel and left n-exact sequence dually. 


Definition 2.15. Let C be an additive category. An n-exact sequence in C is a 


complex 


a) a: m—1 nm 
ROSS ae, eer (2.3) 





in Ch"(C) such that (d°, ...,d"~') is an n-ker of d”, and (d', ...,d”) is an n-coker of 


d°, The sequence (3.1) is called n-exact if it is both right n-exact and left n-exact. 


Theorem 2.16. Let A,B, {Bi|i € I}, {Aj|j € J, Jis finite} be modules over a 


ring R.. Then there is isomorphisms of abelian groups: 
1. HomR(A, |],e; Bi) = [<7 Homer(A, Bi). 
2. Homr(®jesAyj, B= @jesHompR(Aj, B). 


Theorem 2.17. Let A, B, {B;|i € I} be modules over a ring R. Then if I is finite 
there is isomorphisms of abelian groups: Homr(A, Sier Bi) = Pier Homp(A, Bi). 
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Proposition 2.18. A direct product of R-modules |]; Ji is injective if only if 
J; is injective for every i € I. 

Corollary 2.19. Let R be an integral domain and let K the field of fractions of 
R. Then K is an injective R-module. 

Corollary 2.20. Let {M)} yea be a family of R-modules. If A is finite and My 
is injective for every \ © A, then ®ycea My is also injective. 


Theorem 2.21. Let M be an R-module. Then M is injective if and only if for 


every short exact sequence 0 —> A 4B ale C —> 0 of R-modules, 





0 —> Homr(C,M) ~> Homr(B, M) + Homr(A,M) — 0 
is also a short exact sequence, where V(f) = fy and O(f) = fd. 


Proposition 2.22. Let R be a ring. A direct sum of R-modules >) ,-;P; is 
projective if only if each P; is projective. 
Proposition 2.23. Every free left R-module is projective. 


Theorem 2.24. A left R-module P is projective if and only if P is a direct 
summand of a free left R-module. 


Corollary 2.25. 


1. Every direct summand of a projective module is itself projective. 


2. Every direct sum of projective modules is projective. 


Lemma 2.26. Let R be a ring with identity. A unitary R-module J is injective 
if and only if for every left ideal L of R, any R-module homomorphism L —> J 


may be extended to an R-module homomorphism R — J: 
Example 2.27. 


1.Q is an injective Z-module by Lemma (2.26) since for every 
Z-homomorphism f :nZ — Q, where nZ is an ideal of Z forO Ane Z, 
there exists a Z-homomorphism g : Z —> Q defined by g(z) = 2f(n) 80 


g(nz) = {n2) f(r) ") = zf(n) = f(nz) for every nz € Z. 
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2. Note that Z is not an injective Z-module since using the Z-homomorphism 
f :2Z —> Z given by f(2z) = z, there is no Z-homomorphism g : Z —> Z 
such that g(2z) = f(2z) for every 2z € 2Z. Otherwise, 1 = f(2) = g(2) 
2g(1), implying that g(1) = 5. However, since g(1) € Z, this is impossible. 








3 n-injective Module 


Definition 3.1. Let C be an category of R-modules, X* € obj(C) for all0 <i <n, 
and d' for al0 <i<n—lisa morphism in C. An R-module M is n-injective if 


the sequence of R-module in C is left n-exact 


d” 





0 Est wy OS en ntl 


if there is M €C the induced sequence of abelian groups 
n+1 dn n qn-1 
0 — Homce(X"™’, M) — Home (X", M) — 
d! 1 d® 0 
... —> Home(X~,M) — Home (X”, M) 


is right n-exact. 


Proposition 3.2. Let C be an category of R-modules, X* € obj(C) for all0 <i< 
n, and d' for all0 <i<n-—1 is a morphism in C. A direct product of R-modules 
Tlicr Ji is n-injective if only if J; is n-injective for every i € I. 

Proof. Let C be an category of R-modules, X* € obj(C) for all 0 <i <n, and d' 
for allO <i<mn-—1 isa morphism in C. The sequence of R-module in C 


d” 


pee ggg aS en en 





is left n-exact. 
Suppose that |],<; Ji is n-injective. To show that, J; is n-injective for each 
i € I. Now if there is |], Jj the induced sequence of abelian groups this sequence 
is 
in m—1 m—2 
0 —+ Home(X"*1, T] Ji) + Home(X", |] 5) 3 Home(X"}, TT 5) “5 


iel iel iel 
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2 , jo 
. © Home(X?, TI] Ji) + Home(Xt, [] Ji) + Home(X°, J] Ji) 
wel iel iel 
is right n-exact. By Theorem |2.16} (1), 


Home(X", |] Ji) = [| Homa(X’, Ji) 
Keall el 


for each 7 € I. Then this sequence 


0 — [J Home(x"!, J) © [] Home(x”, J) “3 T] Home(x"!, J) © 5 
iel iel eal 


©; T] dome (X?, Ji) © T] Home (X?, J) “+ T] Home(X°, Ji) 
iel i€l i€l 
is right n-exact. Then J; is n-injective for each i € I. 
Conversely, suppose that J; is n-injective. To show that, [],<; Ji is n-injective 
for each i € I. Now if there is J; the induced sequence of abelian groups this 


sequence is 


0 — TJ Home(x"!, i) © [I Home(x”, Ji) “3 T] Home (x"!, Ji) £5 
Kaall iel ie 
© J] Home( x, J) Ss T] Home(x?, Ji) + T] Home (X°, Ji) 
iel wel wel 
is right n-exact. By Theorem [2.16} (1). Then this sequence 
0 —+ Home(X"*1, T] Ji) &s Home(X”, [] Ji) “3 Home(X"1, T] J) ©5 


ie€l all all 
. &s Home(X?, T] Ji) “+ Home(X!, T] J) “+ Home(X°, T] 1) 
iel all iel 
is right n-exact. Then [| 











J; is also n-injective. 





Kall 
Corollary 3.3. LetC be an category of R-modules, X* € obj(C) for all0 <i <n, 


and d' for all0 <i <n—1 is a morphism in C. Let R be an integral domain and 
let K the field of fractions of R. Then K is an n-injective R-module. 
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Proof. By Corollary |2.19} & is injective R-module. Let C be an category of 
R-modules, X* € obj(C) for all 0 <i <n, and d' for al0 <i<n-—-lisa 


morphism in C. The sequence of R-module in C 


d 


4 Ses, 22g 





is left n-exact. By Theorem 


0 —> Home(X""!, K) nui Home(X", M) are 


“sy Home(X}, K) &s Home(X°, K) 


is right n-exact. Then K is n-injective module. 














Corollary 3.4. Let {My} yea be a family of R-modules. If A is finite and My is 


n-injective for every X © A, then ®yen My is also n-injective. 


Proof. Let C be an category of R-modules, X* € obj(C) for all 0 <i<n, and d' 


for allO <i<mn-—1 isa morphism in C. The sequence of R-module in C 


d 


xo Fy xt Gye Po yn yt 





is left n-exact. 
Suppose that M) is n-injective. To show that, @)e,M) is n-injective for each 
» € A. Now if there is M) the induced sequence of abelian groups this sequence 


is 
0 —+ @reaHome(X"*1, My) “+ @ye,Home(X", My) li 


n—2 72 
®yeaHome(X""1, My) pas uae cia ®yeaHome(X?, My) 


rl 70 
, @ye, Home(X!, My) 3 @yeq Home(X°, My) 


is right n-exact. If A is finite by Theorem [2.17] 
Home(X", ®yea My) = ©yea Homp~(X*, My) 
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for every \ € A, Then this sequence 
0— Home(X"*1, ®yea My) sl Home (X", ®yea My) 


qr-1 A qr-2 a2 2 
— Home (X »,Prea My) > fen Home (X ,PveaMy) 


71 70 
, Home(X!, ®xeaMx) 2+ Home(X°, yea My) 











is right n-exact. Then ®@)¢,M) is also n-injective. 





Proposition 3.5. Every R-module injective is not n-injective. 


Proof. Let {M)}yea be a family of R-modules. If ®ye, My) is injective, then My 
is injective for every A € A but ®)ye, My, is not n-injective for every 4 € A and 











then, M) is n-injective for every 4 € A. 





Definition 3.6. Let C be an category of R-modules, Y* € obj(C) for all 0 <i< 
n+1, and f? for all 0 < i < nis a morphism in C. An R-module P is n-projective 


if the sequence of R-module in C is rightt n-exact 





yRy Fyfe yn £ yon 
if there is P € C the induced sequence of abelian groups 


0, f° ae ie re 
Home(P,Y°) — Home(P, Y*) — Home(P,Y*) — 


fe n—-1 ae ae id n+1 
.. —> Home(P,Y" ") — Home(P, Y") — Home(P,Y"™) — 0 


is left n-exact. 


Proposition 3.7. Let C be an category of R-modules, Y* € obj(C) for all0 <i< 
n+1, and f* for all0 <i<1n is a morphism in C. A direct sum of R-modules 


Diet P; is n-projective if only if P; is n-projective for everyi € 1 and I is finite. 


Proof. Let C be an category of R-modules, Y* € obj(C) for allO <i <n+1, and 


f? for all 0 <i <n-—1 is a morphism in C. The sequence of R-module in C 


yyy My Aye 
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is right n-exact. 

Suppose that @jc,P; is n-projective. To show that, P; is n-projective for each 
21 € I. Now if there is @j¢7P; the induced sequence of abelian groups this sequence 
is 

0 fe 1 fi 8 f2 
Home (®ietFi, Y") —> Home (Gia Fi, Y") —> Home (Gia Pi, Y“) — 
n-1 mn, 
<< Home (ier Fi, Y”) = Home(iaF;, Y"*') — 0 
is left n-exact. If I is finite by Theorem |2.16| (2) 
Home (GietPi, ¥) = GieHome(P,, Y") 

for every 7 € I, Then this sequence 


0 fe 1 fi a9 f2 
OictHome(P;, Y") — PictHome(P;, Y°) —> Home(P;, Y*) — 


n-1 mn 
Es @ienHome(P,, Y") 2 GicrHome(P;, Y"t1) — 0 


is left n-exact. Then P; is n-projective for each i € I. 
Conversely, suppose that P,; is n-projective. To show that, @jc,P; is 
n-projective for every i € I and I is finite. Now if there is P; the induced sequence 


of abelian groups this sequence is 


SicerHome(P;, Y°) Eu ®ierHome(P;, Y') zu Home(P;, Y”) ca 
ws BierHome(P;, Y”) aa BicrHome(P;, Y"t') — 0 
is left n-exact. If I is finite by Theorem |2.16| (2) 
Home (ie P;, Y°) © OierHome(P;, Y4) 
for every 7 € I. Then this sequence 


0, f° wn fi o f? 
Home (®iet Pi, Y°) —>+ Home (ier P;, Y°) — Home (Piet Pi, Y“) 


n-1 mn 
Home(@iaP,, Y") £> Home(@ieaP;, Y"t) — 0 











is left n-exact. Then @je7P; is also n-projective. 
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4 One Open Problem 


Using the following definitions, can we prove the following theorems about 


n-projective module and free R-module. 

Proposition 4.1. Every free left R-module is n-projective. 

Proposition 4.2. 
1. Every finite direct summand of a n-projective module is itself n-projective. 
2. Every finite direct sum of n-projective modules is n-projective. 

Definition 4.3. Let n be a positive integer. An n-abelian category is an additive 

category C which satisfies the following axioms; 

(AO) The category C is idempotent complete. 


Al) Every morphism in C has n-ker and n-coker. 


( 
(A2) for every monomorphism f° : X° —+ X? in C and, for every n-coker 
(f°, ft,...,f" +) of f°, the following sequence n-exact: 

Bas 


xoR yl A yn FY ynti, 








A2°?) for every epimorphism g” : X" —>+ X"*+! in C and, for every n-ker 
y ep p g 


0 1 po) 


(Qo iad of g”, the following sequence n-exact: 


xoP, xt 2, yn yn, 








Now one can investigate a divisible modules in n-additive abelian category. 
Next one can obtain all of the result of them as we obtained in this paper, and it 


is an open problem. 
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Abstract 


In this paper, we consider a new class of hemivariational inequalities, 
which is called the trifunction bihemivariational inequality. We suggest and 
analyze some iterative methods for solving the trifunction bihemivariational 
inequality using the auxiliary principle technique. The convergence analysis 
of these iterative methods is also considered under some mild conditions. 
Several special cases are also considered. Results proved in this paper can 


be viewed as a refinement and improvement of the known results. 


1 Introduction 


Variational inequalities theory introduced in 1964 by Stampacchia can be 
viewed as a novel and significant generalization of the variational principles. 
The origin of the variational principles can be traced back to Euler, Newton, 
Lagrange and Bernoulli’s brothers. These variational principles have emerged 
as a powerful tool to investigate and study a wide class of unrelated problems 
arising in industrial, regional, physical, pure and applied sciences in a unified and 
general framework. Variational inequalities have been extended and generalized in 
several direction using novel and new techniques. Panagiotopoulos introduced 
the hemivariational inequalities by using the concept of the generalized directional 


derivatives of nonconvex and nondifferentiable functions. This class has important 
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applications in structural analysis and nonconvex optimization. It has been 
shown [7] that, if a nonsmooth and nonconvex superpotential of a structure is 
quasidifferentiable, then these problems can be studied in the general framework 
of hemivariational inequalities. The solution of the hemivariational inequalities 
gives the position of the state equilibrium of the structure. We would like to 
point out that the hemivariational inequalities include the problem of finding the 


difference of two monotone operators, which is itself an interesting problem, see 


[S! [28]. 


Noor and Oettli [16] introduced triequilibrium problems and have shown 
variational inequalities, fixed-point problems, Nash equilibrium problems and 
saddle-point problems can be studied in the framework of triequilibrium problems. 
Thus it is clear that hemivariational inequalities and equilibrium problems 
are different generalizations of variational inequalities. Noor and Noor [I7| 
investigated the trifunction hemivariational inequalities, which can be viewed a 
significant extension of variational inequalities and hemivariational inequalities. 
We would like to emphasize that hemivariational inequality theory provides us 
with a simple, natural, unified, novel and general framework to study an extensive 
range of unilateral, obstacle, free, moving and equilibrium problems arising 
in fluid flow through porous media, elasticity, circuit analysis, transportation, 


oceanography, operations research, finance, economics, and optimization. 


Convexity theory is a branch of mathematical sciences with a wide range 
of applications in industry, physics, social, regional and engineering sciences. 
The general theory of the convexity started soon after the introduction of 
differential and integral calculus by Newton and Leibnitz, although some 
individual optimization problems had been investigated before that. It is worth 
mentioning that variational inequalities represent the optimality conditions for 
the differentiable convex functions on the convex sets. The convex sets and 
convex functions have been extended and generalized in several directions using 
innovative ideas to consider completed problems. See an excellent book by 
Cristescu and Lupsa [3]. Inspired by the research work going on in this field, 
Noor and Noor introduced and and considered a new class of 
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nonconvex sets and nonconvex functions with respect to an arbitrary bifunction. 
This class of nonconvex set is called the biconvex set and the noncovex function 
is called biconvex function. functions is called the biconvex functions. Noor et 
al have studied some basic properties of the biconvex 
functions. It have been shown that the biconvex functions have characterizations 
as the convex functions enjoy. In particular, it have been shown that the 
optimization conditions of the differentiable biconvex functions are characterized 
by a class of variational inequalities, called the bivariational inequalities, see 


and references therein. 


Variational inequalities and hemivariational inequalities have witnessed 
an explosive growth in theoretical advances, algorithmic developments and 
applications across almost all disciplines of engineering, pure and applied 
sciences. There are several methods for solving variational inequalities and 
bivariational inequalities. Due to the nature of the hemivariational inequalities, 
projection and resolvent methods can not be applied for solving hemivariational 
inequalities. In recent years, the auxiliary principle technique is being used to 
suggest and analyze some iterative methods for solving variational inequalities 
and equilibrium problems. Glowinski, Lions and Tremolieres [5] used this 
technique to study the existence problem for mixed variational inequalities, 
whereas Noor and Zhu et al.[32| have used this approach 
to suggest and analyze some iterative methods for solving various classes of 
variational inequalities and equilibrium problems. In this paper, we again use 
the auxiliary principle technique to suggest several new iterative schemes for 
trifunction bihemivariational inequalities. We also prove that the convergence 
of these methods require either pseudomonotonicity or partially relaxed strongly 
monotonicity. These are weaker conditions than monotonicity. As a special 
case, we obtain new iterative schemes for solving bihemivariational inequalities, 
variational inequalities and optimization problem. The comparison of these 


methods with other methods is a subject of future research. 
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2 Preliminaries and Basic Results 
Let H be a real Hilbert space, whose inner product and norm are denoted by (., .) 
and ||.|| respectively. Let kK be a nonempty set in H. 


We now recall some concepts of biconvex sets and biconvex functions, which 
are mainly due to Noor et al. [24]. 


Definition 2.1. The set Kg in H is said to be biconvex set with respect to an 
arbitrary bifunction 6(- — -), if 


u+ rAB(v — u) € Kg, Vu, uv € Kg, € (0, 1]. 


The biconvex set K’g is also called $-connected set. If 6(v — u) = v — u, then 
the biconvex set Kg is a convex set, but the converse is not true. For example, 
the set Kg = R— (-4, 5) is an biconvex set with respect to 3, where 

v-—u, for v>0,u>0 or v<0,u<0 
B(v—u) = 
u-—v, for v<0,u>0 or v<0,u<0. 


It is clear that Kg is not a convex set. 


Remark 2.1. We would like to emphasize that, ifu+ 8(u—u) =v, Vu,v € Kg, 
then B(v—u) = v—u. Consequently, the 3-biconvex set reduces to the convex set 
k. Thus, Kg C K. This implies that every convex set is a biconvex set, but the 


converse is not true. 


Definition 2.2. The function F on the biconvex set Kg is said to be strongly 


biconves, if 
F(ut+rAB(u—u)) < (1-A)F(u) +AF(v) 
—vX(1—2)||B(v — u) ||", Vu, v € Kg, € (0, 1]. 
Note that every convex function is a biconvex, but the converse is not true. 
If \= 5; then the function F’ satisfies 


2u+ B(v—u) 





F( )< j{Flu) + F@)}-v Fle —wIP, Yuu e Ke, 
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which is called Jensen biconvex function. 


If vy = 0, then Definition (2.2) reduces to 
Definition 2.3. The function F on the biconvex set Kg is said to be biconvez, if 
F(u+AB(vu—u)) < 1-A)F(u) +AF(v) Vu,u € Kg,A € (0, 1]. 
We now consider the biconvex function on the interval Iz = [a,a+ 6(b—a)]. 


Definition 2.4. Let Ig = [a,a+ 6(b—a)]. Then F is a biconvex function, if and 
only if, 


1 1 1 
a x a+B(b—a)|>0; a<ax<at+(b—-a). 
F(a) F(#) F(b) 


One can easily show that the following are equivalent: 


1. F is a biconvex function. 











B(b—a) 
F(«x)—F(a) F(b)—F(a) 
3 r—a = B(b—a) 
4 Fla) F(a) FO) <9, 


(Bb—a)))(a—a) + a)—BR—a)(a—a) T BP—a)(e) = 
where = a+ AB(b—a) € [a,a+ B(b—a)]. 





To derive the main results, we need the following assumption regarding the 
bifunction 6(-—-). 


Condition M. The bifunction {(,—, ) is said to satisfy the following assumptions: 


(i). ByBWv—u)) = yewv—u), Vuve Kg, ye R”. 
(iz). Biu—u—yB(v — u)) (l1-yB(v—u), Vu,v € Kg. 


Remark 2.2. Let 6(-—-): Kg x Kg > H satisfy the assumption 


Biv —u) = B(v—z)+Bl(z-—u), Vu,v,2€ Kg. 
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One can easily show that 6(v— u) = 0 u=v, Vu,v € Kg. Consequently 
B(v—u) =0, for v =ue Kg. Also B(u—u)+ B(u—v) =0, Vu,v,z € Kg. This 


implies that the bifunction 6(.—.) is skew symmetric. 


Let f : H — R bea locally Lipschitz continuous function. Let Q be an open 
bounded subset of R”. First of all, we recall the following concepts and results 


from nonsmooth analysis [2]. 


Definition 2.5. Let f be locally Lipschitz continuous at a given point « € H 
and v be any other vector in H. The Clarke’s generalized bidirectional derivative 
of f at x in the direction B(v — u), denoted by f°(x, B(v — u)), is defined as 

f(a +th+tB(u—u))— f(at+h) 


0 ‘ 
_ = | : 
f ey uy) ae a t 





If 6(v — u) = v, then Definition (2.5) reduces to the following concepts which 
are mainly due to Clarke [2]. 


Definition 2.6. Let f be locally Lipschitz continuous at a given point « € H 
and v be any other vector in H. The Clarke’s generalized bidirectional derivative 
of f at x in the direction v, denoted by f°(x,v), is defined as 

f(c+h+tv) — f(xt+h) 


0 : 
z,v) = lim su . 
f ( : ) toot sae t 





The generalized gradient of f at x, denoted Of(x), is defined to be 
subdifferential of the function f°(x;v) at 0. That is 


Of (x) = {we H: (w,v) < fo(a;v), Vu € H.}. 


If f is convex on K and locally Lipschitz continuous at x € K, then Of (x) coincides 
with the subdifferential f’(a) of f at x in the sense of convex analysis , and 
f°(a;v) coincides with the directional derivative f’(x;v) for each v € H, that is, 
f?(a;v) = (f"(@),v), Vue H. 


For a given nonlinear trifunction F(.,.,.) : kgxKgx Kg —> H and a nonlinear 
continuous operator T’: Kg —+ H, consider the problem of finding u € Kg such 
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that 
F(u,Tu, B(v — u)) +f f°(u; B(v —u))\dQ>0, Vue Ke, (2.1) 
Q 


which is called the trifunction bihemivariational inequality. 


Here f°(u; B(v—u)) = f(x, uy; B(v—u)) = f(a, u(x); B(v(@) —u(x))) denotes 
the generalized bidirectional derivative of the function f(z,.) at u(x) in the 


direction u(x) — u(x). 

We now discuss some special cases of the trifunction bihemivariational 
inequalities (2.1). 
(1). If F(u, Tu, B(v—u)) = W(u, B(v—u)), where B(.,.) is a continuous bifunction, 
then problem (2.1) is equivalent to finding wu € Kg such that 


W(u, B(v—u)) + i. f°(u; B(v—u))\dQ>0, Wu € Ke, (2:2) 


which is called the bifunction bihemivariational inequality and appears to be a 


hew one. 


(II). If F(u, Tu, B(v—u)) = (Au, 8(v—wu)), where A is a nonlinear operator, then 
problem (2.1) is equivalent to finding u € Kg such that 


(Au, B(u — u)) +f f°(u; B(v —u))dQ >0, Vu € Ka, (2.3) 
Q 
which is known as the bihemivariational inequality. 
(III). If F(u, Tu, 8(v — u)) = (Au, v — u), where A is a nonlinear operator, then 
problem (2.1) is equivalent to finding u € K such that 


(Au, v — u) +f f(ujv—u)dQ>0, WeEK, (2.4) 
Q 


which is known as the hemivariational inequality introduced and studied by 
Panagiotopoulos {28} [29] in order to formulate variational principles connected 
to energy functions which are neither convex nor smooth. It is has been shown 
that the technique of hemivariational inequalities is very efficient to describe the 


behaviour of complex structure arising in engineering and industrial sciences. 
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(IV). If f is a differentiable convex function, then problem (2.1) is equivalent to 
finding u € Kg such that 


F(u, Tu, B(v — u)) + (f'(u), B(v —u)) > 0, Vu € Kg, (2.5) 


which is known as the mildly nonlinear trifunction bihemivariational inequality 


and appear to be a new one. 


(V). If f =0, then problem (2.1) is equivalent to finding u € Kg such that 
F(u,Tu,B(v—u)) 2>0, We Ke, (2.6) 
which is called the trifunction bivariational inequality. 


In brief, for suitable and appropriate choice of the trifunction, one can obtain 
several classes of bihemivariational and bivariational inequalities. This clear shows 
that the problem (2.1) is more general and flexible and includes the previous ones 


as special cases. 


Definition 2.7. The trifunction F’(.,.,.) and the operator T is said to be: 


(a) jointly bimonotone, if 
F(u, Tu, B(v — u)) + F(v, Tv, B(u—v)) <0, Vu,u € Kg. 
(b) jointly pseudo-bimonotone with respect to fy f°(u; B(v — u))dQ, if 
F(u, Tu, B(v —u)) + [ f°(u; B(v — u))dQ > 0 


=> 
—F(v, Tv, B(u— v)) -[ f°(u; B(v—u))dQ>0, Vu, € Kg. 
Q 


(c) partially relared strongly jointly bimonotone, if there exists a constant 7 > 0 
such that 


F(u, Tu, B(v — u)) + F(v, Tv, B(z — v)) < yI|B(u—z)|I?, Vu,v, 2 € Ke. 


Note that for z = u partially relaxed strongly jointly bimonotonicity reduces 
to jointly bimonotonicity. This shows that partially relaxed strongly jointly 


bimonotonicity implies jointly bimonotonicity, but the converse is not true. 
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Definition 2.8. The function ie f°(u; B(v — u))dQ is said to be partially relaxed 


strongly bimonotone, if there exists a constant a@ > 0 such that 
| f°(u; B(v — u))dQ +f f°(z; B(u— v))dQ < al|B(z—v)|?, Vu,v,2 © H. 
Q Q 


Note that for z = v, partially relaxed strongly bimonotonicity reduces to 


relaxed strongly bimonotonicity. 


3 Main Results 


In this section, we suggest and analyze some iterative methods for solving 
trifunction bihemivariational inequality (2.1) using the auxiliary principle 
technique of Glowinski, Lions and 'Tremolieres involving Bregman distance 
function as developed by Noor [11] {12} [13] [14] [15}, Noor et al. and 
Zhu et al. [32]. 


For the readers convenience, we recall some basic properties of the Bregman 
convex functions [2]. For strongly convex functions f, we define the Bregman 


distance function as 
Biv,u) = f(v) — f(u) — (f'(u),v—u) > alju—ull?, Yuu eK. (3.1) 


It is important to emphasize that various types of function f give different 
Bregman distance function. We give the following important examples of some 
practical important types of function f and their corresponding Bregman distance 


functions. 
Examples 
1. If f(v) = |lv||?, then B(v,u) = ||v — ull, which is the squared Euclidean 
distance (SE). 


2. If f(v) = S37, a; logv;, which is known as Shannon entropy, then its 


corresponding Bregman distance is given as 


Biv, u) = 3 (« log(~*) + Ui — w), 


F 
i=1 : 
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This distance is called Kullback-Leibler distance (KL) and has become a 
very important tool in several areas of applied mathematics such as machine 


learning. 


3. If f(v) = — OL, log u;, which is called Burg entropy, then its corresponding 


Bregman distance is given as 
” UU; 
Biv, u) = log = = 1), 
(v,u) Y (wes +o ) 


This is called Itakura-Saito distance (IS), which is very important in the 


information theory, data analysis and machine learning. 


Remark 3.1. It is a challenging problem to explore the applications of Bregman 
distance function for other types of nonconvex functions such as biconvex, 


k-convex functions, preinvex functions and harmonic functions. 


For a given u € Kg satisfying (2.1), consider the auxiliary problem of finding 
w € Kg such that 


pF(w,Tw,B(v—w)) + (E'(w)— E’(u), B(v —w)) 
=F | f°(w; Blu —w))dQ>0, Vue Kg, (3.2) 
Q 
where p > 0 is a constant and E’(wu) is the differential biconvex function E(w) at 
We Keg. 


We note that, if w = u, then clearly w is solution of the problem (2.1). This 


observation enables us to suggest and analyze the following iterative method for 
solving (2.1). 


Algorithm 3.1. For a given uo € H, compute the approximate solution Unis by 


the iterative scheme 


PF (tns1, Tunsi, B(u — Un4i)) + (E'(uns1) — EB’ (un), 0 — Ung) 
+0 f uns (0 — tings) \d2> 0, Yo € Ko. (3.3) 
Q 
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Algorithm is called the proximal method for solving problem (2.1). In 
passing, we remark that the proximal point method was suggested by Martinet 
[6] in the context of convex programming problems as regularization technique. 


For the recent developments and applications of the proximal point algorithms, 


see and the references therein. 


If F(u, Tu, B(v — u)) = W(u, B(v — u)), then Algorithm collapses to the 
following method for solving the bifunction bihemivariational inequality (2.2). 


Algorithm 3.2. For a given up € H, compute the approximate solution Un+i1 by 


the iterative scheme 


pW (un41, B(v — Un4i)) + (E"(un41) — E"(un), B(v — Un41)) 
+0 | f° (un41; B(v — Un41))dQ > 0, Vu € Ka, 


If F(u, Tu, B(v — u)) = (Au, B(v — u)), then Algorithm [3.1] reduces to: 


Algorithm 3.3. For a given ug € H, calculate the approximate solution un41 by 


the iterative schemes 


(pAtnt + E' (unit) — E'(un), B(v — uns) 
+0 | f° (uns; B(v —tn4i))d2>0, Vu € Kg, 
Q 


is called the proximal point method for solving bihemivariational inequalities (2.3) 


and appears to be a new one. 
If f(z, u) = 0, then Algorithm [3.1] collapses to: 


Algorithm 3.4. For a given up € H, compute the approximate solution Un+1 by 


the iterative scheme 
pF (un41, Tun+1, B(v — Un41)) + (E" (uni) — E’(un), B(v — un41)) > 0, Vu € Kg. 


In brief, for suitable and appropriate choice of the operators and the 
spaces, one can obtain a number of known and new algorithms for solving 


variational-hemivariational inequalities and related problems. 
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We now study the convergence analysis of Algorithm which is the main 


motivation of our next result. 


Theorem 3.1. Let F(.,.,.) and the operator T be jointly pseudomonotone with 
respect to fi. f°(u; B(v —u))dQ. Let E be differentiable strongly biconvex function 
with module uw > 0. Then the approximate solution un+41 obtained from Algorithm 
3.1 converges to a solution u € Kg satisfying (2.1). 


Proof. Let u € Kg be a solution of (2-1). Then 
F(u,Tu, B(v —u)) + iy f(u; B(v—u))\dQ>0, We Kg, 
implies that 
—F(v, Tv, B(u— v)) — [ f(a, u; B(v—u))dQ>0, Wo € Kz, (3.4) 


since F(.,.,.) is jointly pseudomonotone with respect to fg f°(u; B(v — u))dQ. 
Taking v = wu in (3.3) and v = un41 in (8.4), we have 


pF (und; Tun+1, Blu x Un+1)) a (E"(un41) = E"(un), Blu = Un4+1)) 
—0 fF Fun B(u—tn4s))dO (85) 


IV 


and 
—F(Un4i, Tungi, B(u — Un41)) — | f°(u; B(unta — u))dQ > 0. (3.6) 
Q 
We now consider the function Bregman distance function 


Bu, w) 


E(u) — E(w) — (E"(w), B(u — w)) 
u||G(u—w)|/?, (using strongly biconvexity of E). (3.7) 


IV 


where ys > 0 is a constant. 
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Now combining (3.7) and (3.4), we have 


Blu, Un) — B(u,tngi) = E(un+1) — E(un) — (E"(unt1), B(un+1 — Un)) 

+E" (Un+1) — E"(un), B(u — Un+1)) 

}1||B (tent — tn) ||? + (E (unt) — E’(un), B(u — tnt) 
> p||B(tns1 — Un) Il? — oF uns1, Tents, B(u — tn41)) 


) 
-o | 9 (un; B(u — Un+1))dQ 


|B (tin+1 — Un) Il?, 


IV 





IV 


where we have used (3.6). 


If Un41 = Un, then clearly uy, is a solution of the trifunction bihemivariational 
inequality (2.1). Otherwise, it follows that B(u, un) — B(u, Un41) is nonnegative 


and we must have 


lim ||B(un41 — Un)|| =O => Jim unqi = u. 
N+ Oo N—-0o 


Now using the technique of Zhu and Marcotte [20], it can be shown that the 


entire sequence {u,,} converges to the cluster point u satisfying the trifunction 














bihemivariational inequality (2.1). 


It is well-known that to implement the proximal point methods, one has to 
find the approximate solution implicitly, which is itself a difficult problem. To 
overcome this drawback, we now consider another method for solving (2.1) using 


the auxiliary principle technique. 


For a given u € Kg satisfying (2.1), find w € Kg such that 


pF(u,Tu,B(v—w)) + (E(w) — E"(u), B(v— w)) 
+ pf Plu alv —w))dQ, Vue Kg, (3.8) 





where E£’(u) is the differential of a strongly biconvex function E(u) at u € Kg. 
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Note that problems (3.2) and (3.8) are quite different problems.It is clear that 
for w = u, w is a solution of (2.1). This fact allows us to suggest and analyze 
another iterative method for solving trifunction bihemivariational inequality (2.1). 


Algorithm 3.5. For a given uo € H, compute the approximate solution Un4i by 
the iterative scheme 


pF (wn, Twn, Bu —Unti)) + (E' (uni) — E"(wn), Blu — Unt) 
> <p f (uni Blu tnsa))dO, Wu € Kp, (39) 
BF (un; Tun, B(v — Wr)) + (E(wp) — E"(un), B(v — wa) 
> 1 f (uni Bo — wn) AO, Ww € Kg. (3.10) 


Note that for F'(u, Tu, B(v — u)) = W(u, B(v — u)), Algorithm [3.5] reduces to: 


Algorithm 3.6. For a given ug € H, compute the approximate solution uns by 


the iterative scheme 


pW (wn; Biv = Un+1)) 


— 


(E'(un41) — E"(wn), B(v — un+1)) 
=p [ (oon: te —Un+1))dQ, Vu € Kg, 
LW (un, B(v—wn)) + (E'(wn) — E"(un), B(v — wn)) 

> 1 f (um: Btu =a wre iG 


IV 


which is called the predictor-corrector method for solving the bifunction 
bihemivariational inequality (2.3). 


For F(u, Tu, B(v—u)) = (Au, B(v—u)) Algorithm 3.5 collapses to the method 
for solving the bivariational inequalities (2.2). 


Algorithm 3.7. For a given ug € H, compute the approximate solution uns1 by 


the iterative scheme 
ee me ee Cee) -» | bi Bia’ Ks, 
Q 


(Aun + E' (wy) — E'(un), B(v — wn)) > 1 f (uns Blo —Wn))dQ, Vu € Kg, 
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which is called the predictor-corrector method for solving the bihemivariational 
inequalities (2.2). 

If f(.;.) = 0, then Algorithm 3.5 reduces to the following iterative method for 
solving trifunction bivaraiational inequalities (2.5). 


Algorithm 3.8. For a given up € H, compute the approximate solution Un+1 by 

the iterative scheme 

pF (wn, Twn, B(v _ Un+1)) os (E’(un41) _ E" (wn), B(v _ Un+1)) 20, Vue Kg, 
LF (un, Tun, B(v — Wn)) + CE’ (wn) — E' (un), B(v — wn)) > 0, Wo € Kg. 


Similarly for suitable and appropriate choice of the operators and the spaces, 
one can obtain various known and new algorithms for solving hemivariational and 


variational inequalities. 


We now consider the convergence analysis of Algorithm [3.5] using essentially 
the technique of Theorem)3.1} For the sake of completeness and to convey an idea 


of the technique, we sketch the main points. 


Theorem 3.2. Let F(.,.,.) and the operator T be partially relaxed strongly jointly 
bimonotone with a constant y > 0 and let te f°(u;v — u)dQ be partially relaxed 
strongly bimonotone with a constant a > 0. If EF is strongly biconver function 
with modulus 8 > 0 and0 < p < B/(at+y), O0< p< B/(a+7), then the 


approximate solution Un41 obtained from Algorithm [3.5] converges to a solution 


uée Keg of (2.1). 


Proof. Let u € Kg be solution of (2.1). Then 
p{P(u,Tu, (vw) + f f(us Bw w))dO} > 0, Wwe Ky (3.11) 
Q 


{FP (u, Tu, B(v — u)) + i f°(u; B(v — u))dQ} > 0, Vu € Kg, (3.12) 
Q 
where p > 0 and pz > 0 are constants. 


Setting v = Un41 in (3.11) and v = u in (3.9), we have 


AF (u, Tu, Bunt _ u)) + iE foe B(Un+1 _ u))dQ} = 0. (3.13) 
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and 


pF (wn, Twn, B(U—Unt1)) + (B"(unti) — E"(wn), Bu — Un+1)) 
Zz 2th f° (x, wn; B(u— Uny1))dQ. (3.14) 


As in Theorem 3.1] and from and (3.13), we have 


B(u, wn) — B(u, un+1) 
= E(unti) — E(wn) — (E"(un41), B(Un41 — Wn) 
+(B"(Un41) — B'(wn), B(U — Un+41)) 
> pl|B(un41 — Wn) ||? + CE’ (unt1) — E" (wn), B(u — unt1)) 
> pl|B(unga — Wa)||? — PF (wn, Twn, B(u — Un41)) 


2 of f (tn: B(t6 — tiny1))€0 
pl|8(tn41 — wn) 
—p{F (wn, Twn, B(u — tigi )) tr Fu, Tu, 2B (unsa _ u))} 
=p 7 7(20; Blting1 — d+ [ Pei Bt teagan 
1l|B(0in41 — Wn) ||? — pla +7) |B (einer — wn)? 

= {u—platy)}IlBleing1 — wn), 








IV 


IV 


where we have used the fact that F(.,.,.) and de f°(.;.)dQ are partially relaxed 


strongly bimonotone with constants @ > 0 and y > 0 respectively. 


In a similar way, we can obtain 


B(u, Un) — B(u, wn) = {8 — pla +7) } [|B (wn — un) |l?. 


If Unt1 = Wn = Un, then clearly uy is a oe of the ae hemivariational 
inequality (2.1). Otherwise, forO << p<-z oa and0Q<p< oe, it follows that the 
sequences B(u, Wy) — B(u, Un4i) and B(u, un) — B(u, wy) are nonnegative and we 


must have 


Jim, ||B(un+1 — Wn) || =0 and Jim, ||B(wn — un) || = 9. 
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Thus 
Tim, [3 (unga — tn) || = lim, [18mg — wn) || + ian 180m — ten) || = 0 


Now using the technique of Zhu and Marcotte [20], it can be shown that the 
entire sequence {u,,} converges to the cluster point u satisfying the trifunction 
bihemivariational inequality (2.1). 














We now suggest and analyze some new iterative methods for solving 
the trifunction bihemivariational inequality (2.1) using the auxiliary principle 
technique of Glowinski, Lions and Tremolieres [10] without the Bregman distance 


function as developed by Noor [16-24]. 
For a given u € Kg satisfying (2-1), find w € Kg such that 
pF(u,Tu, B(v—w)) + (w—-—u,v—w) 
+ of f°(u; B(u—w))dQ>0, Woe Kg, (3.15) 
where p > 0 is a constant. Problem is known as the auxiliary trifunction 
bihemivariational inequality. We note that if w = u, then clearly w is a solution 


of the (2.1). This observation enables us to suggest and analyze the following 
iterative method for solving (2.1). 


Algorithm 3.9. For a given uo € H, compute the approximate solution Un+1 by 


the iterative scheme 


pF (wn, Twn, B(U—Wn)) + (Unti — Wn, Vv — Un+1) 
+ | f°(u; B(v = Un41))dQ > 0, Vu € Kg (3.16) 
Q 
nF (un, Tun, v— Un) + (Wn — Un, UV — Wn) 


+ nf #(u lv ~ wn) do >0, Wu € Kg, (3.17) 


where p > O and 7 > O are constants. Algorithm is called the 


predictor-corrector method for solving the trifunction bihemivariational inequality 
2). 


We now study the convergence analysis of Algorithm [3.9] 
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Theorem 3.3. Let u € Kg be a solution of and Un+1 be the approximate 
solution obtained from Algorithm . If F(.,.) ts partially relaxed strongly 
monotone with a constant a > 0 and the operator {, f°(u;—, )dQ is partially 


relaxed strongly monotone with a constant y > 0, then 





lJunt1 — ll? < |leun — al? — (1— 2p(a +7))llunsi — wall? (3.18) 
l|eun — GI]? < [lun — GI|? — (1 — 28(a + 7))|lwn — ual”. (3.19) 


Proof. Let u € Kg be a solution of (2-1). Then 
oF(u,Th,B(v—a) + p|f f(a; B(v - Wd >0, Ve Kg (3.20) 
Q 
nFG@TA,S(v— 0) + nf Pla Glv-W)AN>0, Wwe Ky, (8.21) 
Q 


where p > 0 and 7 > 0 are constants. 


Now taking v = un+1 in (3.20) and v = w in (3.16), we have 


pF (t, Tu, Unti — t) + | f°(u; Bunda — u))dQ > 0 (3.22) 
Q 
pF (wn, TWn, u— Wn) + (Un+1 — Wn, — Un+1) 


+ [ f(us (di — ting1))A > O. (3.23) 


Adding and (3.23), we have 


(Unt1 — Wn, U — Un41) 


> -p[F (Wn) T Un, Bll — tn) + F(t, TH, B(tin41 — @)} 
a i. f(t; B(tin41 — @))d0+ [ f(s B(G — tny1)) 4} 
> -(a+7)pllunsi — wall’, (3.24) 


where we have used the fact that F’(.,.,.) is relaxed strongly monotone with 


constants a > 0. 


Recall the following result, 


2(u,v) = |lu + v||? — |Jull? — lull’, Va,b © A, (3.25) 
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Setting u = G — Un4i and V = Un41 — Wn in (3.25), (3.24) can be written as 





lJen¢2 — Gl? < [la — wall? — 1 — 2(a + 7)p)|un4a — wall’, 


the required (3.18). 
Taking v = @ in (3.21) and v = wy in (3.17), we obtain 





nF (8, TH, B(t, —@)) +9 [ f(u; Blwn —2))dX>0 (3.26) 
nF (Un, Tun, B(% — un)) + (wn — Un, ti — Wn) 
+f #%(uni 8(a ~ ty) de> 0. (3.27) 


Adding (3.26), (3.27) and rearranging the terms, we have 
(Wn — Un, U — Wn) = —B(a+7)|lUn — Wnll’, (3.28) 


since F(.,.,.) and fo, f°(u;—))dQ are partially relaxed strongly monotone with 
constants a > 0 and y > 0 respectively. 


Now taking v = wp — Un and u = % — wy in (8.25), (3.28) can be written as 


|e — all? < a — un||? — (1 — 2(a +7) 8) [wn — Un”, 


the required (3.19). 














Theorem 3.4. Let H be a finite dimensional space and let 0 < p< 1/2(a+ 7), 
0< 8 <1/2(a+ 7). [fu © Kg is a solution of (1) and un4i is an approximate 
solution obtained from Algorithm|3.10, then 


nate tin) = 


Proof. Let ti € Kg be a solution of (2.1). Since 0 < p < 1/2(a+7) and 0 < 
B <1/2(a+7), it follows from (3.18)and (3.19) that the sequences {||w, — al} 


and {||t — un||} are nonincreasing and consequently {un} and {wn} are bounded. 
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Furthermore, we have 


co 


So = 2(a + 7)p)llunta — wall? < [lwo — all? 
n=0 

oo 

S51 - 2(a + 7) A) len — Unll? < luo — a)?, 
n=0 


which implies that 


lim ||unz1 — Wn|| =O and lim |/wp — Un|| = 0. 
T—> OO n oO 
Thus 
lim ||tnti— Un|| < lim |lung1 — wall + lim ||wn — unl| = 0. (3.29) 
n—0o n— oo A SS 


Let d@ be a cluster point of {u,} and the subsequence {un,} of the sequence {Up } 
converge to ti € H. Replacing wy, by un, in (3.15), (3.16) and taking the limit 
nj; —> oo and using (3.29), we have 


F(a, Tt, v — a) +f f(a; Blu —%))dQ>0, WeEK, 
Q 
which implies that @ solves the trifunction bihemivariational inequality (2.1) and 
l]ent1 — GI? < [en — all? 


Thus, it follows from the above inequality that the sequence {uy} has exactly one 
cluster point @ and 


lim (u,) =a, 
Mm—> oo 











the required result. 





In recent years, inertial proximal methods have been suggested and 
analyzed for maximal monotone operators associated with the discretizations 
of the differential equations in times, whereas Noor has used the auxiliary 
principle technique to suggest an inertial method for variational inequalities, 


the converges of which requires only pseudomonotonicity, which is a weaker 
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condition than monotonicity. This clearly improves the convergence criteria of 
the inertial proximal method. We again use the auxiliary principle to suggest and 


analyze an inertial proximal method for solving the trifunction bihemivariational 
inequality (2.1). 
For a given u € Kg satisfying (2.1), consider the problem of finding w € Kg 
such that 
pF(w,Tw, Bu — w)) + (w-u-a(u—u),v—w) 
+0 | f°(u; B(v—w))dQ>0, We Kg, (3.30) 
Q 





where p > 0 and a > 0 are constants. 


It is clear that, if w = u, then wu is a solution of (2.1). This fact 
allows us to suggest and analyze an iterative method for solving the trifunction 
bihemivariational inequality (2.1). 


Algorithm 3.10. For a given uo € H, compute the approximate solution Un+1 


by the iterative scheme 


PF (Uns, Tuns1, B(v — Un+41)) 


+(Un+1 — Un — On (tn = tinct): v— scl ) 


+0 | f°(u; B(v — Ungi1))dQ > 0, Vu € Kg. (3.31) 
Q 





For a, = 0, Algorithm reduces to : 


Algorithm 3.11. For a given uo € H, compute the approximate solution Un+1 


by the iterative scheme 


pF (Unt, Pung, ¥ — Unda) + (Ung — Uns VY — Und) 


+0 f fu (0 —unss))d0> 0, Yo e Ko, 
Q 


which is known as the proximal method for solving trifunction bihemivariational 


inequality (2.1). 
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In a similar way for F(u,Tu,v — u) = (Au,v — u), one can obtain a 
number of new and known proximal methods from Algorithm for solving 
bihemivariational inequalities and its special cases. This shows that the 
new methods suggested in this paper are unifying one and more general than the 


previous ones. 


For the convergence analysis of Algorithm |3.11| we need the following result. 
Theorem 3.5. Let u € Kg be a solution of and let Un4i be the 
approximate solution obtained from Algorithm If the trifunction F(.,.,.) 
is pseudomonotone with respect to |, f°(.;—))dQ and the operator i Pas) de 


is monotone, then 


[Una — all < |lun- all? — |[Unta — Un — On (Un — Un—1)|l" 


40g — aif? lfn1 — a]? + 2Ifty — ty l2}. (8.32) 
Proof. Let ti € Kg be a solution of (2.1). Then 
—F(v,Tv, B(u— v)) +f f°(G; B(v-—@dQ>0, We Kz, (3.33) 
Q 


since F(.,.,.) is pseudomonotone with respect to fo, f°(.;.))dQ. 


Taking v = un+1 in (3.33), we have 


F(Un+i,; Tun+1, B(U — Un41)) + A f°(G; B(@ — Unq1))dQ > 0. (3.34) 


Now taking v = 4 in (3.31), we obtain 


pF (Un41; TUn+1; u— Un+1) ae Cth — Un — On (Un =, me u— pei) 


n p fF (ums; B(@— tnss))Ad0 > 0. (3.35) 
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From (23), (24) and using the monotonicity of f, f°(.;.))dQ we have 


nad — Un — An(Un a ici) u— Un+1) 


IV 


— pF (Un41, Tungi, U — Ung) — pI? (un41;& — Un41) 


- _, PteccetG tea 
Q 


IV 


+ [1G Bln — @)aO} > 0, (3.36) 
Q 
which implies that 
(Un+1 — Un, U — Un+1) > An(Un — Un—1, tt — Un + Un — Un41). (3:30) 


Using (3.25) and rearranging the terms in (3.37), one can easily obtain (3.32), the 
required result. 














Theorem 3.6. Let H be a finite dimensional space. Let un+1 be the approximate 
solution obtained from Algorithm |3.9 and tu € Kg be a solution of (2.1). If there 
exists a € (0,1) such thatO<an<a, Wne Nand So, an|lUn — Un—ill? < 


oo, then limy 309 Un = U. 
Proof. Let t. € Kg be a solution of (2.1). First we consider the case a, = 0. Using 
the technique of Theorem [8.3] we can prove that limp,—,oo Un = @. 

Now we consider the case a,, > 0. From (3.32), we have 


Co 
S> ||2en44 — Un — An (Un = Un—1) 7 
n+1 





[oe] 
S |[uo — Gl)? +S {alleen — G||? + 2llun — un-1I?} < 00, 


n= 


which implies that 
im || nt4 — Un — An(Un = Un—1)||? = 0. 


Repeating the arguments as in Theorem one can easily show that 











limn+oo Un = U, the required result. 
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Conclusion 


In this paper, we have introduced and studied the trifunction bihemivariational 
inequalities. Several special cases are discussed as applications of the trifunction 
bihemivariational inequalities. The auxiliary principle technique is used to 
suggest several implicit and explicit iterative methods for solving the trifunction 
bihemivariational inequalities, Convergence criteria of the proposed methods is 
discussed under suitable mild conditions. Results obtained in this paper continue 
to hold for the special cases. Comparison of the proposed methods with other 
methods need further efforts. The ideas and techniques of this paper stimulate 


further research in these dynamic fields 


Acknowledgement. We wish to express my deepest gratitude to our teachers, 
students, colleagues, collaborators and friends, who have direct or indirect 


contributions in the process of this paper. 
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Abstract 


Recent researches have revealed the capability of Machine Learning (ML) techniques to 
effectively detect fraud in electronic banking transactions since they have the potential to 
detect new and unknown intrusions. A major challenge in the application of ML to fraud 
detection is the presence of highly imbalanced data sets. In many available datasets, 
majority of transactions are genuine with an extremely small percentage of fraudulent 
ones. Designing an accurate and efficient fraud detection system that is low on false 
positives but detects fraudulent activity effectively is a significant challenge for 
researchers. In this paper, a framework based on Hidden Markov Models (HMM), 
modified Density Based Spatial Clustering of Applications with Noise (DBSCAN) and 
Synthetic Minority Oversampling Technique Techniques (SMOTE) is proposed to 
effectively detect fraud in a highly imbalanced electronic banking dataset. The various 
transaction types, transaction amounts and the frequency of transactions are taken into 
consideration by the proposed model to enable effective detection. With different number 
of hidden states for the proposed HMMs, simulations are performed for four (4) different 
approaches and their performances compared using precision, recall rate and Fl-Score as 
the evaluation metrics. The study revealed that, our proposed approach is able to detect 
fraudulent transactions more effectively with reasonably low number of false positives. 
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1. Introduction 


E-banking is a form of banking where funds are transferred as exchange of electronic 
signals rather than cash, checks, or other types of paper documents [1]. Over the last few 
decades, E-Banking has redefined the way banking is conducted across the globe and the 
use of electronic payments platforms has continued to experience significant growth. It 
allows customers a 24-hour access to their accounts with the ability to transfer funds, 
perform on-line payments and apply for loans and other financial products virtually [2]. 


Fraud can be defined as any premeditated act of criminal deceit, trickery or 
falsification by a person or group of persons with the intention of altering facts, in order 
to obtain undue personal monetary advantage [3]. Unfortunately, fraud cases relating to 
cyber-crime perpetrated through E-banking resulted in an actual loss of GH¢14.31 
million and therefore presents a unique challenge to individuals and financial institutions 
that offer those services [4]. To address this problem, financial institutions employ 
various fraud prevention tools such as real-time transaction authorization, transaction 
verification codes, transaction alerts, rule-based detection among others. Fraudsters 
however are adaptive, and given time, they devise several ways to circumvent such 
protection mechanisms [5]. There is therefore the need to implement enhanced 
technologies and systems that can detect fraud in real-time effectively in order to 
maintain the viability of these electronic payment systems where fraudsters constitute a 
very inventive and fast-moving fraternity. As preventive technology changes, so does the 
technology of criminals and the way they go about with their fraudulent activities [6]. 
While it is necessary to detect and possibly prevent fraudulent transactions, it is also very 


critical to ensure genuine transactions are executed successfully. 


One of the most important techniques for intrusion/anomaly detection based on 
machine learning is using Hidden Markov Models (HMM) which are machine learning 
algorithms consisting of hidden states and observable outputs for modelling probability 
distributions over sequences of observations. The hidden state layer is a stable Markov 
chain and its state probability and state transition probability are decided from the initial 
state probability vector z and the state transition probabilities. Observable output layer is 
decided from the observed symbols probability matrix which is derived from the 
observed symbols of each hidden state [7]. 


The application of HMMs ranges from speech and image recognition, intrusion/fraud 
detection to motion/action analysis in videos among others and is generally characterized 
by the following [8]; 
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1. The number of hidden states in the model denoted as N. The state at a specific time 
t is denoted by q;. 


2. The number of unique observation symbols denoted as M. 
3. A transition probability between states denoted by a matrix A = la; jl, where: 
aij = P(qta1 = 5) | Ge = Si). (1) 
Also, Vhiay=1 1S5iSN. (2) 
4. An emission probability matrix, B = [bj (k)], where 
By(k) = PVk = ae 15; = 4) (3) 
wed (kK) = 1, 1Sj/ <N. (4) 
5. An initial probability for each state denoted by the vector z = [7], where 


t= Pq =S),. Dian =1. (5) 


In recent decades, many research communities have been working toward HMM- 
based intrusion detection mainly because of its ability to detect new and unknown 
intrusions and usage in real-time applications by processing data streams on-the-fly. 
HMMs also allow for the usage of heterogeneous data sources as input, and visual 
representation of acquired knowledge relative to the other techniques of machine 
learning. 


Over the past few years, the use of Electronic banking platforms has continued to 
experience significant growth and has redefined the way banking or E-commerce is 
conducted across the world [9]. On the other hand, fraudulent Electronic banking and E- 
commerce activities are becoming more and more sophisticated and challenging leading 
to massive financial losses. Effective and efficient detection of Electronic banking fraud 
is therefore regarded as one of the major challenges to all financial institutions, and is an 
increasing cause for concern [2]. 


According to the Bank of Ghana 2019 banking industry fraud report, fraud cases 
relating to cyber-crime perpetrated through electronic banking and mobile banking 
platforms accounts for the highest value of attempted fraud amounting to GH¢ 50.54 
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million with actual loss of GH¢14.31 million [4]. From available literature, majority of 
the works in the area of HMM-based fraud detection in Electronic banking focuses only 
on payments to merchants for goods and services . Transaction amounts are mostly taken 
as observation symbols and the types of items purchased considered as the hidden states 
of the proposed Hidden Markov Models. In related studies conducted by [10], [11], [12], 
[13], [14], and [15], techniques such as Neural Network , Bayesian Network , Dempster- 
Shafer theory, Support Vector Machine etc. are employed which incorporated other 
forms of electronic banking options such as remote funds transfers and deposits. 
However, all these proposed techniques perform classification based on a single 
transaction while relying on domain-expert features without considering a sequence of 
transactions to make a decision hence producing high levels of false positives. 


A large number of false positives may translate into bad customer experience and may 
lead customers to take their business elsewhere. A major challenge in applying ML to 
fraud detection is presence of highly imbalanced data sets. In many available datasets, 
majority of transactions are genuine with an extremely small percentage of fraudulent 
ones. Designing an accurate and efficient fraud detection system that is low on false 
positives but detects fraudulent activity effectively is a significant challenge for 
researchers. 


This proposed research seeks to develop and implement an improved fraud/intrusion 
detection system for both debit and credit transactions in electronic Banking using 
Hidden Markov Models by incorporating the various electronic banking platforms 
employed by customers, transaction amounts and the frequency at which these 
transactions occur. To determine the transaction profile of customers, the Density-based 
Spatial Clustering of Applications with Noise (DBSCAN) which is capable of 
discovering clusters of different shapes and sizes from a large amount of data containing 
noise and outliers was employed. Synthetic Minority Oversampling Technique (SMOTE) 
was also employed to handle the imbalanced class problem typical of Electronic banking 
datasets. 


The rest of the paper is organized as follows: In Section 2, we present a review of 
related works. The methodology adopted for the study is outlined in Section 3. Detailed 
experimental results and discussion to establish the efficiency of the proposed approach is 
presented in Section 4. Finally, we conclude the paper with some discussions in 
Section 5. 
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2. Literature Review 


Fraud Detection in Electronic Banking is understudied in literature perhaps due to 
security and data privacy concerns. We will begin by considering related works in 
electronic banking in general and then consider those specifically related to the use of 
credit cards which has been given considerable attention by researchers. 


[10] presents a fraud detection system for online banking where differential analysis 
is used to obtain local evidence of fraud where a significant deviation from normal 
behavior indicates a potential fraud. The Dempster’s rule of combination is applied to 
these evidences for final suspicion score of fraud. Their main contribution is a fraud 
detection method based on effective identification of devices used to access accounts and 
assessing the likelihood of being a fraud by tracking the number of different accounts 
accessed by each device. However, their system performs poorly for higher number of 
Hidden states and also when users’ transaction patterns changes frequently. 


[16] considered transaction amounts and purchases types as the emission symbols and 
hidden states respectively of the proposed HMM for online banking FDS. The model is 
trained with the normal behavior of an account holder using Baum-Welch algorithm and 
a One-time-Password is sent to the Customers contact number for authorization if an 
incoming transaction violates the behavior sequence. Although, the accuracy of their 
system was close to 72 percent over a wide variation in the input data, False Positives 
was still high especially when the transaction data is highly skewed. A fraudulent 
transaction could still go through if a fraudster has access to a customer’s phone. 


[11] incorporates several advanced data mining techniques for online banking fraud 
detection by building a contrast vector for each transaction based on its customer’s 
historical behavior sequence. A novel algorithm, Contrast Miner, was introduced to 
efficiently mine contrast patterns and distinguish fraudulent from genuine behavior, 
followed by an effective pattern selection and risk scoring that combines predictions from 
different models. Results from experiments on large-scale real online banking data 
demonstrated that the proposed system achieves substantially higher accuracy and with 
lower false positives by incorporating domain knowledge and traditional fraud detection 
methods. 


[12] rather modeled the sequence of operations in online banking transaction 
processing using HMMs and described how it could be used for the detection of frauds. 


The observation sequence length is fixed to two (2) whilst changing sequence length 
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for training i.e., changing dataset length from 10 to 80 with difference of ten. Simulation 
results revealed that, although the complexity of the system also increases for increased 
observation sequence length, the accuracy of the proposed system is close to 60% with 
reduced false Positive rate. 


The work done by [14] employed HMMs and k-means algorithm for detecting fraud 
in online banking transactions. In their proposed model, a variable is used to keep the 
number of transactions within a period of time before and after each transaction as well 
as the quantified amounts as the observation symbols. If an incoming transaction is not 
accepted by the trained HMM with sufficiently high probability, it is considered 
fraudulent. The feasibility of their proposed model is demonstrated through simulation 
experiments using real-world bank transaction data. In the case of enough historical 
transactions, their model performs well for low, medium frequency and amount of user 
groups. An efficient Prior determination of the number of clusters is considered a major 
challenge in their proposed approach. 


Specifically on fraud detection relating to the use of Credit Cards, [17] considered 
purchase types and transaction amounts as hidden states and observation symbols 
respectively in their proposed HMM. In order to estimate the model parameters, the K- 
means clustering algorithm is employed to determine the spending profile of cardholders. 
An incoming transaction is considered fraud if it is not accepted by the HMM with a 
significantly high probability. Experimental results revealed that, their proposed model 
recorded an accuracy close to eighty (80) percent over a wide variation of the data. An 
efficient prior determination of the number of clusters and significant number of false 
positives were considered the major challenges in their proposed approach. 


[18] performed a comparative analysis of intrusion detection models on highly 
skewed credit card data based on Decision Trees, Random Forest, Support Vector 
Machines (SVM) and logistic regression. The original sample was randomly partitioned 
into k-equal sized subsamples where a single subsample is retained as the validation data 
for testing the model, and the remaining k — 1 subsamples used as training data. With the 
four basic metrics employed, namely True positive (TPR), True Negative (TNR), False 
Positive (FPR) and False Negative (FNR) rates, Simulation results using dataset provided 
by ULB machine learning revealed that, Logistic regression and Random forest shows 
the most precise and high accuracy in the area of credit card fraud detection but requires 
very large dataset for training and also suffers from the imbalanced dataset problem even 
after preprocessing. 
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In order to reduce the number of false positives, [19] proposed a model based on 
automated feature engineering to automatically derive behavioral features based on the 
historical data of a credit card associated with a transaction. A total of 237 features for 
each transaction was generated, and a random forest was then employed to learn a 
classifier. One important feature of their proposed model is that, it also utilizes the 
distance between two locations transactions on an account has occurred and whether they 
occurred in person or remotely is established. The proposed model was tested on data 
from a large multinational bank and compared to existing solutions and revealed that, on 
an unseen data of 1.852 million transactions, false positives was reduced by about 54%. 
However, since their models Perform classification based on a single transaction there 
was a performance degradation when transaction pattern of users changes frequently. 


3. Methodology 


There is generally a very limited number of public datasets on electronic banking for 
research purposes mainly due to personal and security concerns. In this research, a 
Kaggle provided dataset of simulated mobile based transactions is adopted. As detailed in 
Table 1, the dataset is highly imbalanced due to the fact that only 8,312 transactions out 
of the almost 6 million transactions are labeled as fraud. 


Table 1: Details of the Paysim Dataset adopted for the study. 


























Transaction # of Genuine # of Fraudulent Total 
Type Transactions Transactions 

TRANSFER 528812 4097 532909 
CASH-OUT 2233384 4116 2237500 
CASH-IN 1399284 0 1399284 
DEBIT 41432 0 41432 
PAYMENTS 2151494 0 2151494 
TOTAL 6354407 8213 6362620 














‘CASH IN’ and ‘CASH OUT” represents an increase in account balance of a 
customer as a result of cash inflow and a decrease in account balance as a result of cash 
outflow respectively. ‘TRANSFER’ refers to movement of money between users whilst 
‘PAYMENT? represents the settlements made for goods and services to merchants. 
‘DEBIT’ as used in this context signifies the sending of money from a mobile service 
(electronic wallet) to a bank account. 
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3.1. Data pre-processing 


To effectively evaluate the performance of our proposed models on the highly class 
imbalanced dataset, Synthetic Minority Oversampling Technique (SMOTE) is employed 
to generate virtual training records by linear interpolation for the fraudulent transactions 
by randomly selecting one or more of the k-nearest neighbors for each specific fraudulent 
transaction. After the oversampling process, the data is reconstructed and then the 
proposed Hidden Markov Models is applied on the processed data. Specifically, the 
sampling rate is set to 73000 %. 


The proposed SMOTE technique as adopted in this study is presented in Algorithm 1. 
Algorithm 1: The SMOTE algorithm 
Procedure SMOTE (f,R,k) 


Input: Number of Fraudulent Transactions (f); Amount of SMOTE R %; Number of 
nearest neighbors k 


Output: (R/100) * f synthetic Fraudulent Transactions 
1: The number of Fraudulent Transactions is set to f 


2: For each y € f, the k-nearest neighbors of y are obtained by calculating the 
Euclidean distance between y and every other sample in set f. 


3: The sampling rate R is set according to the imbalanced proportion. 


4: For each y € f, R examples are randomly selected from its k-nearest neighbors, and 
they construct the set f;. 


5: For each example y; € f; (t = 1,2,3,...,R), generate a new sample as in equation 
6 below: 


y' =y t+rand(0, 1) * ly — yl. (6) 
3.2. Identifying transaction profile of customers 


For optimal training of the our proposed Hidden Markov Models, a modified Density 
Based Spatial Clustering of Applications with Noise (DBSCAN) and the K-means 
clustering algorithms are executed on each customer’s previous transactions by 
considering the amount and frequency of transactions. K-means is an unsupervised 
learning algorithm for grouping a given set of data based on their similarity where the 
numbers of clusters are fixed a priori. The grouping is performed by minimizing the sum 
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of squares of distances between each data point and the centroid of the cluster to which it 
belongs to [20]. The DBSCAN clustering technique however filters out outliers and 
discovers clusters of arbitrary shapes [21]. We modified the DBSCAN algorithm by 
adding a step that computes the centroid of each cluster later to be used to dynamically 
convert an incoming transaction into an observation symbol in the fraud detection 
process. 


The proposed DBSCAN technique as adopted in this study is presented as in 
Algorithm 2. 


Algorithm 2: The DBSCAN algorithm. 
DBSCAN (dataset, d, minpts) 


Input: A set of points, dataset, distance threshold d, and the minimum number of 
points required to form a cluster, minpts. 


Output: A set of clusters representing the various observation symbols 
l:n = 1, #initialise the cluster index to 1 

2: For each unvisited point pt in dataset, mark pt as visited 

3: Find the neighboring points, N of p 

: If |N| >=minpts then N = NUN’ 


: if p’ is not a member of any cluster, Mark it as noise. 


Nn Nn £ 


: Compute the centroid f of each cluster using (7) below 
ee De 7 
= X;, 
= Da 4 (7) 
xjECi 
where n; is the number of points in cluster c;. 


Spending profiles of accountholders are determined at the end of the clustering step. 
Let w; be the percentage of total number of transactions of an accountholder, then, the 
spending profile p of an account holder, J is determined as in (8): 


p(8) = arg Max(‘p;). (8) 


The cluster number to which most of the transactions of the account holder belongs 
to represents the spending profile of the account holder. The computed centroids are used 
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to generate the observation symbol for a new transaction @ (denoted by @,,,) is defined as 
in (9). 


dm = Varg; min|m — nj|. (9) 


The i" transaction on account A, denoted as Po is suspected to be an outlier if it 
does not belong to any cluster in the set C’ where y refers to the frequency of transaction. 
If the average distance of the amount p of an outlier transaction Pe from the set of 


existing clusters in C’ is Wg, then its level of deviation 0, is given as in (10): 


Wa-é€ 
—— if |N < MinPt. 

os = ft if INe(p)| < MinPts es 
0 otherwise. 


The key idea of the modified DBSCAN algorithm is that for each point p in a cluster 
C;, there are at least a minimum number of points (MinPts) in the e-neighborhood of that 
point p denoted as N,(p) ie., the density in the e-neighborhood has to exceed some 
threshold. The proposed K-means algorithm as adopted in this study is presented in 
Algorithm 3. 


Algorithm 3: K-Means Clustering Algorithm 

1: Specify the number of clusters to assign 

2: Randomly initialize k centroids 

3: Repeat 

4: Expectation: Assign each point to its closest centroid 

5: Maximization: Compute the new Centroid (Mean) of each cluster 
6: Until the Centroid Positions do not change 


Specifically, for this research, the set C = {low-frequency low-amount, low- 
frequency medium-amount, low-frequency high-amount, medium frequency low-amount, 
medium-frequency medium-amount, medium-frequency high-amount, high-frequency 
low- amount, high-frequency medium-amount, high-frequency high-amount} denotes the 
clusters. 


The set R = {transaction_amount, frequency_of_transaction} represents the set of 
attributes used to generate these clusters. 


To compute the probability of an observed sequence, 0 = (9,01, 02,.--,O7—1) With 
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respect to our Hidden Markov Model A, where v = (vo,Uj,U2,...,Ur—1) represents the 
various hidden states, the definition of the emission transition Matrix is defined as in 
(11); 
P(a|v,A) = byo (G9) by1(%) «- byr-1(Or-1)- (11) 
The Initial transition vector, x and State Transition Matrix, A are also defined as in 
(12) and (13); 
P(v|A) = TyAy yt ++ AyT-2,vT-1 (12) 
and 


P(o,v|A) = P(oluv,A)PQA). (13) 


By summing over all possible state sequences, (14) through (16) is obtained 


P(olv) = YP, v|A)x (14) 
= YPColv, A)PWIA)y (15) 
= YT voPyo(Fo)Avv1 oo Ayr-2vT-1Dyr-1(Or = VG: (16) 


The probability of the observation sequences denoted as e;(i), where the system is in 
state q; at time f is defined in (17). 


er (i) = Py, 02, 03, ... Of, Ue = GilA). (17) 

e, (i) is then calculated recursively as in (18) to (20): 
1. Let eg (i) = 1b; (00) (18) 

2. Fori = 0,1,2,...,N —1andt = 1,2,...,T — 1, we compute 
ec(i) = LjZ0 Ler-1 ayi] bi (oe) (19) 
3. Equation (20) is obtained from (19) 
P(olA) = Yo er-1 0. (20) 
3.3. Training the proposed HMMs 


The transaction amounts are categorized into a Low (l) = (0; 100], Medium 
(m) = (100; 500], and High (h) = (500; Transaction Limit] values. The frequency at 
which these transactions occur on a particular is also categorized into a Low (Less than 5 
times a month), Intermediate (Between 5 and 10 times a month), and High (at least 10 
times a month) are also considered by our proposed model. For example, if an 
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accountholder performs about seven (7) transactions with the month with an average 
value of say 300, then the corresponding observation symbol is medium-frequency 
medium-amount (mm). 


The various transaction types are considered the internal states whilst the transaction 
amounts combined with the frequency at which they occur denoted as 
{lL lm, lh, ml,mm, mh, hl, hm, hh} represents the observation symbols of our proposed 
Hidden Markov Model. 


After formulating the hidden states and observation symbol, a hybrid optimization 
algorithm as presented in Algorithm 4 comprising the Baum-Welch, Particle Swam and 
Genetic Algorithms is used to effectively train the proposed models. 


Algorithm 4: A hybrid algorithm for optimizing the parameters of the proposed 





HMMs 

1. Initialize the parameters (A, B, 7r) using the spending profile of the customer 

2.a,() = joo [ee-1 G) ajil bi (Oz) (21) 

3. BO = Yio [Bet ais] bj (O¢41) (22) 
~\ _ OBO 

4. ye) = Es (23) 
Oe (i)aijbj(Or4+1) Bt4+1G) 

5.7.) =—~“Foa (24) 

6.For0 < ij < N-1 

7.1; = Yo) (25) 
=VST-2y () )/ST-2y. (j 

8. aij = Le=0 Ve G@I)/Xt=0 Ve O (26) 

9. bk) = Yeeto,1,..7-1 Vt Q/yiety- @O (27) 

Op=k 
10. Go to 6 


11. After 100 iterations, each solution becomes a chromosome for the genetic 
procedure and the fitness function P(O|A) is applied. 

12. A multiple point crossover and mutation is performed to select the best 50 
solutions for the next generation which are then positioned as particles in a search space 


using the PSO technique. 
13. The position and velocity of each particle during each iteration is updated using; 
14. X,(t +1) = X;(t) + V;(t) (28) 
15.V,(t +1) = wV,(t) +7 UC((0,1]J) (CX +i @ — XH) 
+r, U([0,1]) ( %(0)- X:(0). (29) 
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16. Compare the best solution each particle and the best position of the entire group 
and make appropriate adjustments 

17. Termination Criteria Reached? If yes go to 18, otherwise go to 13 

18. Output A, B and x 


3.4. Fraud detection 


To effectively classify an incoming transaction as fraudulent or otherwise, sequence 
of observation symbols, say 0 = 01,02,...,0, are extracted from the training data of an 
account holder and its probability of acceptance, 0, is computed by the model as in (30) 
by employing (20) 


01 = p(04, 02, 03,...,0,|A). (30) 


An incoming transaction occurring at time is converted to an observation symbol 
denoted as o;, 1 using (9) is used to replace the first observation symbol, o, and its 
probability of acceptance by the model denoted as 02 is also computed as in (31). 


02 = P(o2, 03,04, +» Or44|A). (31) 


The newly generated transaction is classified as fraud and if the difference between 0, 
and 0, denoted as Ad is above a predefined threshold (¢) as in (32). 


Ad/d, > €. (32) 


A genuine transaction is added to the sequence permanently to contribute to 
determining the validity or otherwise of the next transaction since transaction behavior of 
an accountholder could be dynamic. Otherwise, the transaction is declined, and the 
symbol is discarded. 


Due to the highly imbalanced nature of the dataset, precision(p), recall(R) and F1- 
scores (F) as presented in equations (33) to (35) respectively are used as evaluation 
metrics (Wedge et al. [19]). Tp, Fp, Fn represent True Positives, False Positives and 
False Negatives respectively. 


p =Tp/(Tp + Fp) (33) 
R=Tp/(Tp+Fn) (34) 
F =2+*(P*R)/(P+R). (35) 


Precision quantifies the number of correct positive predictions made whilst Recall 
refers to the number of correct positive predictions made out of all possible positive 
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predictions. F-Measure however provides a way to combine both precision and recall into 
a single measure. 


4. Simulation Results and Discussion 


For different number of hidden states, four (4) sets of simulations were performed in 
two (2) stages using Python programming and their performance compared. For all the 
four sets of experiments, the proposed hidden Markov models were executed in the 
second stage. In the first stages of the first and second set of experiments, K-means and 
the modified DBSCAN algorithms were executed respectively. In first stage of the third 
set of experiment, both SMOTE and K-means techniques were employed whereas 
SMOTE and the modified DBSCAN clustering techniques were executed during the last 
set of experiments. The dataset was loaded and divided into two, 80% of it is used for 
training and evaluation whilst the rest is held back for validation. 


4.1. Precision comparison 


The precision of the four approaches is presented for different number of hidden states 
and presented in Figure 1. It is very clear that our proposed approach 
(SMOTE+DBSCAN+HMM) performed better for the various hidden states. Applying 
only the modified DBSCAN clustering technique with Hidden Markov Models 
performed relatively better than that of employing K-Means. It is also worth noting that, 
relatively higher values of precision scores were recorded when the SMOTE technique is 


adopted. 
1.20 
1.00 
0.80 
z —_ * 
S —e— PROPOSED 
wn 
2 0.60 
oq —<—HMM+KMEANS 
* 0.40 HMM+DBSCAN 
=—=@— HMM+SMOTE+KMEANS 
0.20 
0.00 


2 3 4 5 
NUMBER OF HIDDEN STATES 


Figure 1: Comparison of the precision of the four (4) different approaches for different 
number of hidden states. 
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4.2. Recall rates comparison 


A comparison of the Recall rates of the four (4) different approaches for different 
numbers of hidden states are presented in Figure 2. It is evident that approaches that 
employed the SMOTE technique appear to perform relatively better. Similarly the 
modified DBSCAN clustering technique performed better as compared to the K-means. 


It can also be observed that, for higher values of N, all the approaches performed well 
except for those that employed Only K-means and Hidden Markov models without 
handling the class imbalance classification. 


1.20 


i an 


wy 0.80 
= —— PROPOSED 
= 0.60 
a: == HMM+KMEANS 
oO 
~ 0.40 HMM+DBSCAN 
—= HMM+SMOTE+KMEANS 
0.20 
0.00 


2 3 4 5 
NUMBER OF HIDDEN STATES 


Figure 2: Recall rates of the four (4) approaches for different number of Hidden states. 
4.3. F-measure comparison 


In Figure 3, the Fl-score for the four approaches are presented various Hidden states. 
It is observed that, higher Fl-scores are obtained when the modified DBSCAN clustering 
technique is used as compared to using the K-means. Also, approaches that incorporated 
the SMOTE technique performed better. Employing both SMOTE and the modified 
DBSCAN clustering algorithms appears to perform relatively better than the other. All 
approaches performed relatively better when the number of hidden states is 3. 
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Figure 3: Fl-scores of the four (4) different approaches for various number of Hidden 
states. 


5. Conclusion 


In this research, an improved electronic banking fraud detection framework based on 
Hidden Markov Models (HMM) and modified Density Based Spatial Clustering of 
Applications with Noise (DBSCAN) is proposed and implemented. The Synthetic 
Minority Oversampling Technique (SMOTE) is also employed due to the highly class 
imbalance nature of the dataset adopted. With different numbers of hidden states, 
simulations were performed two stages for four (4) different approaches in Python and 
their performance compared. For all the four sets of experiments, the proposed hidden 
Markov models were executed in the second stage. In the first stages of the first and 
second set of experiments, K-means and the modified DBSCAN algorithms were 
executed respectively. In first stage of the third set of experiment, both SMOTE and K- 
means techniques were employed whereas SMOTE and the modified DBSCAN 
clustering techniques were executed during the last set of experiments. 


Generally, our proposed approach (SMOTE+DBSCAN+HMM) performed relatively 
better for all the various hidden states in terms of precision, recall and Fl-Scores. 
Employing the modified DBSCAN clustering technique to determine the spending profile 
of customers and subsequently performed relatively better than using the K-Means 
algorithm since it filters out most of the easily recognizable fraudulent transactions before 
the proposed HMMs are applied. It is also evident from the simulation analysis that, the 
SMOTE technique effectively handles the class imbalance classification necessary to 
achieve improved performance. 
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Abstract 


In this paper, we introduce the generalized Oresme sequence and we 
deal with, in detail, three special cases which we call them modified 
Oresme, Oresme-Lucas and Oresme sequences. We present Binet’s formulas, 
generating functions, Simson formulas, and the summation formulas for these 
sequences. Moreover, we give some identities and matrices related with these 


sequences. 


1 Introduction 


The Oresme sequence, {On}n>0, was introduced by Nicole Oresme (1320-1382) 
in the 14-th century. Oresme found the sum of the rational numbers formed by 
the terms 0, 5, ‘, 3 + ss ne .. 37 These numbers form a second order sequence 
and are defined by the recurence relation 





1 1 
On+2 = On41 _ qOm Oo = 0, O71 = 2° 


In [4]. Horadam presented a history and obtained an abundance of properties 
of these numbers. Oresme numbers have many interesting properties and 
applications in many fields of science (see, for example, (12 ]3[4)9)). 

The purpose of this article is to generalize and investigate these interesting 
sequence of numbers (Oresme numbers). First, we recall some properties of 


Fibonacci numbers and its generalizations, namely generalized Fibonacci numbers. 
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The Fibonacci numbers, Lucas numbers and their generalizations have many 
interesting properties and applications to almost every field such as architecture, 
nature, art, physics and engineering. The sequence of Fibonacci numbers {Fy }n>0 
is defined by 

Fy, =Fyit Fro, n>2, Fo =0, = 1, 


and the sequence of Lucas numbers {L,,}n>0 is defined by 
Ln = In-1 + In-2, n= 2, [Io = 2, L) = 1. 


The generalization of Fibonacci sequence leads to several nice and interesting 
sequences. The generalized Fibonacci sequence (or generalized (r, s)-sequence 
or Horadam sequence or 2-step Fibonacci sequence) {W,(Wo, Wi; 7, 5)}n>o (or 
shortly {W,}n>0) is defined (by Horadam (6)) as follows: 


Wr =TWr_-1 + sWn_2, W=a,wWw,=b, n>2 (1.1) 


where Wo, W, are arbitrary complex (or real) numbers and r,s are real numbers, 
see also Horadam and Soykan [12]. 

For some specific values of a,b,r and s, it is worth presenting these special 
Horadam numbers in a table as a specific name. In literature, for example, the 
following names and notations (see Table 1) are used for the special cases of r, s 


and initial values. 


Table 1: A few special case of generalized Fibonacci sequences. 























Name of sequence W,(a, 0; r, 5) Binet Formula OEIS/10] 

(4)" | (54)" 

2 2 

Fibonacci W,,(0,1;1,1) = Fr A000045 

( ) Vi ; 
Lucas Wii hi, (2335) sh (54) 000032 

1++/2)" — (1—2)” 
Pell Wr(0,1;2,1) = Py (1+ v2) -(- v2)" a oooiag 
2/2 

Pell-Lucas Wr(2,2;2,1)=Qn (14+ V2)" +(1- V2)” 4002203 
Jacobsthal W,(0,1;1,2) = Jn ae" A001045 
Jacobsthal-Lucas W,,(2,1;1,2) = jn 2” + (—1)” A014551 
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Here, OEIS stands for On-line Encyclopedia of Integer Sequences. 
The sequence {W,,}n>0 can be extended to negative subscripts by defining 


r 


1 
4g == Way + We 
WN: 3 —n-n + [Wa(n-2) 


for n = 1,2,3,... when s £ 0. Therefore, recurrence (1.1) holds for all integer n. 
Now we define two special cases of the sequence {W,,}. (r,s) sequence 
{Gn(0,1;r,s)}n>0 and Lucas (r,s) sequence {H,(2,r;7,5)}n>0 are defined, 


respectively, by the second-order recurrence relations 


Gn+2 = rGn+1 Ae sGn, Go = 0, G1 = I, (1.2) 
An+2 = rAnsy Tr 5H, Ho = 2, Ay =?7, (1.3) 








The sequences {Gn}n>0, {Hn}n>o and {En}n>so can be extended to negative 
subscripts by defining 


r 1 
Gua Ge n— —G_ n—2)) 
ae) ho 

r 1 
A_, = ——A_(n- = ty) 
gi (n-) + 5 4-(n-2) 


for n = 1,2, 3, ... respectively. Therefore, recurrences (1.2)-(1.3) hold for all integer 
n. 
Some special cases of (r,s) sequence {G,,(0,1;7,5)}n>0 and Lucas (r,s) 


sequence {H,,(2,7;1T, )}n>o0 are as follows: 
1. G,(0,1;1, 1) = Fy, Fibonacci sequence, 
2. H,(2,1;1,1) = Ly, Lucas sequence, 
3. Gn(0,1; 2,1) = P,, Pell sequence, 
4. H,(2,2;2,1) = Qn, Pell-Lucas sequence, 
5. G,(0,1;1,2) = Jn, Jacobsthal sequence, 


6. H,(2,1;1,2) = jn, Jacobsthal-Lucas sequence. 
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The following theorem shows that the generalized Fibonacci sequence W,, at 


negative indices can be expressed by the sequence itself at positive indices. 


Theorem 1. For n € Z, for the generalized Fibonacci sequence (or generalized 
(r, s)-sequence or Horadam sequence or 2-step Fibonacci sequence), we have the 


following: 


(a) 
Win = (1)°*1s-"(W, — HaWo) 
= (=1)""s""(W, — H.W): 


(b) 
(—1)"+15-” 


W_n= 
—W? + sWe + rWoW, 


((2W, — rWo)WoWrsi — (W? + sW@)W,). 

















Proof. For the proof, see Soykan Theorem 3.2 and Theorem 3.3]. 


The following theorem presents sum formulas of generalized (r,s) numbers 


(generalized Fibonacci numbers). 


Theorem 2. Let x be a real (or complex) number. For all integers m and j, for 
generalized (r,s) numbers (generalized Fibonacci numbers), we have the following 


sum formulas: 





(a) If (—s)™x? — cH» +140, then 


((—3)"" 2 — Hig) Winns 
sci + (-8)" 0" Winntj—-m + Wj — (-8)™ tWj—-m 
er . (1.4 
2 Wink+j (—s)™x? —rH, +1 a 











(b) If (—s)"x? — 2H», +1 = u(x —a)(x — b) = 0 for some u,a,b € C with u 4 0 
anda #b, t.e., x =a orx=b, then 
(x(n +2) (—s)" — (n+ 1)Hn)2"Wy mn 
~ —s)™ 1)" Wrn+j—-m — (—s)"” j—m 
So 2 Wonks suet (—s)™ (n+ 1)2"Winsj (-s)"”" W; . 
k=0 





2(—s)" 2 — Hm 
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(c) If (—s)™x? — cH» +1 = u(x —c)? = 0 for some u,c € C with u $0, ie, 





x=c, then 
(n +1) ((—s)™ (n + 2)a” — nz” Hm) Wnt; 
a: anti T\ilaaV" 2? Winton 
Soong = (n +1) (=8) - 
= 2 (—s) 














Proof. It is given in Soykan Theorem 4.1]. 


Note that (1.4) can be written in the following form 
((—s)”” so Fy a Wes + 3)" WN tea 


Sl ekw e + (Hm — (—8)" #)W3 ~ (=8)"" tWj-m 
k=1 “~— ( OM tH, + 1 








co 
We give the ordinary generating function }> W,x” of the sequence {W,,}. 


n=0 


co 
Lemma 3. Suppose that fw, (2) = >> W,x” is the ordinary generating function 


n=0 


co 
of the generalized Fibonacci sequence {W,}n>0. Then, >> Wna” is given by 
n=0 





us Wot+(W, —rW 
S> Wr” = at te ome (1.5) 
c=. l—-—raz—-sx 











Proof. For a proof, see {12} Lemma 1.1]. 





Binet’s formula of generalized Fibonacci sequence can be calculated using its 


characteristic equation (the quadratic equation) which is given as 
zg? —r2—s=0. (1.6) 


The roots of characteristic equation are 


ate ps (1.7) 








where 
A=r’+4s 


Earthline J. Math. Sci. Vol. 7 No. 2 (2021), 833-3867 


338 Yiiksel Soykan 





and the followings hold 


a+ 8 


r, 


Q 
iy 
I 


= 8; 


(a—f)? = (a+) —4a8 =r? 4+4s. 


1.1 Binet’s Formula for the Distinct Roots Case 


In this subsection, we assume that the roots a and 6 of characteristic equation 
(1.6) are distinct. Using these roots and the recurrence relation, Binet’s formula 


can be given as follows: 


Theorem 4 (Distinct Roots Case). Binet’s formula of generalized Fibonacci 


numbers 1s 


bra” b 3” = bya” _ bb” 





Wr = | = 1.8 
a=——p° p=a a—p (1.8) 
where 
by = WwW, — BWo, bo = Wi = aWo. 
(1.8) can be written in the following form: 

W,, = Aja” + AgB” (1.9) 

where 
A - WwW, = Wo. A _ WwW, — oWo 

a—p B-a 

Note that 


(W? = sWe = rW,Wo) 


A,A 
I: 2 —(r2 + 4s) ’ 





Aj tA. = Wo. 


We next find Binet’s formula of generalized Fibonacci numbers {W,,} by the 


use of generating function for Wp. 
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Theorem 5 (Binet’s formula of generalized Fibonacci numbers). 
dja” = dz 3” 


WB) Bem 


(1.10) 





where 


dy = Woa + (Wy = rWo), 
dg = Woh + (Wi —rWo)Z. 

















Proof. For a proof, see Theorem 1.2]. 


Note that from (1.8) and (1.10) we have 
Wi-8Wo = Woa+(Wi-—rW), (1.11) 
W,-aW, = Wob+ (Wy = rWo) 8. (1.12) 


For all integers n, (r,s) and Lucas (r,s) numbers (using initial conditions in 
(1.8) or (1.10)) can be expressed using Binet’s formulas as 
qQ” nr 
+ p ; 
(a—B) (B-a) 
Hy, = ga iar 





G, = 


respectively. 


1.2 Binet’s Formula for the Single Root Case 


In this subsection, we assume that the roots a and 6 of characteristic equation 


(1.6) are equal, i.e., a = 8. So (1.6) can be written as 
a? —ra—s=(2—a)* =2* -2ar+07 =0. 


Note that in this case, 


. 
a= -, 
2 
r = 2a, 

2 

eee 

S$ = -a@ 1? 
re+4s = 0. 
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Using the root a@ and the recurrence relation, Binet’s formula can be given as 


follows: 


Theorem 6 (Single Root Case). Binet’s formula of generalized Fibonacci numbers 


as 


Wr = (D, + Don)a” (1.13) 
where 
dD, = Wo, 
1 
Dy = (Wi —aWo). 














Proof. For a proof, see Soykan [13]. 
Note that (1.13) can be written as 


1 
(Wo + — (Wi — aWo) n)a” 





Wn = — 
a 
= (nW,-—a(n—1)Wo)a”™! 
r\n-l 
= (nW - 5 (n-1)W) (5) 
We also see that 
ow, — 
DiDo = Wo(2W, wo). 
r 
D,+ Dy = M1 


For all integers n, (r,s) and Lucas (r,s) numbers (using initial conditions in 


(1.8) or (1.10)) can be expressed using Binet’s formulas as 


G, = na”, 


Ayn = 2a”, 
respectively. 
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2 Generalized Oresme Sequence 


In this paper we consider the case r = 1,5 = —j. A generalized Oresme sequence 


{Wr}n>0 = {Wn(Wo, W1)}n>0 is defined by the second-order recurrence relations 
1 
Wr = Wn — GWn-2 (2.1) 


with the initial values Wo = cp, W1 = c, not all being zero. 


The sequence {W,,}n>0 can be extended to negative subscripts by defining 
W_n = 4W_(n—1) — 4W_(n_2) 


for n = 1,2,3,.... Therefore, recurrence (2.1) holds for all integers n. 
Eq. (1.13) can be used to obtain Binet formula of generalized Oresme numbers. 


Binet formula of generalized Oresme numbers can be given as 


Wr = (D, + Don)a” (2.2) 
where 
dD, = Wo, 
1 
Dy = (Wi awe): 
(a4 
ie., 


1 
Wr = (Wo + mi (Wi _ aWo) nja”. 


Here, a= 6 = 5 are the roots of the quadratic equation 


1 
w—2+7=0. (2.3) 


ie. the roots of characteristic equation (2.3) are equal. Note that 


a+6 = 1, 
1 

aB a q’ 
a-P = 0. 
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and 


1 
WwW, = (Wo + 2 (m -_ 5H) n) x=. 


The first few generalized Oresme numbers with positive subscript and negative 


subscript are given in the following Table 2. 


Table 2: A few generalized Oresme numbers. 








n Wr W_n 

0 Wo Wo 

1 Wi 4Wo — 4W, 

2 W, — 4Wo 12Wp — 16W, 

3 3W, — Wo 32W  — 48W, 

4 5W, — Wo 80Wo — 128W 
5 2Wi-3Wo 192Wo — 320W; 
6 4mM-AWo 448Wo — 768W, 
7 GWi-AWo 1024Wo — 1792, 
8 gWi- Wo —-2304Wo — 4096, 
9 Wi- g@Wo —-5120Wo — 9216, 
10 52gWi— qaqgWo 11264Wo — 20480, 





Now we define three special cases of the sequence {W,,}. Modified Oresme 
sequence {Gn}n>o0, Oresme-Lucas sequence {Hy}n>9 and Oresme sequence 


{On}n>o are defined, respectively, by the second-order recurrence relations 


1 

Gni2 = Gnit—FGn, Go = 0,41 =1, (2.4) 
1 

Hn42 = Hnti~GHn, Ho = 2, Hi = 1, (2.5) 
1 1 

On+2 = On+1 — qOnm Oo = 0, O71 — 2 (2.6) 





The sequences {Gn}n>0, {Hn}n>o0 and {On}n>0 can be extended to negative 
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subscripts by defining 


Gn = 4G_(m-1) — 4G_(n_2), 
Ayn = 4H_(n-1) — 4H_(n-2), 
O_n = 40 _(n-1) = 40_(n—2); 


for n = 1,2,3,... respectively. Therefore, recurrences (2.4)-(2.6) hold for all 
integers n. 
Next, we present the first few values of the modified Oresme, Oresme-Lucas 


and Oresme numbers with positive and negative subscripts: 


Table 3: The first few values of the special second-order numbers with positive 


and negative subscripts. 





























n 0 2 3 4 5 6 7 8 9 10 11 
Gr 0 1 I i 3 18 is i i6 256 256 i034 
G4 4 16 48 128 320 768 1792 4096 9216 20480 45056 
H, 2 1 1 1 1 1 1 1 1 1 1 1 
2 4 8 16 32 64 128 356 512 To24 
A_n 4 8 16 32 64 128 256 512 1024 2048 4096 
GO @ 2 1 3 1 5 3 7 1 9 5 1 
2 2 8 4 32 32 128 32 512 512 2048 
O_n 2 8 24 64 160 384 896 2048 4608 10240 22528 





For all integers n, modified Oresme, Oresme-Lucas and Oresme numbers (using 


initial conditions in (2.2) can be expressed using Binet’s formulas as 








_ n-1 n 
n Ras na a 9Qn—-1? 
1 
= mo 
n a 2a = gn—-1? 
nr 
— nm 
On = no = Qn? 


respectively. 


Earthline J. Math. Sci. Vol. 7 No. 2 (2021), 833-367 


344 Yiiksel Soykan 








Note that 
G,. = 20% 
Ce = sat = "Hn, 
On = ar = pin. 


[oe] 
Next, we give the ordinary generating function )> W,,2” of the sequence 


n=0 


co 
Lemma 7. Suppose that fw, (x) = >> W,x” is the ordinary generating function 


n=0 


co 
of the generalized Oresme sequence {Wn}n>0. Then, >> Wyx” is given by 





n=0 
3) Wat =4x HOt aoe (2.7) 
= (x — 2) 





Proof. In Lemma|3| take r= 1,5 = — 








1 
a 





The previous Lemma gives the following results as particular examples. 


Corollary 8. Generated functions of modified Oresme, Oresme-Lucas and 


Oresme numbers are 


n ~— ’ 
= (« — 2) 
CO 
4 
ae = ag) 
n= 


[oe] 
2 
as 


respectively. 


Proof. In Lemma [7| take W, = G, with Go = 0,G, = 1, Wy, = Ay with Hp = 
2,4, =1 and W, = Op, with Op = 0,O;1 = 5; respectively. 
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3 Simson Formulas 


There is a well-known Simson Identity (formula) for Fibonacci sequence {F;,}, 
namely, 
Raha 7? =(1 


which was derived first by R. Simson in 1753 and it is now called as Cassini 


Identity (formula) as well. This can be written in the form 


Phat fy 


=(-1)". 
Fy Fn-1 








The following theorem gives generalization of this result to the generalized Oresme 


sequence {W,,}n>0.- 


Theorem 9 (Simson Formula of Generalized Oresme Numbers). For all integers 


-(i) 


Proof. For a proof of Eq. (3.1), see Soykan [11], just take s = —+, 
7 


n, we have 
Wi Wo 


Wasi Wn 
an ; (3.1) 
Wo W-1 


Wh Wr 1 























The previous theorem gives the following results as particular examples. 


Corollary 10. For all integers n, modified Oresme, Oresme-Lucas and Oresme 


numbers are given as 








Gantt Gn | — -1 
CG, Gaa..| ~~ 2 

| Hnit Hn | _ 
Hy, Hn-1 j 
Onti On | _ -1 
On Ona | ~ a? 








respectively. 
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4 Some Identities 


In this section, we obtain some identities of generalized Oresme, modified Oresme, 
Oresme-Lucas and Oresme numbers. First, we can give a few basic relations 


between {W,,} and {Gy}. 
Lemma 11. The following equalities are true: 


Wn = 16(2Wo —3W1) Gris —4(5Wo — 8W1) Gnas, (4.1) 
Wy = 4(8Wo —4W1) Gn43 — 4(2Wy —3W1) Gnso, 

Wr = 4(Wo- Wi) Grta — (83Wo — 4W1) Gn4i, 
Wi 
Wi 


WoGn+1 — (Wo — Wi) Gn, 
1 
= WG, - qWoGn-1, 





and 
(Wo —2Wi)?Gn = 64(Wo — 3W1)Wr4ya — 16(3Wo — 8W1)Wris, 
(Wo —2W1)*Gn = 16(Wo — 4W1)Way3 — 16(Wo — 3W1) Wate, 
(Wo —2Wi)*Gn = —-16WiWri2 — 4(Wo — 4W1) Wrst, 
(Wo —2Wi)*Gn = —4WoWni1 +4WiWn, 
(Wo —2W1)?Gn = —4(Wo — W1)Wn + WoWn-1. 


Proof. Note that all the identities hold for all integers n. We prove (4.1). To show 


(4.1), writing 


Wr =ax Gnia tb x Gris 


and solving the system of equations 


Wo = ax G4+6x G3 
Wi 


axGs5+bx G4 


we find that a = 16 (2Wo — 3W,), b = —4(5Wo — 8W). The other equalities can 


be proved similarly. 
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Note that all the identities in the above Lemma can be proved by induction 


as well. 


Next, we present a few basic relations between {H,,} and {W,}. 


Lemma 12. The following equalities are true: 














(Wo —2Wi)Hn = —32Wrz4 + 16Wr4s, 
(Wo —2Wi)Hn = —16Wrz3 + 8Wr4o, 
(Wo —2Wi)Hn = —8Warso+4Wrst, 
(Wo —2Wi)Hn = —4Way1 +2Wa, 

(Wo —2Wi)Hn = —2Wn +Wr-1. 


Now, we give a few basic relations between {W,,} and {O,}. 


Lemma 13. The following equalities are true: 


Wr = 32(2Wo — 3W1) Ona — 8(5Wo — 8Wi) Onis, 
Wr = 8(3Wo — 4Wi) Ons3 — 8 (2Wo — 3W1) On+2, 
Wn = 8(Wo —W1) On+2 — 2(3Wo — 4W1) Ons, 
Wr = 2WoOn+1 a2 (Wo = W) On, 
1 
Wn = 2WiOn — 5WoOn-1, 
and 
(Wo —2W1)? On = 32(Wo — 3W1)Wrsa — 8(3Wo — 8W1) Was, 
(Wo —2W1)?O, = 8(Wo —4Wi)Wr+3 — 8(Wo — 3W1) Wr +2, 
(Wo —2W,)? On = —8WiWp+2 —2(Wo — 4W1)Wrst, 
(Wo —2W,)?O, = —2WoWn41+2WiWn, 
1 
Wo — 2W, On = —2(Wo —-W1)W, + =<WoWn-_1- 
2 


Now, we give a few basic relations between {G;,}, {Hy} and {O,}. 
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Lemma 14. The following equalities are true: 


and 


and 


and 


and 

















He. =: W668 cx: 
Ha. = SCAG Ss: 
Fie = AC 5 = 9G 4, 
Hy = 26.4 =¢.. 
= Gn 5Gn-1, 

Os. Be SG a Eben. 
Of, = “8Cysg HEC HS, 
Ox, i 0G SOs oe 
On = 5G 

Gn = —96On44 + 640n43, 
Gy =" 820944 20,55, 
Gh = “80s 48075, 
Gn = 20n, 


= 320 n+44 = 160 p43, 
= 160 n+3 — 80n+2; 


An, 
An, 
A, = 80n4 
Hy, 
An, 








12 — 40n41, 
= 4On+1 20n,; 
= 207 = One4; 

n 
n = on = nHp, 

nr 
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We now present a few special identities for the generalized Oresme sequence 


Theorem 15 (Catalan’s identity of the generalized Oresme sequence). For all 
integers n and m, the following identity holds: 


m2 


= aoe (WO 2W,)°. 


Wra+mWn—m a. WwW = 


Proof. We use the identity 


1 1 
Wr = (Wo + 2 (m _ 5H) n) Xx an" 














As special cases of the above theorem, we have the following corollary. 


Corollary 16. For all integers n and m, the following identities hold: 


(a) GnimGn—m _ G? = ee 


(b) HnimHn—-m — H2 =0. 


n 


(€) Ont Ones SO? = 28. 


Note that for m = 1 in Catalan’s identity of the generalized Oresme sequence, 


we get the Cassini identity for the generalized Oresme sequnce. 


Theorem 17 (Cassini’s identity of the generalized Oresme sequence). For all 


integers n, the following identity holds: 


1 


ean WO 2W,)°. 


WriiWn-1 — W? = 
As special cases of the above theorem, we have the following corollary. 


Corollary 18. For all integers n, the following identities hold: 
(a) Gin41Gn—-1 — G2 = — siz. 
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(b) HnsiHn_1 — H? =0. 
(c) On+10n-1 = O? = — sir. 
The d’Ocagne’s, Gelin-Cesaro’s and Melham’ identities can also be obtained 
by using 
1 
Wr= (Wo +2 (m _ 5H) n) Xx an" 
The next theorem presents d’Ocagne’s, Gelin-Cesaro’s and Melham’ identities of 


generalized Oresme sequence {W,, }. 


Theorem 19. Let n and m be any integers. Then the following identities are 


true: 


(a) (d’Ocagne’s identity) 











m—n 
WmtiWn —WmWn+ = un = a (Wo — 2Wi)? 
(b) (Gelin-Cesaro’s identity) 
Wr+2WntiWn—-1Wn-2-— Wr 
1 2 2 6 2 2 2 
oS (4(5n*—4)W7+(5n*—10n+1)W5 —4(5n*—5n—4)W1 Wo) (Wo — 2W1)°. 


(c) (Melham’s identity) 


Wrt1Wnt2Wn+6—W2,3 = (2(7n+15)W—(7n+8)Wo) (Wo — 2W1)?. 


1 
~ 98n+9 





Proof. Use the identity W,, = (Wo + 2 (Wi - $Wo) n) x oat 











As special cases of the above theorem, we have the following three corollaries. 
First one presents d’Ocagne’s, Gelin-Cesaro’s and Melham’ identities of modified 


Oresme sequence {G,,}. 


Corollary 20. Let n and m be any integers. Then the following identities are 


true: 
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(a) (d’Ocagne’s identity) 


mn 
Gm+iGn - GmGn+1 = ——— 
(b) (Gelin-Cesaro’s identity) 
5n? — 4 
Gint+2Gn+1Gn—-1Gn_2 = Gh = aces 
(c) (Melham’s identity) 
7n+ 15 
Gin41Gn42Gn46 — G43 = — 


Second one presents d’Ocagne’s, Gelin-Cesaro’s and Melham’ identities of 


Oresme-Lucas sequence {H;,}. 


Corollary 21. Let n and m be any integers. Then the following identities are 


true: 


(a) (d’Ocagne’s identity) 
AmtiAn _ Am Ani = 0. 


(b) (Gelin-Cesaro’s identity) 





An42An41An, 1An, 9 —- Hi=0. 


(c) (Melham’s identity) 


An+1Hn42Hn+6 — 343 = 0. 


n 


Third one presents d’Ocagne’s, Gelin-Cesaro’s and Melham’ identities of 


Oresme sequence {O,}. 


Corollary 22. Let n and m be any integers. Then the following identities are 


true: 
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(a) (d’Ocagne’s identity) 





m—n 
Om4+10n — OmOn+1 = — 
(b) (Gelin-Cesaro’s identity) 
bn? — 4 
On+2On410n-10n—2 Or, = ( 94n . 
(c) (Melham’s identity) 
7n +15 
On+10n420n+6 — OF 43 = aa 


5 On the Recurrence Properties of Generalized 


Oresme Sequence 


Taking r = 1,s = —} in Theorem |1| (a) and (b), we obtain the following 


Proposition. 


Proposition 23. For n € Z, generalized Oresme numbers (the case r = 1,8 = 


—5) have the following identity: 


Wen = (yt (=F) Wy = Hao) 





From the above Proposition, we have the following corollary which gives 
the connection between the special cases of generalized Oresme sequence at the 
positive index and the negative index: for modified Oresme, Oresme-Lucas and 
Oresme numbers: take W, = G, with Go = 0,G,; = 1, take W, = A, with 
Ay = 2,H, = 1 and W, = O, with Op = 0,O, = 5, respectively. Note that in 
this case H, = Hn. 


Corollary 24. For n € Z, we have the following recurrence relations: 
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(a) modified Oresme sequence: 
Gin =—4"Gn = —n x 2071, 

(b) Oresme-Lucas sequence: 
B= a ok, 


(c) Oresme sequence: 


O_» = —4"O, = —n X 2”. 


6 The Sum Formula )°;_, EW saws of Generalized 


Oresme Numbers 


In this section, we present sum formulas of generalized Oresme numbers. 


Theorem 25. Let x be a real (or complex) number. For all integers m and j, for 


generalized Oresme numbers, we have the following sum formulas: 


(a) If 2-2"2? — cHm +1404 0, then 


n ‘ hs DN ah gs ze W; = DO, ae Waa 
yee Wink+j 4 : 
k=0 





1 
2-2my72 — 7 Hy, +1 On 
(b) If 2-2" 2? — cHm +1 = u(x —a)(x — b) = 0 for some u,a,b€ C with u £0 

anda#b, t.e., x =a orx=b, then 
(a(n + 2)2-2™ — (n+ 1)Hm)2"Wi4mn 


3 ky ne. + Dame: + 1)? Wiig an -_ 2 Ga 
Be ae Q-Imt1y — 
k=0 we 





(c) If 2-72? — tHm +1 = u(x — c)? = 0 for some u,c € C with u £0, ie, 
x=c, then 
(n + 1) (2-2 (n + 2)a” — na” Hm) Winn; 
ss kw. ki = + n(n + 1 ales cee, eee en 
mk+j = : 
k=0 





Q-2m+1 
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Proof. Take r=1,s = —F and H, = H,, in Theorem [2| 


Note that (6.1) can be written in the following form 


(2-2, _ Hye Wes ae De WY india 
3 Woes; = + (Hm — 2-2" 2)W; — 2-?"2Wj_m 
a nes 2-2my2 —_ x Hy, +1 : 





As special cases of m and j in the last Theorem, we obtain the following 


proposition. 
Proposition 26. For generalized Oresme numbers, we have the following sum 
formulas: 
(a) (m=1, 7 =0) 
If ; (x — 2)? #0, t.e., x #2, then 


3 a*W, _ (=) a W, +2" Wr-1 + 4Wo + 4(W, — Wo)a 
k=0 Ge 2 





and 


if 1 (a — 2)? =0, i.e. 2 =2, then 





— (n+ 1)((2 — 4)n + 22)a"—W,, + n(n + 1)a"- Whi 
Ss Ww, = 5 : 
k=0 


(b) (m= 2, j =0) 
If + (a — 4)? #0, t.e., c #4, then 


3 Wo, — 28 Wan + 2" "Wana + 16Wo + 4 (4Wi — 30) @ 
2k: F 
k=0 = 4)? 





and 


ft (a =e = 0, t.e., x =4, then 





3 2*Woe = (n+1)((@ — 8)n + 22)2"-1 Won + n(n + 1)2"-!Won_a 
5 . 
k=0 
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(c) (m=2, 7 =1) 
If + (a — 4)? #0, t.e., c #4, then 


ii k (a _ 8)" Wons1 + ttl Won_1 + 16W, + 4(W, _ Wo)x 
ae Wor+1 = (x Zz 4)? ) 





and 


fe rr aoe ay? = 0, te, 2 =4, then 





ed (n+1)((x — 8)n + 22)2"—-1 Won + n(n t+ 1)2"-1Won_1 
ye Wor = 5 : 





(a) (m=-1, 7 =0) 
If (2a — 1) #0, te, 04 5, then 


as i 4a? HAW ni +4(2 —1) a" 'W_,»z+ Wo - 401 
yews 








(Q¢ —1)? 
and 
if Qa - 1)? =0, ie, 2 = 4, then 
k n(n + 1a" 1Wenst + (n+ 1)((e = 1)n + 20)0" Wen 
ye W_k ; 


(e) (m= —2, j = 0) 
If (4a — 1)? £0, i.e, oF q, then 


3 x W. 2k = 162"*!W_2n+2 +8 (2x =, 1) oases a =p Wo —4 (4W, = Wo) Be 
- (4a — 1)? ; 





and 


if (4a — De =0,4¢e,2= =i then 


ys kwon, = 2n(n + La” !W_onye + (n + 1)((22 — 1)n + 4x)2"—-!W_on 
x La = A ‘ 
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(f) (m= —-2, j =1) 
If (42 — 1)? £0, i.e, oF q, then 





3 kw vere 162"t!W_aon+3 +8 (2a _ 1) 21 W_on41 +W,-4 (3W, = Wo) x 
k=0 = (4a — 1) 
and 


if (4a — 1)? =0, i.e, c= Z, then 





Ae 2n(n + 1)2"-!W_ony3 + (n +:1)((Qe — 1)n + 4x)2"-!W_anst 
Sox W_orsi = ri : 
k=0 


From the above proposition, we have the following corollary which gives sum 
formulas of modified Oresme numbers (take W,, = G, with Go = 0,G, = 1). 


Corollary 27. Forn > 0, modified Oresme numbers have the following properties: 


(a) (m=1, 7 =0) 
If ; (x — 2)" #0, t.€., ux #2, then 





o) 


3 CO, = (29 —4)a™t 1G, +a" 1G,_14+ 42 
b= 
ae (a — 2) 


and 


if i (a — 2)" = 0, t.e., x = 2, then 





) b> ‘ 
2 
k=0 


(b) (m= 2, j =0) 
If (a — 4)? #0, i.e, cx #4, then 





P 


3 PGs (2 — 8)a"*!Gon + 2" Gon—2 + 16x 
r= (a — 4)? 
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and 


if = (x — 4)? = 0, te, x =4, then 


” (n +1)((x — 8)n + 2x)2""!Gon + n(n + 1)2"!Gon_2 


So 2* Gor = 9 


(c) (m= 2, j = 1) 
If 3 (e- 4)° £0, i.e, « #4, then 


3 *G _ (@ — 8)a"t Gong + £"*1Gon—-1 + 40 +16 
© Gak41 = 1 5 ; 
k=0 (x — 4) 








and 
if i (a — 4)? = 0, i.e, x =4, then 


a _ (n+ 1)((a — 8)n + 22)2"“" Gong + n(n + 1a" "Gani 
Nee Gorui = 5 : 
k=0 





(d) (m=—1, j = 0) 
If (22-1)? £0, ie, «oF 5, then 


se a*G_,= 4a"*1G_ nai +4 (a - z) OG = 4 
k=0 (2x — 1) 





us n(n + 1)2"-1G_ngi + (n+ 1)\((2 — 1)n + 22)2"!G_y 


Ss 'Gs a 5 


(e) (m= —2, j = 0) 
If (42 — 1)? £0, i.e, oF q, then 





Seka = 162"*1G_onso + 8 (2x — 1) 2"*1G_on — 16x 
—2k — ’ 
= (42 — 1)? 
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and 
if (4x — 1)? =0, i.e., 2 =1, then 

— , Qn(n + 1a”! G_onse + (n+: 1)((22 — 1)n + 42)2"—-!G_on 
Dae G_ox, = A » 
k=0 


(f) (m= —2, 7 =1) 
If (4a — 1)? £0, ie, oF ;, then 


3 a*G k+l = 162" *'G_mn+3 +8 (2x — 1) t"1G_ony1 + 1-122 
= 1 ee , 
as (42 — 1)? 





and 


if (42 — 1)? =0, ie, x = = , then 





2n(n + 1)2"1G_ons3 + (n+ 1)((2z — 1)n + 42)" 41G_ong1 
Sox G_oK+1 = ml , 


Taking W, = AH, with Hp = 2,H, = 1 in the last proposition, we have the 
following corollary which presents sum formulas of Oresme-Lucas numbers. 
Corollary 28. For n > 0, Oresme-Lucas numbers have the following properties: 
(a) (m=1, j =0) 

If + (x — 2)" # 0, t.e., u #2, then 


¥ = (29 —4) a", + 21 Hy_1+8—4e 
eG (a — oy 





and 


if t(a- 2)? =0, i.e., c =2, then 





See (n+ 1)((@ — 4)n + 22)a""'H, + n(n + 1)c? Ha 
5 
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(b) (m= 2, j =0) 
If i (a — 4)? #0, t.e., x #4, then 


y a* Hox = (x — 8), + 2°"! Hono + 32 — 82 
k=0 (= 4)? 








fe age (n+1)((a — 8)n + 22)2"—-1 Ho, + n(n + 1)2"—! Hon_2 
S 2k = ; 
2 


(c) (m= 2, j= 1) 
If + (a — 4)? £0, ie, « £4, then 





3 bay, — @—8)at Hones +2" Han +16 — 40 
UT fob +1 = 5 ; 
k=0 (x —4) 





and 


if TH (a — 4)? = 0, i.e, x =4, then 





ie _ (n+ 1)((@ — 8)n + 22) 0"~" Hong + n(n t+ 1)2"~1 Aon-1 
Soa Aloxy1 = 5 : 
k=0 





(d) (m=—1, j = 0) 
If (22-1)? £0, ie, cF 5, then 


3 oH, = 4a? iy +4 (2 —1) cane +2—4¢ 
= (22 — 1) 





and 


if (2x2 — ie =0,4¢6,2= 5; then 


3 okH_, = hat Var Han + (nt I)((@— Dnt 29)? an 
,= d : 
k=0 
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(e) (m= —2, j = 0) 
If (4a — 1)? £0, i.e, oF z, then 





ri 


—— 162”41 AT onse + 8 (2a —1)a"t! Ho, +.2- 82 
~— Ho = 2 
oar (4a — 1) 


and 


if (4a — 1)? =0, i.e, c= Z> then 


- fe 2n(n + 1)2”-A_onse + (n +1)((2r — 1)n + 4x)2”-1A_on, 
a i : 





(f) (m= -2, j=1) 
If (4a —1)? £0, i.e, a4 q, then 





2 


9 a* H. 2k+1 = 16a"** Hants + 8 (2x — 1) aaa « aero +1— 4a 
= — 2 
k=0 (4a — 1) 


and 


if (4a — De =0,4¢6,2= ip then 





oa 2n(n + 1)2""1A_ons3 + (n +:1)((Q2 — 1)n + 42)" Aon 
ae: A_on41 = 1 ; 
k=0 


From the above proposition, we have the following corollary which gives sum 
formulas of Oresme numbers (take W,, = O,, with Op = 0,0, = 5). 


Corollary 29. Forn > 0, Oresme numbers have the following properties: 


(a) (m=1, 7 =0) 
If + (x — OV? # 0, t.e., x #2, then 





3 20, = (¢@ —4)a"O, + 2°t1O,_-1 + 2¢ 
: (x — 2) 


p) 
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and 
if t(e- 2)? =0, i.e., c =2, then 
(n+ 1)((2 — 4)n + 22)a"-10, + n(n + 1)2"-1On-1 


S20; = 5 


k=0 





(b) (m= 2, 7 =0) 
If + (a — Ay? #0, i.e, x #4, then 


so x" Oox = (x = 8)2"tOon a 2”*1Oon—9 + 8x 
i=, (x — 4)? 





and 


if i (a — 4)? = 0, i.e, x =4, then 


5 a8On _ (n+ 1)((x — 8)n + 22)2"-1Oon + n(n + 1)2"-!Oon_2 
k=0 











2 
(c) (m= 2, 7 =1) 

If + (a — 4)? #0, i.e, x #4, then 

3 ose (x — 8)2"*1Oon41 + 2"t1Oon—-1 + 22 + 8 
2k+1 — ; 

ar (a — 4)? 

and 

if i (x — 4)? = 0, te, x =4, then 





ae (n+1)((x — 8)n + 22)2"-1Oons1 + n(n + 1)2"-lOon_1 
Soe Or = 5 ; 
k=0 


(d) (m=—1, j = 0) 
If (22-1)? £0, ie, «F 5, then 





3 0, = Ag" HO ay +4 (a — Z glu, = ou 
a (22 — 1) 
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and 


if (2a —1)? =0, ie., e =4, then 


n(n +1)2"-1O_yii + (n +1)((e — 1)n + 22)2""!O_y 


Yo 0-4= 5 





(e) (m= —2, j = 0) 
If (42 — 1)? £0, i.e, oF 4, then 


_ 162+? O_on+2 + 8 (22 — 1) 2"*1O_on — 8a 


k 
3 ee (42 —1)? 





and 


if (42 — 1)? =0, ie,  =1, then 


2n(n + 1)2"O_anpe + (n + 1)((Qe — 1)n + 42)2"-!O_on 


n 
ye "O_o = 4 
k=0 





(f) (m= —2, j =1) 
If (4a — 1)? £0, i.e, oF q, then 





oak 32¢"t10_on43 + 16 (22 —1)2"410_on41 + 1-122 
Sox O_2k+1 = 2(4 1 2 ) 
k=0 (4a — 1) 

and 
if (4a — 1)” =0, ie., a = 4, then 


(n+ 1)2"—-!O_ons3 + (n+ 1)((22 — 1)n + 42) 2"! O_ons1 


” k an 
ye O_o¢41 = A 
k=0 


Taking x = 1 in the last three corollaries we get the following corollary. 





Corollary 30. For n > 0, modified Oresme numbers, Oresme-Lucas numbers and 


Oresme numbers have the following properties: 
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1. 

(a) Theo Ge = -3Gn + Gn_1 + 4. 

(b) Spo Gar = 5(—7Gan + Gan—2 + 16). 

(c) Typo Garg = §(—7Gan41 + Gan—1 + 20). 

(d) Sop=0 G—r = 4(G_n4i — 1). 

(e) \op—o0 G-2k = 8 (2G_2n+2 + G_on — 2). 

(f) Sopp G-2r41 = 5(16G ant+3 + 8G_an41 — 11). 
2. 

(a) peo He = -3. An + An-1 t+ 4. 

(b) Soyo Har = $(—7Hon + Hon—2 + 24). 

(c) ae, Aopy = $(—7Han41 + Hon; +12). 

(d) Yhap Hee = 2(2H_n4i — 1). 

(e) hig H-2n = 2(8H_-2n+2 + 4H_2n — 3). 

(f) Sopp Hee = §(16H_2n+3 + 8H_2n+1 — 3). 
3. 





(a) Rug Ok = —3On + On—1 + 2. 

(b) S "po Oak = 5 (—7O2n + Oon—2 + 8). 

(c) hep Ores = §(—7Oan41 + Oon—1 + 10). 
(d) Soh=0 O-k = 2(20_-n41 — 1). 

(e) hig O-2k = §(2O0-2n42 + O_2n — 1). 


(f) Soyo O-2n41 = 75 (320_2n43 + 16O_2n41 — 11). 
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7 Matrices Related with Generalized Oresme 


Numbers 


We define the square matrix A of order 2 as: 


os ie 
1 0 


such that det A = i Then, we have 


Wa 
( H) a 


and 


If we take W,, = Gy, in (7.1) we have 


Gn4i _ 1 —4 
Gh a 


We also define 


and 


Theorem 31. For all integers m,n, we have 
(a) B, = A” 
(b) CyA” = A"C, 


(c) Cram = CnBm = BmCn- 














Proof. Take r=1,s = —F in Soykan Theorem 5.1.]. 
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Corollary 32. For all integers n, we have the following formulas for the modified 


Oresme, Oresme-Lucas and Oresme numbers. 


(a) Modified Oresme Numbers. 
Mi PE \ el Gat iG | 
1 0 Gn =7 Gani 
(b) Oresme-Lucas Numbers. 
Ane 1 —; _ (n+ 1)An41 —4nHy, . 
1 0 nn, —7(n —1)Hp-1 
(c) Oresme Numbers. 
goat —\ =f 20m  -40n 
“a0. fh 8On = =20n J 
Proof. (a) It is given in Theorem 31] (a). 


(b) Note that, from Lemma [14| we have 


Using the last equation and (a), we get required result. 
(c) Note that, from Lemma |[14| we have 
Gn = 20n. 


Using the last equation and (a), we get required result. 














Theorem 33. For all integers m,n, we have 


1 
Wratm = WrGm4+1 = 4 melon (723) 














Proof. Take r=1,s = —F in Soykan Theorem 5.2.]. 
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By Lemma [IT] we know that 
(Wo — 2W1)? Gm = —4WoWinga + 4WiWm, 
so (7.3) can be written in the following form 
(Wo — 2W1)? Whim = Wr(4(Wi—Wo) Wins 1+tWoWm)+Wr—1(WoWm+1—WiWm). 
Corollary 34. For all integers m,n, we have 


1 
Gnim = GnGm+1— -<Gn-1Gm, 





4 
1 
Hy+m — AG = qiin-1Gms 
1 
Onin — OnGm+1 _ qOr-1Gm: 
and i 
Ones = 20nOm+1 — 5 On-10m- 
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Abstract 


In this paper, we give some characterizations for proper f—biharmonic 
curves in the para-Bianchi-Cartan-Vranceanu space forms with 3-dimensional 


para-Sasakian structures. 
1 Introduction 


As a natural generalization of biharmonic curves, the concept of f—biharmonic 
curves was introduced by Lu in [4]. Since this paper, many authors studied 
f-biharmonic curves in several spaces: Ou considered f—biharmonic curves 
on a generic manifold and gave a characterization for them in n—dimensional 
space forms [6]. Guvenc and Ozgur studied f—biharmonic Legendre curves 
in Sasakian space forms [2]. Karaca and Ozgur investigated f—biharmonic 
curves in Sol spaces, Cartan Vranceanu three-dimensional spaces and homogenous 
contact three-manifolds [3]. Dua and Zhang examined f—biharmonic curves in 
Lorentz—Minkowski spaces [1]. 

On the other hand, in a very recent paper [5], Lee constructed 
the para-Bianchi-Cartan-Vranceanu model with 3-dimensional para-Sasakian 
structure and found the necessary and sufficient conditions for biharmonic Frenet 


curves. 
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In this paper, we investigate f—biharmonic curves in this 3-dimensional 
para-Sasakian manifolds. We obtain some characterizations with respect to the 
special situations of curvature and torsion functions of these curves. Throughout 
the paper, all geometric objects (curves, manifolds, vector fields, functions etc.) 


are assumed to be smooth. 

2 Preliminaries 

2.1 Para-Sasakian manifolds 

We recall fundamental ingredients of para-Sasakian manifolds from [5]. A (2n + 


1)—dimensional differentiable manifold M is said to be an almost paracontact 


manifold if it admits a (1,1)-tensor field y, a vector field € and a 1-form 77 satisfying 


g=I-n@é, n(é)=1. 


For an almost paracontact manifold M, we have y€ = 0 and no yp = 0. 
If a (2n + 1)—dimensional manifold M with almost paracontact structure 


(vy, €,7) admits a compatible pseudo-Riemannian metric such that 


9(pX, PY) = —G( X,Y) + 0(X)n(Y), (2.1) 


then we say M is an almost paracontact metric manifold with the paracontact 


metric structure (y,€,7,g). Putting Y = €, we have 
(X) = g({X, &). (2.2) 
If the compatible pseudo-Riemannian metric g satisfies 
dn( X,Y) = g(X,¢Y), 


then 7 is a contact form on M, € the associated Reeb vector field, g an associated 


metric and (M, y, €,7,g) is called a paracontact metric manifold. 
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For a paracontact metric manifold M, an almost paracomplex structure J on 
M x R is defined by 


d d 


where X is a vector field on M, t the coordinate of R and f a function of M x R. 
If the almost paracomplex structure J is integrable, then the paracontact metric 


manifold M is said to be normal or para-Sasakian. 


Proposition 1. [8/ An almost paracontact metric manifold (M,y,€,n,g) is 


para-Sasakian if and only if 
(Vxy)¥ =—-g(X,Y)E+n(V)X, (2.3) 


for any vector fields X,Y on M, where V is Levi-Civita connection of g. 


2.2 Frenet-Serret equations 


Let y: 1 > M be a unit speed curve in a three-dimensional Lorentzian manifold 





M such that 7’ satisfies g(7/,7’) = €1 = +1. The constant ¢; is said to be the 
causal character of y. A unit speed curve is called spacelike or timelike if its causal 
character is 1 or -1, respectively. A unit speed curve is called a Frenet curve if 
g(y", 7") #0. A Frenet curve has an orthonormal frame field {T = y’, N, B} along 


y. Then the Frenet-Serret equations are given by 


Vrl _ cok N, 
VrN = —€,KT — €37TB, 
VrB = egTN, 


where Kk = Vy | is the geodesic curvature and 7 is the geodesic torsion of ¥y. 
The vector fields T,N and B are called tangent vector field, principal normal 
vector field and binormal vector field of 7, respectively. 

The constants €2 and €3 are defined by g(N, N) = €2 and g(B, B) = €3, and 
called second causal character and third causal character of 7, respectively. The 
equation €,€2 = —e€3 holds. 


A Frenet curve ¥ is a geodesic if and only if « = 0. 
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Proposition 2. Let {T,N,B} are orthonormal frame field in a Lorentzian 
3-manifold. Then, 


TA,N=e3B, NA, B=e1T, BALT =e0N. 


2.3 f—Biharmonic maps 


A map ¢: (Mm,g) > (Nn, h) between two pseudo-Riemannian manifolds is called 


harmonic if it is a critical point of the energy 


1 
= 5 | lao? ao, 
2 


where Q is a compact domain of M,,,. The tension field 7(@) of ¢ is defined by 


7(@) = tr(V%d@) = a0 (V2 do(e:) — db(Ve,ei)), 


w= 1 


where V® and {e;} denote the induced connection by ¢ on the bundle ¢*T'Ny. A 
map ¢ is called harmonic if its tension field vanishes. The bienergy E2(@) of the 


1 
= 5 | IrCodIP dey, 
Q 


and @¢ is called biharmonic if it is a critical point of the bienergy, where 


map ¢ is defined by 


Q is a compact domain of M,,. Clearly, all harmonic maps are biharmonic. 
Non-harmonic biharmonic maps are called proper biharmonic maps. ‘The 
bitension field T2(@) of ¢ is defined by 


= Die((VEVE, — VE.,6,)7(6) — RY (r(4), dole) dole), (2-4) 


where R™ denotes the curvature tensor of N,,. A map ¢ is called biharmonic if its 
bitension field vanishes. 
A map ¢ is called f—harmonic with a function f : M > R, if it is a critical 


point of the energy 


1 
= 5 | fllasl? dey, 
Q 
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where 2 is a compact domain of M,,. The f—tension field T¢(@) of ¢ is given by 


T(%) = fr(d) + do(gradf) 
see [7]. The f—bitension field 72,(@) of ¢ is defined by 


72,5(@) = fro(b) + Afr($) + 2V% 00/7 (0): (2.5) 


A map @¢ is called f—biharmonic if its f—bitension field vanishes (see [I], 
[4]). Non-harmonic and non-biharmonic f—biharmonic curves are called proper 
f—biharmonic curves, and if f is constant, then an f—biharmonic curve turns 


into a biharmonic curve [4]. 


3 f—Biharmonic Curves in Para-Sasakian Space 


Forms 


Lee introduced the concept of para-Bianchi-Cartan-Vranceanu model with 
3-dimensional para-Sasakian structure in |5| as follows: 
Consider the set 


D = {(e,y,2) € B®: 1+ 5 (0? +9") > 0}, 
where ¢ is a real number. Remark that if c > 0, then D is the whole R3(z, y, z). 


On the region D, the contact form 7 is taken as 


ydx — xdy 


=dz+ ; 
: 1+ $e? +9") 





Then, the characteristic vector field of 7 is € = Z. 
Next, the Lorentzian metric is equipped as 
_ —dzx? + dy? 
14 §@ +P )P 


The Lorentzian orthonormal frame field (e1,e2,e3) on (D,g_) is given by 


a. oO 
dz az 


G22 Se 


_ ee OE, ND, 
T+ 5? +9)’ 





a={1t5(e’ +x} 








= Cc, 2 by, 0 
eat s@ tye ter, 
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Then the endomorphism field y is given by 


p(e1) = €2, ple2) =e1, vle3) = 0. 


The Levi-Civita connection V of (D,g.) is described as 


Ve.e1 = —cye2, Vee2 = —cyei + €3, Ve,e3 = —€2, 
Ves@1 = —creg — €3, Ve ye2 = —cre1, Vere3 = —€1, 
Ve3€1 — —€2, Ve,€2 = 61; Ve,e3 = 0. 


The contact form 7 on D fulfills 
dn(X, Y) 7. gc X, yY), AYE x(D). 


Furthermore the structure (gc,y,€,7) is para-Sasakian. The non-vanishing 


components of the curvature tensor R of (D,g.) is given by 


Rey, €2)e2 = —{3 + C(x" — y”) ter, R(e1, €3)e3 = €1, 
R(eg,ei)er = {3+c7(a* —y’)}e2, R(e2,e3)e3 = e2, 
R(e3,e1)e1 = —e3, R(e3,e2)e2 = e3. 


For the sectional curvature K of (D,g-), we have 
K(eo, €3) =-l= —K(e3, €1), 


and 
K(e1, €2) = R(e1, €2, €1,€2) = —{38 + C(x? — y*)}. 


So, (D, gc) is of holomorphic sectional curvature H = —{3 + c?(a? — y?)}. 

For the case c = 0, the holomorphic sectional curvature H equals —3, thus the 
space D becomes para-Sasakian space form. In the next, we will deal with the 
case c= 0. 

Now, suppose that y : I > (D,gc) is a curve parametrized by arc-length 
and {T, N, B} is an orthonormal frame field tangent to D along y, where T = 
Tye, + Teeg + T3e3, N = Nye, + Noeg + N3e3 and B = Bye; + Boeo + Bge3. 

The f—biharmonicity condition for curves on (D,g.) is obtained in the 


following theorem. 
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Theorem 1. Let y: I > (D,g-) be a curve parametrized by arc-length. Then y 


is f—biharmonic if and only if the following relations are satisfied: 


Sex! f + 2x7 f' = 0, 
Kf" +26! fl + flix” + egn9 + e1K7? + Ke(€3 — 4(B)")] = 0, es 
—Qer f! — f(2n'r + K7’) — 4eikfn(N)n(B) = 0. 
Proof. Let y = y(s) be a curve parametrized by arc-length. We use formula (2.5). 
From [5], we have 
T(y) =aiVrT = —e3kN, (3:2) 
R(T, N,T, N) = €3 — 4B3, (3.3) 
R(T, N,T, B) = —4e,N3Bs3, 


T2(7) = Begnn'T + €9(K" — egn(eyK? +6377))N +1 (26'7 +7’) B+ egkR(T, NJT. 


(3.4) 
Moreover, from [i], we have 
Voraag™(9) = £'Vr(VrT) = enf'WN +x(—e1kT — est BY, og 0 
Aft(y) = f’VrT = fl eokN. 
Therefore, combining the equations (3.2), (3.4) and (3.5), we obtain 
(3.6) 
To¢(7¥) = 8egKn' fT + eof(K” — egk(€,K? + €377))N + €1f (2n'7 + xr')B 


+eofxR(T, N)T + conf” N + 2e9f' [KN + K(—e1KT — €37B)]. 


If we take inner product of equation (3.6) with T, N and B, respectively and use 


the equations (3.3), we get (3.1). 


Proposition 3. Let y: I > (D,gc) be an f—biharmonic curve parametrized by 














arc-length. If & is a non-zero constant, then y is biharmonic. 


Proof. Under the assumption « is a non-zero constant, from the first equation in 
(3.1), obviously we get f’ = 0. So, y is a biharmonic curve. 














Proposition 4. Let y: I > (D,gc-) be an f—biharmonic curve parametrized by 


arc-length. If rT is a non-zero constant and n(N)n(B) = 0, then y is biharmonic. 
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Proof. Under the assumption 7 is a non-zero constant and 7(N)n(B) = 0, using 
the first and third equations in (3.1), we get 


Km . 2f 
and oy 
K 
Th + 7? =0. (3.8) 


Putting equation (3.7) in (3.8) shows that f is constant, therefore 7 is a biharmonic 


curve. 














Proposition 5. Let y: I > (D,g-) be an f—biharmonic curve parametrized by 
F fine) 
arc-length. If 7 is a non-zero constant, then f =e T : 


Proof. Under the assumption 7 is a non-zero constant, if we use the first and third 
equations in (3.1), we obtain 


Kf! 
a ~3F (3.9) 
and 
— Qer fl — 2fr'r — dere fn(N)n(B) = 0. (3.10) 











Setting equation (3.9) in (3.10), we get the result. 


Proposition 6. Let y : I > (D,gc) be a non-geodesic curve parametrized by 





arc-length and suppose that T = 0. In this case, y is f—biharmonic if and only if 


the following equations are valid: 


Pre=d, (3.11) 
(fr)" = —fr(egn” + €9(€3 — 4n(B)”)), (3.12) 
n(N)n(B) = 0, (3.13) 


where c, € R. 


Proof. Under the assumption 7 = 0, if we use equations in (3.1) by integrating 











first equation, we deduce the results. 
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Proposition 7. Let y : I > (D,g-) be a non-geodesic curve parametrized 
by arc-length and suppose that T and & are non-constants. In this case, y is 


f—biharmonic if and only if the following equations are valid: 


Perse, (3.14) 
(fx) = —fr(egn” + e177 + €2(€3 — 4n(B)”)), (3.15) 
freer = ge (3.16) 


where c, € R. 


Proof. Under the assumption 7 and « are non-constants, if we use equations 


in (3.1) by integrating first and third equations, we obtain (3.14), (3.15) and 


(3.16). 














From the last two propositions, we can give the following theorem. 


Theorem 2. An arc-length parametrized curve y : I — (D,gc) is proper 
f—biharmonic if and only if one of the following situations is true: 


(i) 7 =0, f =cn~*/? and the curvature « solves the equation below: 


3(«’)? — 2am” = —4n? [e307 + €o(€3 — 4n(B)*)). 


f= 4eyn(N)n(B) 
te i 


(i)7 #0, = 2 , f = aK */? and the curvature & solves the 
1 





equation below: 


i _ Bele) 


3(K')? — Qen” = —4K2[e3K2(1 — €9 A ) + e2(€3 — 4n(B)?)}. 
1 





Proof. (i) The first equation of (3.1) gives 
f=an?/. (3.17) 


By replacing the above equation into (3.12), we obtain the result. 
(ii) From the first equation of (3.1), we have 


fauna (3.18) 
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Setting the above equation in (3.16), we get 
4eyn(N)n(B) 
oF _ 4e10 i n 


; 
= : 3.19 
- 2 (3.19) 


And finally putting equations (3.18) and (3.19) in (3.15), we obtain the result. 

















Consequently, we can express the following corollary. 


Corollary 1. An arc-length parametrized f—biharmonic curve y : I > (D, gc) 


with constant geodesic curvature is biharmonic. 
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Abstract 


This paper presents a new generalization of the Topp-Leone distribution called the Topp- 
Leone Weibull Distribution (TLWD). Some of the mathematical properties of the 
proposed distribution are derived, and the maximum likelihood estimation method is 
adopted in estimating the parameters of the proposed distribution. An application of the 
proposed distribution alongside with some well-known distributions belonging to the 
Topp-Leone generated family of distributions, to a real lifetime data set reveals that the 
proposed distribution exhibits more flexibility in modeling lifetime data based on some 
comparison criteria such as maximized log-likelihood, Akaike Information Criterion 
[AIC = 2k — 2log(L)], Kolmogorov-Smirnov test statistic (K—S) and Anderson 
Darling test statistic (A*) and Crammer-Von Mises test statistic (W*). 


1. Introduction 


Lifetime distributions are statistical models used for analyzing real life problems 
based on survival time. Many statistical distributions have been proposed to model 
lifetime data and the Topp-Leone distribution introduced by [14] is one of such 
distributions. In practice, most lifetime data sets encountered exhibits a bathtub hazard 
rate property and the one parameter Topp-Leone distribution happens to be the simplest 
distribution with such hazard rate property, but being a single parameter defined on a unit 
interval, its flexibility is limited in handling lifetime data sets. 


Several generalizations of the distribution have been introduced to address the 
aforementioned drawback. These generalizations are found in the works of [1-5,9,10]. In 
this paper, we introduce a new generalization of the Topp-Leone distribution which 
serves as an alternative distribution among the Topp-Leone generated family of 
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distributions. We shall call the proposed distribution, ““Topp-Leone Weibull Distribution 
(TLWD)”. 


[14] proposed a J-shaped univariate distribution with cumulative distribution function 
defined by 


F(x) = (2) (2-2). 0<x<b<w, 0<a<1 (1) 


and the corresponding density function given by 


fon = (2) (1-2)(2-2)". Q) 


The survival function and the hazard rate function of Topp-Leone distribution are 
obtained using equations (1) and (2) as 


s(x) = 1- F@) =1 -(2) (2-4) (3) 


and 


H(x) = F(x) = 2a(2)""(1-2) 


[11] improved on the mathematical properties of the Topp-Leone distribution by deriving 





(4) 


higher moments for the distribution and a general usefulness of the distribution. 
Motivations for generalizing the Topp-Leone distribution arose after the work of [11]. 
One of such motivation is that, suppose a random variable x follow the Topp-Leone 
distribution then the random variable x can have either a finite support (0 < x < b) or 


an infinite support (0 <x <b< oo). 


Suppose we fix the parameter b = 1 in equations (1) and (2), then we obtain the 
cumulative distribution function and the density function of the one parameter Topp- 
Leone distribution defined on a unit interval respectively as 


F(x) = x%(2-—x)*, 0<x<1l,a>0 (5) 
and 
f(x) = 2ax*1(1 — x) (2-—x)*1. (6) 


Figure | shows the graphical illustration of the density function and the hazard rate 
function of the Topp-Leone distribution for varying values of the shape parameter 
defined in the interval 0 < a@ < 1 and a fixed value of the scale parameter b = 1. 
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Figure 1: Density and hazard rate functions of the Topp-Leone distribution. 


The plots in Figure 1(a), clearly shows that the density function of the TLD exhibits a 
reversed-J shape for different values of the shape parameter defined in the interval 
0 <a <1 anda fixed value of the scale parameter b = 1, while the plots in Figure 1(b) 
indicates that the hazard rate function of TLD exhibits a bathtub shape whenever 
0 <a <1 and a non-decreasing shape for a => 1. The remaining sections of this paper 
are organized as follows: Section 2 presents some mathematical properties of the 
proposed distribution which include; the probability density function, cumulative 
distribution function, Hazard rate function, Survival function, Quantile function, 
Moments, Moment generating function, Renyi entropy, and distribution of ordered 
statistics. The model parameter estimation and simulation study on the maximum 
likelihood estimates of the proposed distribution are given in Section 3. Finally, in 
Section 4, we applied the proposed distribution to a real data set and compared its fit with 
the fit of some existing Topp-Leone generated family of distributions. 


2. Mathematical Properties of the Proposed Distribution 


2.1 The density and cumulative distribution functions of the proposed distribution 


[13] introduced a new class of the Topp-Leone generated family of distributions with 
the cumulative distribution function defined as 


F(x) = [G@)]*[2 - G@)]*, (7) 
and the corresponding density function given by 
F(x) = 2ag(x)(1 — G(x) )[G@)]** [2 - Gx). (8) 
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Taking G(x) as the distribution function of the Exponentiated exponential 
distribution developed by [8], the authors obtained the cumulative distribution function of 
the Topp-Leone Generalized Exponential Distribution as 


a 
F(x) = [1-e*]/"[2- (1-e*)] , (9) 
and the corresponding density function given by 


f(x) = 2apae*[1 — oP (1 - (1 -e-*)*) 


x[2-(a-eF (10) 


Using a similar approach of generating new distributions defined in equations (7) and 
(8), we assume the random variable X to follow the 2-parameter Weibull distribution with 
cumulative distribution function and probability density function respectively defined by 
G(x) = 1-78" (11) 
and 
g(x) = aOx%1e- 9x", (12) 


Inserting equations (11) and (12) into (7) and (8), we define the cumulative 
distribution function and the probability distribution function of the Topp Leone Weibull 
distribution (TLWD) respectively as 


Fx) ={-e*)[(2-(- e7Oxy 17 
={1- en 20xey 6) 
and 
fla) = Pada ter" —(1—e [ea ae 
= 2nada®en2"(1 = 0-2) We 
The series representation of equation (14) can be obtain as 


wx a= = 5 so, 
(1—e-26x*)"" = ey ae ‘) (—1)) e202) 
A-1 


f@) = 2220 Y-( ; \ Opiate teen, (15) 
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The graphical plots of the probability density function are shown in Figure 2. 





f[x) 








Figure 2: Density function of the TLWD for fixed value of the parameters. 
The plots in Figure 2 clearly indicate that the density function of the TLWD exhibits 
a left-skewed, right-skewed, reversed J-shape and symmetric unimodal shapes. 
2.2 Survival function and hazard function of the proposed distribution 


Using equations (13) and (14), the survival function and the hazard rate function of 


the Topp Leone Weibull distribution are given by 
S(x) = 1- F(x) 





=1— {t= en 2020) (16) 
and 
— FS 
H(x) = 1-F(x) 
(17) 


2AaO xt 1e-20x% 


- {(a-e-20x%) (1-7262")| 
The graphical plots of the hazard rate function of the TLWD for different choice of 


the parameter value is shown in Figure 3. 
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h(x) 
4 





Figure 3: Hazard rate function of the TLWD for fixed value of the parameters. 


2.3 Quantile function of the proposed distribution 


The quantile function of a random variable x is obtained by solving the system 
equation of F(x) = p. Thus, given the cumulative distribution function F(x) defined in 


equation (13), the p“” quantile function of the TLWD can be obtain as 
i= en 20xey = p 
1 —e- 26x" = pla 
Pe ae pla 
taking the natural logarithm of both sides we have 
—20x* = log (1 = p/2) 
x* = —s,log (1 a pla) 
x ={-Slog (1- pay", O<p<1. (18) 


The median of the TLWD can be obtained by substituting p= A in equation (18) 
which yields, 


median = {-—log (1- 05)"2)} (19) 
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2.4 Moments of the proposed distribution 


Let X be a continuous random variable with density function f(x), then the r‘” 
moment about the origin of X is defined by 


E[X"] = f°..x" f)dx (20) 


inserting the density function of the TLWD into equation (20), the r“” moment about the 


origin of the random variable X is defined by 


co a a A- 
E[X"] = f° 2daOxT*4-1e-26x7(4 — g26x" 1" gy (21) 


aA- = ; Bs 
{1 — e-20x7V"" = ey ) (—1)) e262) 


so that equation (21) now becomes, 
A-1 


ELX"] = 2208 S720(" 


) (- 1)/ i xr ta-1 p-20x% (14) dy (22) 


evaluating the integral part of equation (22), we have 


co i a a 7 
is xtta-1, 20x At) dx 


/ 
let y = 20x%(1 + j), which implies that, x = aos i 


VWq-1 


oe tf ay 
so that, dx = aap Pacem y 


1/ rta-1 1/ -41 
co y a = 1 y a 
Jy {lee ee Fecem a 


ib ylati-YVatVa-1 -~yq 
0 af2e(atp]/att-VatV/a yey 


co ylaeY PCV # 1) 
fo a[20(1+j)]"/a*4 ey a[2o(1+j)]'/ats * 


Thus, the r“" moment about the origin of the TLWD is given by 


z ae fae ;  TC/at1) 
E[X"] = 218 Di20( j ) (1) ee. (23) 


The first four raw moments of the TLWD are obtained from equation (23) as 
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220 B-0(" 5°) (1 pean 


Jj 20(1+/)) /att’ 


= 
ie 
ll 


r(2/q+1) 
20(1+ sy /att’ 


i ‘) (—1yi —TCla*t) 


zr 
os 
| 
N 
SS 
D 
M 

mm 18 
Oo 

"oN Fi aa. 


i 
w 
II 
N 
SS 
D 
M 
78 
oO 


20(14f)|/att’ 


r(A/at) 
20(14f)|*/att * 





Furthermore, using the first four raw moments defined above, the measures of 
skewness (S;) and kurtosis (K,) are obtained as 


— H3-3ugMt2u3 ei | Wee H4a—4U3Ut+6 MyM? —3y* 
: (uh-n2)'/2 . (Ha —H?)? 


Tables 1 and 2 show the theoretical moments of the TLWD for different values of the 
parameters. 


Table 1: Theoretical moments of TLWD for fixed value of the parameter (a = 2). 








ul. (0=2,A=3) (0=2A=5) (0=4A=3) (6=4,A=5) 
nA 0.6452 0.7310 0.4562 0.5169 
us 0.4583 0.5708 0.2292 0.2854 
ii 0.3542 0.4741 0.1252 0.1676 
a 0.2951 0.4173 0.0738 0.1043 
o 0.0420 0.0364 0.0211 0.0182 
Sy 0.4980 0.5143 0.4581 0.5078 
K, 3.2886 3.3859 3.5931 3.4248 





Table 2: Theoretical moments of TLWD for fixed value of the parameter (a = 5). 








ul. (0=2,A2=3) (0=2,A=5) (6=4A=3) (0=4,4=5) 
ul 0.8288 0.8750 0.7215 0.7617 
ne 0.6986 0.7740 0.5294 0.5866 
uh 0.5982 0.6921 0.3946 0.4566 
i 0.5199 0.6254 0.2986 0.3592 
o? 0.0117 0.0084 0.0088 0.0064 
Sx -0.1406 0.2517 -0.1393 0.0333 
K, 3.3840 0.9577 4.4313 5.0395 
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From Tables | and 2, we observed that the TLWD can be right skewed (S; > 0), left 
skewed (S; <0) and approximately symmetric (S; ~ 0). Also, at some fixed values of 
the parameters, the distribution can be leptokurtic (K, > 3), platykurtic (K, < 3) as well 
as mesokurtic (K, ~ 3). This claim was clearly shown in Figure 2. 


2.5 Moment generating function of the proposed distribution 


Let X be a continuous random variable with density function f(x), then the moment 
generating function of X is defined by 


co 


Mx(t) = f_,e*f(x)dx (24) 


inserting the density function of the TLWD into equation (24), we obtain the moment 
generating function of the TLWD as 


(oe) a ayA-1 
My) = J, e2Aaex® e778 (1 = e028") dx (25) 
using the series representation, 


wax = . ae 
{1 — en26x"}"" = ye Cs :) (—1)/ e~ 20% 


(tx)™ 


m! 





eft = peters 


equation (25) becomes, 
My(t) = 2200 S729 Siao( 5 (CII ame te 20 dx (26) 


evaluating the integral part of equation (26), we have 


co a. ay, a 7 
J xmta-1, 20x +N) dx 


co y Va oe -y 1 y yg 
mer 20(14)) C * 204) L264) y 


= ic yat1-YatVa-1 -Yq 
~ 40 a[20 (+f) att Vat Ta © y 











= co ylae-Y 
~ Jo a[2e(1+j)]"/att dy 


oO mt+a-1,—20x*(14+)) = PCat 
(ee e dx azo att (27) 


upon substituting equation (27) into equation (26), the moment generating function of the 
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TLWD is given by 


Mx(t) = 228 Seo Eeao(’ ‘et ee (28) 


m! [20(1+j)]"/atl 


2.6 Renyi entropy of the proposed distribution 


[12] defined an entropy of a random variable X as a measure of variation of 
uncertainty associated with the random variable X. The Renyi entropy of X with density 


function f(x), is defined by, 
tRY) = qalog ff’ @)dx, y>0, y#1 (29) 


inserting the density function of the TLWD into equation (29), we obtain the Renyi 
entropy of the random variable X following the TLWD as 


T(V) = sylog|(2aay” f° x@-Dy 9 28x* Y{1- e202 MY Gy] (30) 


using series representation, 


f1— en20xtyt? S(, 4) (—1)) e-20x*%/ (31) 


substituting equation (31) into equation (30), we have 
as Y (a-1)y p-20x%(y+j) 
= log (2Aa0)” y\F-0 (- 1) f x Ye Y*Ddx| (32) 
again, evaluating the integral part of equation (32), 
{ x (@-DYV 9 -20x"(rt)) dy 


VWa-1 


90 y Ha\ 1 y 
== 7 Pree s 
Jo {Leseaal C 0+) cna AY 


yvVatYa-1 


al2a(y+p~ "lat Va 


e Ydy 


=" 


(oe) yY + AMy md ey 


0 ary, AY 


a[2e(1+j)]”t 


(1-y) 
[ x@-Dr e264 0 dy = te a). (33) 
: alzocy+syyyt a 


upon substituting equation (33) into equation (32), the Renyi entropy of the random 
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variable X following the TLWD is given by 


r(y + S “Y/q) 


ee ES 4 
ala(y+ pt Va 


ta) = zalog | eraoy Ey0(“7 7”) -1y 


2.7 The distribution of the ordered statistics of the proposed distribution 


Suppose that Yy.. < Yon <-+*+ <Yp.n is the order statistics of a random sample 
generated from TLWD, then the probability density function of the k“” order statistics, 
say X = Y,.n is given by 


h(x) = Sapa FOAL - FOF) (35) 


substituting the cumulative distribution function and the density function of TLWD 
defined in equations (13) and (14), into equation (35), we have 


2AaOn! oo (oe) n—k A-1 we a . 
h(x) = ee | )( j )Cnite 20x%(1+j) 


itk-1 


a A 
x x1 [fa — @ 20x") (36) 
using the series representation, 


panty Mitk-1) reo i+k— peek 
{1 — e7 26x") eR) Dian aa << D) (—1)™e~20x%m 


so that equation (36) becomes, 


eae ties pe 5a @ : 2) a ') ae A - 1)) 


x (—1)!ti tm @-1 g—20x%(1+j+m) (37) 


Hence, the s‘” moment of the k“" order statistics from the TLWD is defined by 


E(X§) = J, xShi(x)dx (38) 
_ __ 24a8n! nw po po (Nn—k\(A-1\(AG+k—-1) 
iz Gabe i= dj=0 Zieo( i ) ( r ) ( m ) 
x (—1)'ti+m [ce yxsta-1_-20x(1+j+m) (39) 


evaluating the integral part of (39), 


i xSta-1_-20x%(14+j+m) 
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lety = 20x%(1+j +m), which implies that, x = Pocesesl 
Vqy-1 


dy 


1 


so that, dx = a20(1+j+m) Pacers 


co 1/ Sta-1 1/ -1 
Cllactaal | aatamlatenl 
0 20(1+j+m) a20(1+j+m) L20(1+j+m) 


ee ylatt-Vatta-1 
0 @[2e(itj+m))/at1-/at*/a 
ie ylae-Y a = T(S/q@+1) 
0 @[20(1tjampart OY a[20(14+j+m)]/att 


e dy 


equation (39) now becomes 


[=a ealeu gp a 


U m 


_4yitjtm __PC/at1) 
os a [20(1+j+m)]}/a*? * (40) 


3.0 Parameter Estimation of the Proposed Distribution 


3.1 Maximum likelihood estimation 
Let (%1,%2,..,X,) be random samples from the TLWD, then the log-likelihood 
function of the TLWD is defined by 
£(x,6) = XL, log [2aaox%te-262°f4 — eee | ob = (a,6,Aa)* 
= nlog(2Aa8) + (a — 1) EP, log(x;) — 20D, xF 
+(A—1) X%, log(1 — e729"). (41) 


On differentiating the log-likelihood function with respect to the parameters, we obtain 
the score function as, 


20x& 


l i i 
So = B+ DhLy log(x;) — 26 Diy xf log (xi) + 202 - 1) EE ee 


= = <+ YL, log(1 — 729%: ) 
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ae on Vise. Bde we xe 20x 
30 a 2 a1 % + 2(A - Ndi feat) 


The maximum likelihood estimator ¢ of @ can be obtained by solving the system of 


oe 


equation ae 0. This equation cannot be solved using direct method since the system of 


equation is non-linear. However, an iterative method of solution such as the Newton 
Raphson method can be used. The Newton Raphson iterative scheme has the form: 


b= by —H*(bx)U(bx), & = (a,0,a)" (42) 


where U(@,) is the score function and H(@;) is the Hessian matrix, which is the second 
derivative of the log-likelihood function. The Hessian matrix is defined by 


aze 2k ate 
dada dadd dAadA 
aze 2k ate 
H = 
(Px) a00a 0000 a0aA 
aze ae 2B 
adda OAdO DADA 









































where, 

-20x% a .-20x% 
ae on ae 2 n xillog(xi)le"20* (1-20xf-27 7°) 
“~=-4 20 7, x7 [log (x) ]° + 208A —-1) oe, =e ety’ 

-20x% a .-20x% 
ae R n yd n “loge U(1-20xf-e ') 
9a00 ~ ‘d0aq ~~ 2 di=1%t EEO cee —— ieee = 
a2e 7K n xpe 20% 
a00A  aaae 3 ye 2887) 

GF se Ohta ae xflog(xije 2 
dada sé aS i {1-e-26*7} 
ae nm 
dada AP 

a 
ae on n pee 7% 
aoag oz 4A Mier e aety 


3.2 Interval estimate 


The asymptotic confidence intervals (CIs) for the parameters of TLWD (a, 6,A) are 
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obtained according to the asymptotic distribution of the maximum likelihood estimates of 
the parameters. 


Suppose @ = (@,6,A) is MLE of ¢, then the estimators are approximately bi-variate 


normal with mean (a, 6, A) and the Fisher information matrix is given by: 


I(x) = —E(H(4x)). (43) 


The approximate (1 — 6)100 CIs for the parameters a, 8 and A respectively, are 


at Z5),yvar(@) Ot 25), [var(6) and A + 25), var(A) 


where var(@), var(@) and var(A) are the variance of @,@ and A which are given by the 
diagonal elements of the variance-covariance matrix 1~1(,) and Zs D is the upper (9/ 2) 


percentile of the standard normal distribution. 


3.3 Simulation Study 


In this section, we investigate the asymptotic behaviour of the maximum likelihood 
estimates of the parameters of the Topp-Leone Weibull distribution (TLWD) through a 
simulation study. A Monte Carlo simulation study is repeated 1000 times for different 
sample sizes n = 50, 75, 100, 200 and parameter values (a2 = 1.0,6 = 0.3,A = 0.5), 
(a = 1.0,6 = 0.3,A = 0.5) and (a2 = 1.0,0 = 0.3,A = 0.5). [6] suggested an algorithm 
for the simulation study employing the following steps. 


1. Choose the value N (i.e. number of Monte Carlo simulation); 


2. choose the values @p = (@, 99, Ag) corresponding to the parameters of the 
TLWD (a, 6, A); 


3. generate a sample of size n from TLWD; 

4. compute the maximum likelihood estimates ¢, of @,; 
5. repeat steps (3-4), N-times; 

6. compute the: 

(i) Average Bias = va, = bo), 


(ii) Mean Square Error (MSE) = —vialh, a bo) and 
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(iii) Coverage Probability of the 95% confidence interval of the estimates @;given by 


N 
1 a R a Ae 
CP(o) = ADI (¢- Z5),,\var($) <do< bit Zep, var()) 
t=1 





where /(-) is an indicator function and | var(¢) is the standard error of the estimate ¢;. 


The coverage probability computes the proportion of times the confidence interval 


contains the true value of the parameter do. 


Table 3: Monte Carlo simulation results for average bias of the MLE. 














Parameter N Average Bias (a) Average Bias (0) Average Bias (1) 
a= 1.0 50 0.2970 0.0314 0.0994 
0=0.3 75 0.1357 0.0293 0.0724 
A=0.5 100 0.1364 0.0082 0.0362 

200 0.0507 0.0068 0.0157 
a=1.5 50 0.4154 0.0341 0.1622 
0=0.4 75 0.2859 0.0141 0.0740 
A= 0.6 100 0.1714 0.0140 0.0540 
200 0.0721 0.0074 0.0242 
a= 2.0 50 0.5309 0.0403 0.2115 
0=0.5 75 0.2740 0.0087 0.0883 
A=0.8 100 0.2340 -0.0114 0.0319 
200 0.0735 0.0104 0.0513 








Table 4: Monte Carlo simulation results for mean square error (MSE) of the MLE. 








Parameter N MSE (a) MSE (6) MSE(A) 
a= 1.0 50 0.6165 0.0727 0.2892 
A0=0.3 75 0.2463 0.0422 0.1285 
A=0.5 100 0.2105 0.0398 0.0916 

200 0.0598 0.0146 0.0271 
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a=1.5 50 1.5398 0.0995 0.4345 
0=0.4 75 0.8736 0.0511 0.1777 
A=0.6 100 0.3915 0.0392 0.3917 

200 0.1208 0.0197 0.0477 
a= 2.0 50 3.4213 0.0975 0.6521 
0=0.5 75 0.9752 0.0568 0.2400 
A=0.8 100 0.5449 0.0408 0.1602 
200 0.2143 0.0260 0.1023 

Table 5: Monte Carlo simulation results for coverage probability of the 95% CIs of the 

MLE. 

Parameter n CP (a) CP (0) CP (A) 
a= 1.0 50 0.8200 0.8133 0.8333 
0=0.3 75 0.8800 0.8933 0.8933 
A=0.5 100 0.9200 0.8867 0.8967 

200 0.9300 0.9400 0.9467 
a=1.5 50 0.8667 0.8000 0.8167 
0= 0.4 75 0.9200 0.8967 0.8867 
A=0.6 100 0.9433 0.9233 0.9033 
200 0.9500 0.9200 0.9133 
a= 2.0 50 0.9167 0.8900 0.8700 
0=0.5 75 0.9567 0.9133 0.8900 
A=0.8 100 0.9700 0.9067 0.8800 
200 0.9300 0.9167 0.9233 








Tables 3, 4 and 5 present the Monte Carlo simulation results for the average bias, 


mean square error and coverage probability of the 95% confidence interval of the 


parameter estimates of the TLWD at different choice of parameter values. Clearly from 


Table 3, we observe that the parameter a is positively biased, parameter / is positively 


biased, while parameter 0 could either be negatively or positively biased. Also from 


Table 4, as the sample size n increases, the values of the mean square error of the 


parameter estimates decreases (tends to zero). Finally, from Table 5, we observe that the 


coverage probabilities of the CIs are quite close to the nominal level of 95%. 
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4. Application of the Proposed Distribution 


In this section, we fit the proposed distribution to a real data set alongside with some 
well-known lifetime distributions from the Topp-Leone Generated family of distributions 
with density functions given by; 

(1) Topp-Leone Inverse Weibull Distribution (TLIWD) 


a-1 


ed oe 


(11) Topp-Leone Bur XII Distribution (TLBXIID) 





= -2q 14-1 
fa) = 2a0x9 (1 4x9) OPP [1 — (24 x9] 
(111) Topp-Leone Exponential Distribution (TLED) 
f(x) = 2ade~*4*(1 _ en 2ax)e™*. 
(iv) Topp-Leone Linear Exponential Distribution (TLLED) 
Ax? Ax? Gee 
f(x) = 2a(6+ ane 20+) f a ee) : 
(v) Topp-Leone Nadarajah-Haghighi Distribution (TLNHD) 
f(x) = 2aOA(A + Ox) 4e2(4-4460") [4 — erase) 


The dataset consists of the tensile strength, measured in GPa, of 69 carbon fibers 
tested under tension at gauge lengths of 20mm reported in [7]. Table 6 presents the 





dataset. 
Table 6: Tensile strength, measured in GPa, of 69 carbon fibers. 
1.312 1.314 1.479 1.552 1.700 1.803 1.861 
1.865 1.944 1.958 1.966 1.997 2.006 2.021 
2.027 2.055 2.063 2.098 2.14 2.179 2.224 


2.240 2.253 2.270 2.272 2.274 2.301 2.301 
2.359 2.382 2.382 2.426 2.434 2.435 2.478 
2.490 2.511 2.514 2.535 2.554 2.566 2.57 
2.586 2.629 2.633 2.642 2.648 2.684 2.697 
2.726 2.770 2.773 2.800 2.809 2.818 2.821 
2.848 2.88 2.954 3.012 3.067 3.084 3.090 
3.096 3.128 3.233 3.433 3.585 3.585 
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The comparison criteria considered in this paper includes, the estimates of the 
parameters of the distribution, Log — lik, Akaike Information Criterion AIC = 2k — 
2log (L)], Kolmogorov-Smirnov Statistic (K —S) and Anderson Darling test statistic 
(A*) and Crammer-von Mises test statistic (W*). Where, k is the number of estimated 
parameters and L is the value of the log-likelihood function evaluated at the parameter 
estimates. 


Table 7: Summary statistics for the tensile strength of 69 carbon fibers dataset. 





Zk 


Models Estimates Log-Lik AIC K-S Ww A 


* 


(p-value) (p-value) (p-value) 





ae 0.0395 0.0154 0.1435 
TLWD 6 = 0.0194 —~48.8598 103.7195 ; ; 

eee: (0.9999) (0.9996) (0.9990) 

C= E10 0.1216 03164 2.1117 
TLBXID 6@=1.8205  —59.9967 125.9933 

Saar ok e (0.2591) (0.1215) (0.0799) 

& = 0.5469 0.1176 0.2617 1.7344 
TLIWD 6=34.8899  —58.0304 122.0608 

pete (0.2960) (0.1741) (0.1294) 

= 88.1467 0.0950 0.1601 ‘1.1210 

TLED = 88.1467 46201 113.2403 

A = 1.0186 (0.5617) (0.3608) (0.2994) 

cece 0.0660 (0.0627 0.4406 
TLLED  6=-0.0826 —50.5881 107.1762 

eee (0.9240) (0.7985) (0.8072) 

Gee 0.0652 0.0604 ~—-0.4235 
TLNHD 6=0.0243 -50.4692 106.9383 

Bearer (0.9311) (0.8129) (0.8246) 





Discussion of Result 


The superiority of a model over another can be characterized by the model having the 
maximized loglikelihood and the least value in terms of AJC, K-S Statistic, Anderson 
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Darling test statistic (A*) and Crammer-von Mises test statistic (W*). Table 7 clearly 
reveal that the proposed distribution (TPGLD) having the maximized loglikelihood and 
the least values in terms of AJC, K-S Statistic, Anderson Darling test statistic (A*) and 
Crammer-von Mises test statistic (W*), demonstrates superiority over the Topp Leone 
Inverse Weibull Distribution (TLIWD), Topp Leone Bur XII Distribution (TLBXIID), 
Topp Leone Exponential Distribution (TLED), Topp Leone Linear Exponential 
Distribution (TLLED) and Topp Leone Nadarajah-Haghighi Distribution (TLNHD) in 
modeling the lifetime dataset under study. This claim was further supported by graphical 
illustration of the density and cumulative distribution fit of the distributions for the real 
lifetime dataset as displayed in Figure 4. 



















































Density fit for the Data Set Cumulative distribution fit for the Data Set 
(va) og 
o 7 — TLWD = 
---- TLBXIID 
=. TLIWD: 2 | 
(7 eam cL Ua TLED 2 
ad --- TLLED 
_— TLNHD © | 
P= ue 
7 
5 37 a 
a on 
N os 
o | way N 
ra Yeh ° é 
7s sae . --- TLLED 
o | o | ~~ TLNHD 
o oO 
I T T T 1 T T T T T 
15 2.0 25 3.0 35 15 2.0 25 3.0 35 
data data 


Figure 4: Density and cumulative distribution fit of the distributions for the dataset. 


Conclusion 


In this paper, we proposed a new member of the Topp-Leone generated family of 
distributions called the Topp-Leone Weibull Distribution (TLWD). The mathematical 
properties of the proposed distribution such as; the density function, cumulative 
distribution function, survival function, hazard rate function, moments and related 
measure, moment generating function, quantile function, Renyi entropy and the 
distribution of ordered statistics were derived. The method of maximum likelihood 
estimation was used in estimating the parameters of the proposed distribution. 


Finally, an application of the proposed distribution to a real lifetime data set 
alongside with the Topp-Leone Bur XII distribution, Topp-Leone Inverse Weibull 
distribution, Topp-Leone Exponential distribution, Topp-Leone Linear Exponential 
distribution, and Topp-Leone Nadarajah-Haghighi distribution, reveals that the proposed 
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Topp-Leone Weibull distribution, demonstrates superiority and offers more flexibility in 


modeling the lifetime data set under study. 
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Abstract 


The aim of this article is to initiating an exploration of the properties of bi-univalent 
functions related to Gegenbauer polynomials. To do so, we introduce a new families 
Ts, ¢,4,0,9,3,t,6) and S;(o,n, 8,3, t,5) of holomorphic and bi-univalent functions. 
We derive estimates on the initial coefficients and solve the Fekete-Szeg6 problem of 
functions in these families. 


1. Introduction 


“In [20] Legendre studied orthogonal polynomials comprehensively. The importance 
of orthogonal polynomials for contemporary mathematics as well as for a wide range of 
their applications in physics and engineering, is beyond any doubt. It is well-known that 
these polynomials play an essential role in problems of the approximation theory. They 
occur in the theory of differential and integral equations as well as in mathematical 
statistics. Their applications in quantum mechanics, scattering theory, automatic control, 
signal analysis, and axially symmetric potential theory are also known [7,12]. In 
practical, the Gegenbauer polynomials is special case of orthogonal polynomials. They 
are representatively related with typically real functions Tp as discovered in [19]. 
Typically, real functions play an important role in the geometric function theory because 
of the relation Tp = COSp and its role of estimating coefficient bounds, where Sp 
indicates the family of univalent functions in the unit disk with real coefficients and COSp 


denotes the closed convex hull of Sp.” 
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On this subject in geometric function theory, the so-called Fekete-Szegé type 
inequalities (or problems) which estimate some upper bounds for |as - yas | for 
holomorphic univalent functions. Its origin was in the disproof by Fekete and Szeg6 [16] 
conjecture of Littlewood and Paley that the coefficients of odd univalent functions are 
bounded by unity. 


Let A stand for the collection of functions f have the type: 


f@=z+ > Gans (1.1) 
n=2 


which are holomorphic in the open unit disk U = {z € C : |z| < 1}. 


Further, symbolized by S the subfamily of A consisting the functions that are 
univalent in U. 


According to the “Koebe One-Quarter Theorem” [13] each function f from S has an 


inverse f~+, which fulfills 


f'F@)=z (ev) 


and 


1 
fF) =m, (wl <ro(f), m6) =F) 
where 
g(w) = f-1(w) = w — aqw? + (2a% — a3)w? — (5a? — 5Sana3 +a4)w* +. (1.2) 


A function f € A is named bi-univalent in U if together f and f~+ are univalent in 
U. Let 2 indicate the family of bi-univalent functions in U satisfying (1.1). Beginning 
with Srivastava et al. pioneering work [36] on the subject, the large number of works 
associated with the subject have been presented (see, for example [1,2,4,5,8,9,10,11,14, 
17,18,21,22,25,28,29,30,31,32,33,34,35,37,38,39,40,41]). We see that the set 2’ is not 
empty. We see that the functions 

Z 1 14+2z 
aay =log(—) and —log(1 —z) 


are in the set X with the corresponding inverse functions 
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respectively. But the functions 
z? Zz 
LZ 2 and eg ‘ 
are not a member of the set 2’. 
The problem to find the bound of |a,|,(m = 3,4,...) of functions f € 2 is still an 
open problem. 


The fundamental distributions like, the Pascal, the Binomial, the Poisson, the 
Logarithmic, the Borel have been partially considered in the “Geometric Function 
Theory” from a theoretical point of view (see for example [6,15,24,26,43]). 


We say that the discrete random variable x have a beta negative binomial distribution, 


B48) 7 BM+6,A41) 1606 if 


if it has the values 0,1,2,3,... with the probabilities BGA)? BG) ?2 


1) aa .. , respectively, where 7, 6,4 are named the parameters. 
a ee ais See _T@+Drq+ Ort orm +a) 
aaa i ae a Bay) ttT@)lm+otat+0r@ra) 


(err 


ee ET ae 
(n+)9(0+n+2,1!’ 


where (@),, is the Pochhammer symbol defined by 


@,_= = Sone 


T(a) = la(at1)...(a+n-1) (nEN). 


Recently, Wanas and Al-Ziadi [42] studied the following power series whose 
coefficients are probabilities of the beta negative binomial distribution: 


(1) 6(®)n-1 AQ) n-1 


£2 (@)a2+) eee a, ZeU, 
= 419 + oO ++ Yn-1M— V! 


where 7,3,0 > 0. We see that, by making use of ratio test we deduce that the radius of 
convergence of the above power series is infinity. 


Now, we consider the linear operator BF : A — A which is defined as follows: 


zEU, 


= ep ae) 
Bf @) = Ba(@) « f@) = 24 YO nt aga, 
n=2 





Earthline J. Math. Sci. Vol. 7 No. 2 (2021), 403-427 


406 Abbas Kareem Wanas 





6699 
* 


where indicate the convolution of two series. 


“For the functions f and g be holomorphic in U. We say that the function f is said to 
be subordinate to g, if there exists a Schwarz function w holomorphic in U with w(0) = 
0 and |w(z)| < 1 (z € U) such that f(z) = g(w(z)). This subordination is indicated by 
f<g or f(Z) < g(Z)(zE VU). It is well known that (see [23]), if the function g is 
univalent in U, then f < g if and only if f(0) = g(O) and f(U) c g(U).” 

Recently, Amourah [3] studied the generating function of Gegenbauer polynomials 


H5(z,t) that are given by the following recurrence relation: 


1 


Net) (1 — 2tz + 22)9’ 


where 6 € R \ {0}, t € [—1,1] and z € U. For fixed t, the function Hs is holomorphic in 
U, so it may be expanded in a Taylor-Maclaurin series as note that if t = cos 6, where 


B €(—4,5), then 


co 


a 1 = 6 n 
Hs(z,t) = (i 2tz +225 ~ 2% (t)z", 


where G.°(t) is Gegenbauer polynomial of degree n. 


Clearly, Hs generates nothing when 6 = 0. Thus, the generating function of the 
Gegenbauer polynomial is set to be 


Ho(z,t) = 1—log(1 — 2tz+ 27) = y Ge(t)z". 
n=0 


oe ; rane 1s 
Furthermore, it is worth to mention that a normalization of 6 to be greater than — 5 is 


desirable [12,27]. Also, Gegenbauer polynomials can be introduced by the following 


recurrence relations: 


1 
Gn) = 5[2t(n + 6 — DGn-1() — (n + 26 —2)Gn-1 0], 
with the initial values 
Got) =1, Ge(t) = 26t and G8(t) = 26(6 + 1)t? -6. (1.3) 


Remark 1.1. Choosing the special values of 5, the Gegenbauer polynomial G°(t) 
reduces to the following well-known polynomials: 





http://www.earthlinepublishers.com 


New Families of Bi-Univalent Functions Governed by Gegenbauer Polynomials 407 





1) Taking 6 = 1, we have the Chebyshev Polynomials. 


2) Taking 6 = > we obtain the Legendre Polynomials. 


2. Main Results 

This section start with defining the families Ts(v,¢,u,,7,0,31,t,6) and 
Ss(o,7,9,4,t, 6) as follows: 

Definition 2.1. ForO <y<1,0<¢@<1,0<syw<1,teE G. 1] and 6 is a nonzero 


real constant, a function f € 2 is called in the family T;(y, ¢, u, , 0,3,t, 6) if it fulfills 
the subordinations: 


’ ’ ny 4 
( (88,f@) ) [i _ pi eiv@) , ( 28is@) } 1 





< 
(1 - 2tz + 22) 


Bf) Brf@ (B°,f@) 
and 
YY 1, mW gp 
Be Be Be 
Ww ( ing) (1 = x) Ww ( aaa (w)) 4 ii 1 i Ww ( nr ) 2 1 — 
B,39(w) Brg (w) (8°, 7 (w)) (1 — 2tw + w?) 


where the function g = f~? is given by (1.2). 


fo: 1 ; ; 
Definition 2.2. ForO <o0 < 1,t € G. 1 and 6 is a nonzero real constant, a function 


f €% is called in the family S;(o,7, 6,4, t, 6) if it fulfills the subordinations: 


(Bi.7@) + oz? (B8,F@). (1 — 2tz + z2)6 
and 
Wy (88,9(w)) + ow? (88,90) (1 —2tw +w2)8’ 


where the function g = f~? is given by (1.2). 


In particular, if we choose ¢=0 and y=1 in Definition 2.1, the family 
Ts(v,¢,4,n,9,3,t,6) reduces to the family Ts(y,0,3,t,5) of bi-starlike functions 
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which fulfills the conditions: 
2(Bif@) _ 1 
Bi Ff) (1 — 2tz + 27)? 
and 
w (84 ,9(w)) 1 
0 7 2)6’ 
BrgW) (1 —- 2tw+w?) 
where the function g = f~? is given by (1.2). 


If we choose o = 0 in Definition 2.2, the family S;(0,7,0,4,t,6) reduces to the 
family §»5(y, 6,4, t, 6) of bi-convex functions which fulfills the conditions: 


B Eve). a eee 
(wif) | 2a +2) 
and 
(B%,g0n) A 2a Wye 
where the function g = fo ON 


Theorem 2.1. ForO<y<1,0<¢<10<u<1te G1] and 6 is a nonzero 


real constant, let f € A be in the family Ts(y, ¢, yn, 8,4, t, 6). Then 


2/d|tl +4444 ITO) 2/dlt 





la2| < 
502(y + o(u + 1)) F2 + O20 + 1) 2(n +3) 
5+ DO(v+ou+)) POG+Orat+yra+y) |p. 
-06707,¢,47,6,)' (7+ Ory +Il2H+ 6414+ II(H)ra) 
and 


eae 4l?(n +9444 1)F72()F72(0)62t? 
pe A2(y+ out 1) T2@y + 9)T2(4 + 1)F2(n +3) 
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Fe 2°(n+6+3+ 2) ~M)TQ)ld|t 
66 +1)(yt+¢(2u+ DY) + OFAt 2M +3)’ 


where 


2(6 + 1)(y + O(2u + 1))TQ + 2) 


a 6,1) = 
(v,¢,47, 9,4) lq +0 +342) 


ara +OI704 DF +D[yv -D+oH4+ D2y + - D+ D)-2(y + oBu+ D)] 
l’mM7+04+14+DF™M)rQ) 





. (2.1) 


Proof. Let f € T;(v,¢, 4, 9,3,t, 6). Then there exists two holomorphic functions 
u,v:U — U given by 
u(z) = uz t+Ugz?7 +u3z27 ++ (z EU) (2:2) 
and 
vw) = y4wtv,w2+v3wi t+ (we), (2.3) 


with u(0) = v(0) = 0, |u(z)| < 1, |v(w)| < 1, z,w € U such that 


[: ere) nee (Bi) | ; ( 2 etsy) 
Bil (Z) 





Braf@) (88.7) 
=1+G?(Qu(z) + 6} (Qu? (z) + (2.4) 
and 
WV ! " 0) 
w (8.90) ae (85,90) ae (8i,90¥)) 
Big) Big W) (B°, g(w)) 





=14+G62()v(w) + G2 (Hv?) +> . (2.5) 
Combining (2.2), (2.3), (2.4) and (2.5), we obtain 


ny ' Wy 7P 
[ (®/@) aad (Bf) a [ a (Biaf@) } 
Br af) Bi af) (Bo f (z)) 


= 14 Gi(t)u,zt [G2 (uz + Go (tuz|z? +-- (2.6) 
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and 


YY ' Wy 7? 
(" (#2900) [i yt (#2900) 7 ( 0 (Bing) } 
Br,9W) Br 9(w) (w°, g(w)) 


=14+ GP Ovw + [G2 Ov. + G2 Ov{]w? + . (2.7) 
It is quite well-known that if |u(z)| < 1 and |v(w)| < 1, z,w € U, then 
Ju;,|}<1 and |v;|)<1 forall iEN. (2.8) 
Equating the coefficients in (2.6) and (2.7), we deduce that 


O(v+o+D)M+OTO+ DIM +2) 


ly +604+44+1l@Mra) az = GP (tur, (2.9) 


66 +1)(y t+ o(2u+ DY) + O)FQt 2) +9) 
Tq +6+3+20MmrQ) a3 


ere +60)F704+ DPq+[yvy-D+ oust D(2v + @-DUt+D) -2(v¥ + 6Bu+ DY ae 
2F4(7n +6 +44 1)F*(@)r2Q) 


= GP (tu, + G2 (tus, (2.10) 


Oy t+ GH+ DINO + OO + 1PM +3) 


Tn +64+34+DF()rQ) a, = GP(t)y, (2.11) 


and 


66+ V(yt+ o(2u+ D/P + rat 2rm+a 


lq +0+1+2@rQ@) (2a3 — as) 


pee (Qt+ Ort DI2M 4+ Diy -— 1) + $+ D(2v + @-Dw4+1)) -2(yv + But DY oe 
2F4(n +6 +44 1)F*@)r2Q) 


= GP (t)v2 + G2 (v7. (2.12) 
From (2.9) and (2.11), we conclude that 

uy =v; (2.13) 
and 


267(v+ d(ut 1)” T27+ 9204+ 12 + » 2 


r2(7 +6 +34 Dl2(m)F2(3) 2=(G? () (ug + vz). (2.14) 
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Adding (2.10) to (2.12), yields 


er(n + OT] +3) 
r@)rQ) 


where Q(y, ¢, Hn, 8, 4) is given by (2.1). Consequently, we have 


OY, P, w)az = GP (t)(uz + va) + GF) (uz + vf), 











6l(n + OT +3) 262(y + (w+ 1)) T2(y + O20. + 12 +: GE OO], 
Fapra) 27%, ¢ ~) - 
rq +0 +34 DIAm)T2Q) (G9), 
= G?(t)(uz + v2). (2.15) 


Further computations using (1.3), (2.8) and (2.15), we obtain 
2|6|tl7 +6444 D)°Q)rQ)y2/6t 


502(y + o(u +1)) 2 + O20. + 12 (n +3) 


_, | 66+ D627 + oH+D) 12M + O04 Dre +3) 
-08°°0,¢, OF +0rq+9r2q7+64+144+ DI™Mra) 


|az| < 





Next, if we subtract (2.12) from (2.10), we deduce that 


200 + 1)(y+ (2u+ 1D) + OFO+2)r~q +2) 
(y+ 6+3+2)(™MrQa) 


= G2 (O)(uz — v2) + GF (uz — v7). (2.16) 
In view of (2.13) and (2.14), we get from (2.16) 


(a3 — a3) 


r2(n +0 +34 ITAn)I2Q) (GPO). 
~ 262(y + b(t 1)) T2(y + O20 + D2 +2) 


Py t+64+34+2)TMragio 
* 360 +1)(yt+¢o2u+D)rq+OrQt 2m +a 


(uj + v7) 
(uz — V2). 


Thus applying (1.3), we obtain 
4Al?2(n +6 +34 1)7()I72()62t? 
62(y + out 1)) F2(y + O20 + DIZ +) 
2°(n +0414 2) (TQ) |[dlt 
+56 +1(yt+¢o(2u+D)'~q+Ora+2rq +2) 





la3| < 
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Putting @ = 0 and y = 1 in Theorem 2.1, we demonstrate the next outcome: 


Corollary 2.1. For t € G. 1| and 6 is a nonzero real constant, let f € A be in the 
family Ts (y, 6,4, t, 6). Then 
2|6|tl(y7 +04+44+ DI@)ra)y 2[6lt 
66272(n + O24 1)P2(y +.) 


la2| < 





_9 6(6 + 1)67T72(n + 6)F2(00 + 1) 24m +4) Je 
—9626(n, 0,3) + O(n +27 +9 $34 D(A) 
and 
2 4T2(n +6434 1)F2(y)I2(9)82t? 2°(n + 0414+ 2)T()rQ)ldlt 
las] 64r2(n+O)F2004+1274+2) 00641 +OTO+ 204+)’ 
where 


2(0+1)FQ+2) —-20r~™m+6)r?204+1)r@ +3) 


S09.) = Tey eex4a) > Tq +0 4941 FOO) 


Theorem 2.2. ForO0 <o0 <1,t€ G. 1] and 6 is a nonzero real constant, let 


f € A be in the family Ss(o,n, 6,3, t, 5). Then 


Ty +64+344+ DrMPra)|d|ty dle 
662(o + 1)2F2(y + O)F20 + 1)2(y +) 
26(6 + 1)62(o + 1)2F2(y + 6204 + 1)T2(9 +0) 2 
—0657Y(0,n,9,)T(n + OTH +92 +6 +44 IIMA) 


laz| < 





and 


l2(7n +6444 1)F?()P2()67¢? 
~ 62(¢6 +:1)2F2(n + ONT 200 + 1)T 2 (n +: 3) 
2°(n + 6414+ 2)T™@)Ir(a)|dlt 
30(6+1)(20 + 1)I(n + OA)TOA+ 2) (n +)’ 


la3| < 


where 


3(06 + 1)(20 + 1)TQ 4+ 2) 


¥(o,n,6,3) = 
(0,7, 8,4) ly) +0 +342) 
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4060 +1) + O72 +1) +9) 


r4y7+64+44+1)r()ra) ey) 


Proof. Let f € Ss(o,7,6,3,t,6). Then there exist two holomorphic functions 
u,v:U—>U 


mr 


Zz (B8.f@)) + (20 + 1)z? (B8f@) +023 (@9,f@)) 
Zz (B2f@) oz (B8sf@)) 
=14G?(t)u(z) + G2(t)u?(z) ++ (2.18) 


and 


mr 


w (9,9(w)) + (20 + 1)w? (89,90) + ow? (89,9(w)) 
Ww (89,9(w)) + ow? (89,9(w)) 
=1+G62()vw)+Giv?(w) t+. (2.19) 
where u and v have the forms (2.2) and (2.3). Combining (2.18) and (2.19), yield 


mr 


Zz (B8.f@) + (20 + 1)z? (B8.f@) +023 (B9,f@)) 
Z (B8.f@) + oz? (Bef) 


=14+G2(t)uyz + [G2 @uz + GS ()uz|z? + -- (2.20) 





and 


mr 


w (@%,9(w)) + (20 + 1)w? (9,9(w)) + ow? (%,9(w)) 
w (@9,9(w)) + aw? (89,9(w)) 
=14+G%()yw t [G2 Ov. + Gi(Hvi|w2 +. . (221) 
Equating the coefficients in (2.20) and (2.21), we deduce that 
20(06 +1)T + 4)FO+ 1) +4) 


— co 
ln +0434 DMO) a2 = Gi (t)uy, (2.22) 
3016 +1)2o+DIM+6TQ+ 2) +4) 
PQ +e+a+2FQr@ 
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462(6 +1)7T2(y + O20 + DI2(y +3 
rea = GP(t)u, + G3 (uj, (2.23) 


206 +1 + OTO+ II] +9 


— cd 

Tn +0434 DFM)ra) a, = Gy (t)vy (2.24) 
and 
300 + Dot WIM + OTG+ 2+ |, 2 
—__F@+0+34 Dra) 6242 ~ 93) 

467(o +1)?F2(n + OT20 + IT2(y +3 
AEG FO FTE DGG) HOM + HOME (2.25) 
In view of (2.22) and (2.24), we have 
am (2.26) 


and 


862(o + 1)*F2(y + O)T2(4 + 12 +3 
Re = (GP) (uj + v2). (2.27) 


If we add (2.23) to (2.25), we conclude that 


261 (yn + OTC +4) 
r@) ra) 
where Y (a, 7, 6,4) is given by (2.17). 


YC, , 0, das = G2 (t) (uz a V2) + 68 (t) (uz + v2), (2.28) 


By substitute the value of uz + v7? from (2.27) in (2.28), yields 


2orn+ OM D i, og 9) 876+ DP AtOPAt+DPAtYEO] 
,n, 9, 2 
r@ra) P27 +6 +24 1I?()F2Q) (G8). 








= GP (t)(u2 + v2), 


or equivalently 


r2( +0 +34 1P*CH)P2@) (GEC) (uz + v2) 
as = SS OO IO, 
* 26V(0.7,8,90( + OP + F209 + 8 43-4 DEMFO) (ce) 


—862(o + 1)?F2(y + 8)F2(1 + 12 + NGS (4) 


(2.29) 


Further computations using (1.3), (2.7) and (2.29), we obtain 
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lq] +6444 1)@Mra)dlty dle 





la2| < 





602(0 + 1)2T2(n + O)F2(4 + 1)P24(y +3) - 
26(6 + 1)04(0 + 1)7TF4(7 +. 924 4+ 1) 2m +:4) . 
—067Y(0,n,6, MT (n+ OPH +27) + 6414+ 1)FM)rQ) 
Next, if we subtract (2.25) from (2.23), we deduce that 
66(6 + 1)(20 +1) (y+ 6)TO+ 2)TQ +4) 
lyg+64+344+2)rq@ra) 

= GP (t)(u2 — v2) + G2 (O(uz — vf). (2.30) 

In view of (2.26) and (2.27), we get from (2.30) 


(a3 — az) 


r2() +0 +24 DIA@)I2Q) (GP). 
~ 862(o + 1)2F2(@m + O20 + D2 +) 
Pq +4 +3+2rayragr@ 
60(8 + 1)(20 +1) + OO + 2M +3) 


a3 (uf + vz) 


(uz — V2). 


Thus applying (1.3), we obtain 


l4(7n +6 +44 1)F?2(y)P2()67¢? 
62(0 + 1)2F2(n + O)T2(4 + 1) 2( +) 
2°(n +6414 2) Mra) [dle 
30(6+1)(20 + 1) (7 + OA)TA+ 2) (n +4) ° 


la3| < 


Putting o = 0 in Theorem 2.2, we demonstrate the next outcome: 


Corollary 2.2. For t € @ 1 and 6 is a nonzero real constant, let f € A be in the 
family $5(7, 6,4, t, 6). Then 


ry +6 +24 DPM rO|sleyTsle 


la| < 





56202(n + OT 200 + DI2(y +3) 
E 25(5 + 1)62F2(n + O)T2(00 + D2(y +) 2 
—062U(n, 6, OT + OFM +HlF2() +0 +34 DITA) 


and 
T2(n +6444 1)F2()P2(2)67t? 2°(n +6 +314 2)F@)ra)|dle 


Isl <“poreG + Ora + DI +a) | 30+ DI +OlOt DI +d)’ 
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where 


3(0+1)FQ4+2) 467] +6704 1I)F@ +)) 


MOM =Teyera42) Ty+0+1+DF@lQ 


Next theorems, show “Fekete-Szeg6 problem” of the families Ty (y, ¢, u,,7, 6,4, t, 6) 
and Ss(o,, 6,3,t,6). 


Theorem 2.3. For 0<y <1,0<@<1,0<yn<1, te (5,1) &£€ Rand 6 isa 
nonzero real constant, let f € A be in the family Ts(y, ¢, u,, 8,4,t, 6). Then 


lag — Ea3| 
2t|o|T (7 +9 +342) 
fc +1)(vt+ ¢(2u+D)q+ Ort 2)r~q +3)’ 
bs ly+64+14+2)r(@Mra) G 
~ 200+ D(y + o(2u+ 1) + 4)FQ+2)r~™ +3) 








for |§—1| 


507(y + o(u + 1)) T2( + O20. + DF2(n +3) 
= | 55+ )E(v+oH4)) PO+OrE+DrOH+) | ,2 
-08°70,¢6,47,6, 9° (7 + OF(y +27 + 8 4+44 ITM) 
267t7T2(n +0 +14 1)T2()r7Q) 





x 


8t3|632 +0414 DIP()raylé - 1 
562(y + P(u+1)) 12] + O)F2(4 + D2) +) 
Ls | 5(6+162(y + oUt D) P24 O24 D2 +) . 
-0570(7,¢,47,8,)TM+ OPH +27) +6414 1ITQ)rQa) 
T(y+64+34+2)r()ra) 
for |E-1|= x 
2006+ 1)(v+o(2u+ D+ Ort 2r~q+) 








507(y + o(u + 1)) 12 + O20. + DF2(n +3) 


. | 
“5 | 565+ DO%(y + out 1) P+ Ort DIM +9 | 2 
: 08 OA(y, 647,68, DT (yn + OT] +27 4+ 6434+ DI) 
262t272(yn +6 +34 1T2(m)F2() 


Proof. In the light of (2.15) and (2.16), we deduce that 
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Py +6 +34 2TMTOAG?(O (uz — v2) 


--ég2 = (1 -Ha2 | 
az — a5 = (1 02 + 66+ Dy +n + DI + OG+DI +H 


r2(n +8 +34 DTA(m)l20) (GPO) (ue + v2) 
60 + OM + PMO) +8 3+ 107, 6,447,8,3) (G20) - 
202(y + (ut 1)) T2(y + O20 + D2 +: GS(t) 


@ Py + 6 +44 2) MPA)G?P OH) U2 — v2) 
2016 + 1)(v + o(2u +1) P+ OTO+2)rm +2) 


= (1-8) 





2g rn +6+34+2)@)rQ) 
=H (vo - 2010 + 1)(y + ¢(2u+ DI + OTA+ 2) +9) ue 
& ( oo T(7 +6 +44 2)0(y) FQ) ) 
ws 2000+ (y+ o2ut+ DM + OO+ 2) +9 vs 
where 
2 
+6434 DPM) (G20) A-H 
wg) = 


Aq + OPO + MCMPArG] + 6 +1+ YOY, bu, 8,3) (G2) - 
202(y + p(w + 1)) T2(y + O20 + DFA + G8(C) 
According to (1.3), we deduce that 
2t|o|\T (7 +6 +14+2)TM)ra) 
66+ 1I)(yt+ o(2u+ DY) +O)FQ+2)M +3)’ 
Ty +64+34+2)()rQ) 


0< go ee 
la, — &a2| < Ww 200 + 1)(y¥+ out D+ OO +t 20 +3) 
4t/d|lpS)l, 
Tn +04+3+4+2)TQ)ra 
bE) > (7 +O +4+2)TrQ) 


2006 + 1)(y+ O(2n+ 1) + OTO+2)rq +2) 


After some computations, we obtain 
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la3 — éa3| 
2t|6|[[m+6+14+2)TM)ra) 
68+ 1(y+ o2u+ D+ OTO+ 2) +3)’ 
Ty +64+14+2)r™Mra) 








for |F -—1|/< x 
2016+ 1(yt+o(2u+ D+ Ort 2)rq +) 
567(y + P(u+1)) T2] + O)F2(04 + D2) +) 


LS | 56+ NOV +POU+D) PO +OPrO4+DPG+2 | ,2 
08° 0,647.8, 91) + OF +4927 + 6 +44 DIMIQ) 





262t2F2(n +0 +34 DI2(m)r2Q) 





50?(y + (w+ 1)) T2@ + O20. + FA +3) 
Ee | 56+ 1)O(y + PUF1) P+ O24 D242 | 42 
-66707,¢,4.7,6, 9TH + OT +92 +6434 D°M)rO) 
lq +6+44+2)r@)ra) 


for |& — 1| >————_—— *7?r 
20(6 + 1)(y + (2u+ DIM + OAt2rq + 


567(y + P(u+1)) T2 + F204 + D2) +) 
= | 554+ )E(v+oHt)) PO+OrEt+DrO+) | ,2 
-68°07,¢6,47,6,)° (7 + OF(y + )P27 + 8 4+44 ITM) 
262t2F2(n + 9 +14 DI2(m)l2Q) 


: 8t3|63|T2(7 +9 +34 DIZ) )IE — 1 


Putting € = 1 in Theorem 2.3, we demonstrate the next outcome: 
Corollary 2.3. ForO <y<1,0<@<1,0<u<si,teE G. 1] and 6 is a nonzero 
real constant, let f € A be in the family Ts(y, ¢, u,,n, 9,4,t, 0). Then 


2t|6|T@q +8+1+2)r@)ra) 


|a3 — a3 |< 
06 +1)(y + o(2u+ D)P~™ + OPQ + 2) +3) 


Putting @ = 0 and y = 1 in Theorem 2.3, we demonstrate the next outcome: 
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Corollary 2.4. Fort € G. 1], € € Rand 6 is a nonzero real constant, let f € A be 
in the family Ty (n, 6,3,t, 6). Then 


|a3 = Ea3| 
2t|oITq +6+4+2)PQ)rQ) _ 
6(6+1)']+6)TO+2)~ +3)’ 
ly +6414 2)F(M)ra) 


for 8 11 S F564 Drm + Ort DI +) 


66°T7(n + O)F2(4 + 1)F7( +3) 
2g 6(6 + 1)62T2(y + O)F2(4 + 1)72(y +3) 2 
-665°S(n,6,)T (7 + OTH + YF27(7 + 6 +44 1)PM)TO) 
267t7T2(n +0 +44 1)F2()r2Q) 


8t3|63|[2(n +O +4 LA(n)F2(HIE - 1 
602(y + bu + 1)) T2(m + O20 + 1)P2(m +3) 

5 | 5(6 + 1)62(y + Ou + 1)) T2() + O20 + DI2Q] +9) | . 

~9526(n, 0, (yn + Oy + DP2() +9 44 YPM) 
a lg +64+314+2)rQ@)ra) ‘ 
~ 2064+ 1(vy+ o2u+ 1D) +OrO+2r~q+ 
50?(y + o(u+1)) 2 + O20. + 12 (n +3) 
rE | 5(5 + 1)6*(y + (ut. 1)) 12 + O20 + D2 +2) | 2 

~9526(n, 0, (yn + Oy + NP2() +9 44 DPM) 

262t7T2(n +0444 1)F2()rZ2Q) 


for | -1| 


Theorem 2.4. ForO <a <1,te€ = 1}, € € Rand 6 is a nonzero real constant, let 
2 


f € A be in the family Ss(o,n, 6,4, t,6). Then 
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las = Eaj| 
2t|6|\T(m +6 +34 2)TM)ra) . 
30(6 + 1)(20 +1) (+ 8)TO4+ 2)F™ +3)’ 
lq7+@4+344+2)rq@)ra) 


Or SNS 396+ Deat DIG + Ot AIG * 


66*(o + 1)7T2(7 + O)T70 + 1) 2 +9) 
25(6 + 1)67(o +.1)7T7(y + O)F200 + 197 +3) 5 
: ~ 1-062Y(0,7,6, T+ OF + HI2) + 6 34+ DIMA) 
67t772(7n +0 4+44+1)F2()I2Q) 


< 2t3|63|f2 +4434 DF@M)ra)é - 1] 
502(a + 1)212(n + O20 + D2) +) 
E | 26(5 + 1)62(o + 1)2F2(y + O)F2( + 1)F2(y +) : 
—062Y(0,n,9, MT H+ OlH+ F277 +6434 DF)ra) 
(7 +64+344+2)F()rQa) 


for IS 12 3664 Dot DI + OG + DIG +y © 


667(0 + 1)7T2(7 + O)T70 + 1) 2 +9) 
i | 26(6 + 1)02(o6 +1)772(y + O)T274 1927) +) 2 
. —06*Y(6,n,0,)T H+ OPQ +D)F27H +0434 DFMO) 
67t772(7n +0444 1)F2()I2Q) 


Proof. In view of (2.29) and (2.30), we deduce that 


Pn +6 +34 2) M)FO)G2 (t)(u2 — v2) 


a3 — $03 = (1-803 + Fg past Drm + OlAtDIG+tD 


r2(n +8 +24 DTA(m)l2Q) (GPO) (ue + v2) 
26Y(o,n,6, WT (n + O(n + OP2) + 8 434+ DI()LQ) (? (wo) 
—862(o + 1)?F2(y + O)T2(4 + 1) 2H +:1)G9(t) 


Py +6 +34 2)PMTOAGE() (uz — v2) 
60(6 +1)(20 +1) + 4)TO +2) 4+ 3) 


=i lee) 
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_ Go) @) +4 rq +6 +4+2)@ra) 
~~2 \?S2 3664 Dot DI +Orat2r+y) 2 
ry +6 +34 2)@)ra) 
: (ow ~ 360 + D2o+ DI +O lat2rm + 5) v2} 
where 
ji rn +8 +34 DIA@)r2@) (G2) A-H 
Pp = 


6Y(o,9, 8,20) + OP +2 + 6 ++ DOPE) (G50) 
—462(¢ + 1)*F2(7 + O20 + 12] +: GIO) 
According to (1.3), we deduce that 
2t|\o|Tq+6+1+2)r@)ra) 
30(0 +1)(20 + 1)(n + ATO+ 2) +3)’ 
0<|e®|< lq7+0@4+344+2)rq@ra) 


cere 30(6 + 1)(20 + DT + OO + 20 +2) 


2t|d|le@)I, 


io@)) = lq +64+14+2)r@ra) 


30(0 + Do+ Dry t+ OlFAt2I7 +d 


After some computations, we obtain 
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las = Eaj| 
2t|o|T( + 6 +44 2)()FQ) . 
30(6 + 1)(20 +1) (+ 8)TO4+ 2)F™ +3)’ 
Tn +0434 2)0()TQ) 
for |€ — 1| << ——_——  —  —— — — —  _ 
30(8+1)(20+ 107 + 4 FO+ 2) (n+) 
66*(o + 1)7T2(7 + O)T70 + 1) 2 +9) 
“ 25(6 + 1)67(o +.1)7T7(y + O)F200 + 197 +3) 2 
: -06*Y(6,n,0,)T H+ OPQ +D)F27H +0 +34 DFTA) 
67t772(7n +0 4+44+1)F2()I2Q) 


2¢3/63|T2q +6 +34 D2) rz — 1 
§62(o + D2F2q7 + OF20 + DFA) +) 
: 26(6 + 1)62(o + 1)2F2(y + BIZ + DI2(y +.) . 
—062Y(0,n,9, MT H+ OlH+ F277 +6434 DF)ra) 
rn +6434 2)()FQ) 
300+ )@o+DIt+Orat2arn+y 


IA 


for | -—1| = 


667(0 + 1)7T2(7 + O)T70 + 1) 2 +9) 
_ | 25(6 + 1)67(o +.1)7T7(y + 8)F204 + 1972 +3) 2 
: -065*Y(6,n,0,)T H+ OPQ +)F27yH +0434 DFTA) 
67t772(7n +0444 1)F2()I2Q) 


Putting € = 1 in Theorem 2.4, we demonstrate the next outcome: 
Corollary 2.5. For0<a0<1,teé G. 1 and 6 is a nonzero real constant, let f € A 
be in the family Ss(o,n, 9,4, t,6). Then 


2tl6|Tq+68+144+2)TQ)rQa) 
30(6 + 1)(20 + 1) (7 + OA)TO4+ 2) (7 +:9) 





las = a3| < 
Putting o = 0 in Theorem 2.4, we demonstrate the next outcome: 


Corollary 2.6. For t € G. 1], € € Rand 6 is a nonzero real constant, let f € A be 
in the family $5(n,0,3,t,6). Then 
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las a Ea3| 
2tjo|ly+6+14+2)r~M)ra) 
36(6 +1) (7 + 4)FO+2)(y +3)’ 
ly +0@4+14+2)rq@)ra) 


for I$ — S354 pr@ + Oat DI ty 


667F2(n + O)F2(4 + 127) +: 9) 
26(6 + 1)67F2(n + O)204 + D247 +.) ; 
—9067U(n,8, NTH + OH +)F2(7 +6 +34 DIMI) 


- §2¢2T2(n +6 +24 DI2r2) 
< 2¢3|63|[2(n +9 +34 1)04(n)lZ2)é — 1 
607T2(y + @)F20 + 1)F2(y +4) 
aly 26(6 + 1)02F2(n + O)F2(4 + 1)P2(y +3) Je 2 
662U(n, 6TH + OH +HP2(7 +9 +34 IEMA) 
; nee Ty +0+31+2)@) ra) 
or lk 2 S564 Dra + Oratara +a 
667F2(n + O)F2(4 + 1727 +: 9) 
_ | 26(6 + 1)67F2(y + 8)T2(4 + IT727() +: 9) 5 
—057U(n, 9, NT) + OTH + HPA) + 6 +34 DIMA) 
67t772(7 + 6 +44 1)F2(y)T2Q) 
Conclusion 


The primary objective was to create the families Ts(y,¢,u,,7,0,1,t,6) and 
Ss(o,n, 8,4,t, 5) of bi-univalent functions which governed by Gegenbauer polynomials. 
We generated Taylor coefficient inequalities for functions in the families 
Ts(v,¢,4,7,9,3,t,6) and Ss(o,n,6,3,t,6) and viewed the famous Fekete-Szegé 


problem. 
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Abstract 


In this work a novel distribution has been explored referred as Burhan distribution. This 
distribution is obtained through convex combination of exponential and gamma 
distribution to analyse complex real-life data. The distinct structural properties of the 
formulated distribution have been derived and discussed. The behaviour of probability 
density function (pdf) and cumulative distribution function (cdf) are illustrated through 
different graphs. The estimation of the established distribution parameters are performed 
by maximum likelihood estimation method. Eventually the versatility of the established 
distribution is examined through two real life data sets. 


1. Introduction 


In biomedicine, engineering, economics, and other fields of science, statistical 
distribution plays essential role in modelling data. The exponential and gamma 
distributions are popular distributions for assessing statistical data and are regarded as life 
time distributions. Among these distributions, the exponential distribution has one 
parameter and several attractive statistical characteristics, such as memory less and a 
constant hazard rate. Various extensions of these distributions have been made in the 
Statistical literature to lead greater flexibility. Lindley formulated a one parameter life 
time distribution with the following probability density function in (1958). 


2 


rarer Oe 2 ee y>0,0d>0 





fQ,8) = 
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In recent years, researchers have worked extensively upon Lindley distribution, 
formulating one and two parameter distributions. For the modelling of diverse complex 
data, Ghitney et al. [2] conducted an in-depth study on the Lindley distribution, 
demonstrating that the Lindley distribution outranks the exponential distribution in terms 
of waiting times for bank customers. They also demonstrate that the contours of the 
hazard rate function of a Lindley distribution are increasing whereas the mean residual 
life is decreasing. Zakerzadeh and Dolati [3], Nadarajah et al. [4], they extended Lindley 
distribution with addition of new parameters and expounded the performance of the 
extended distribution through data sets. In recent years many authors have made different 
contributions to modify the Lindley distribution. Merovci [5], has introduced transmuted 
Lindley distribution and discussed its several properties. Sharma et al has introduced the 
inverse of the Lindley distribution and propounded its different properties. Shanker et al 
[6] developed a new life time distribution named Akash distribution which were proved 
superior then exponential and Lindley distribution. Sen et al. [7] has proposed a one 
parameter distribution called xgamma distribution and they have proved by an 
application that xgamma distribution provide better fit than exponential distribution. K.K 
Shukla [8], have suggested the Pranav distribution and studied its different properties. 
Aijaz et al. [9] has introduced transmuted inverse Lindley distributions and estimates its 
properties. Aijaz et al. [10] has developed a new distribution and name it Hamza 
distribution and studied its different properties. When it comes to analysing more 
complex data, these distributions have both advantages and limitations. In this study, an 
attempt is made to establish a new three parameter distribution that is considerably more 
pliable and produces better results than previous ones. The newly defined three parameter 
distribution’s probability density function is as follows. 


gBt1 yb 
FOB.) = ape reget ee y>0,a,B,0> 0. (1.1) 


The suggested distribution is referred to as the Burhan distribution, and it has been 
developed by combining two distributions, exponential (@) and gamma (f,8), by 
employing the linear combination method. 


fy, a, B,8) = pgidv,8) + 1 — p)g2, 8,4), 


gBt1 


Where p= aoP+r(By 





B 
910,89) = 6e-%: y>0,0>0; go,9)= RUE y>0,p,0>0. 
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Figures (1.1) and (1.2) expound few layouts of pdf of Burhan distribution for varying 
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Figure 1.1: pdf of Burhan distribution Figure 1.2: pdf of Burhan distribution 
under different values to parameters. under different values to parameters. 


The associated cdf of equation (1.1) is given by 


(apo e~°Y+r(p+1,6y)) 


F(y,a,P,@) =1- : >0,68,0>0. 1.2 
(y, a, B, @) p(aoP+r(B)) y B (1.2) 


Figures (1.3) and (1.4) expounds few layouts of cdf of Burhan distribution for 
varying parameters 
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Figure 1.3: cdf of Burhan distribution Figure 1.4: cdf of Burhan distribution 
under different values to parameters. under different values to parameters. 
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2. Mathematical Properties of Burhan Distribution 
2.1. Moments of Burhan distribution 
h 


Let Y be a random variable follows Burhan distribution. Then r*” moment denoted 


by yu; is given as 


He =EW")=| y"FO.48,8)dy 


=| y'Fo.a8, 8 ay 


00 yh 
[reefer 
0 B 
= «| ye” v4ef yT*B e- 8 dy, 
0 B Jo 


After solving the integral, we have 


: — apokr(r +1) + T(r +B +1) 
Br portBet 


Substituting r = 1,2,3,4 we obtain first four moments about origin as 
, a6? + (6 +1)F(B) 
M1 = oB+2 


, 2008 + (6 +2)(6B +1)F(f) 
Ha = pees 


, 600% + (6 +3)6+2)6+DI(h) 
M3 = BHA 


, 24008 +(B+4)B4+3)6+2B6+ DF) 
= = — = pies = = 





2.2. Moment generating function of Burhan distribution 


Let Y be a random variable follows Burhan distribution. Then the moment generating 


function of the distribution denoted by My (t) is given 


My (t) = E(e™) [ ero. a, B,0) dy. 
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Using Taylor’s series 


=| (1404S +& t...) F008. Ody 


o t’ co 

-Y 5] y’ f(y, a, B, 8) dy 
r=0 — 9 
so rT 








t 
=o ae 
r=0 
tT apo’r(r+Dt+r(r+B+D 
My(t) = y oe Bortb+1 : 
r=0 


3. Renyi Entropy of Burhan Distribution 
If Y denotes a continuous random variable having probability density function 


f(y, a, B, 8). Then Renyi entropy is defined as 


1 foe) 
Ta = z=; 09 [“Prondy, 





where y > 0 andy #1. 
Thus, the Renyi entropy of Burhan distribution (1.1) is given as 


1 g(B+Dy 7 aa ee 
A ee ee “_\ o-vovg 
RY) = 7— y cals +I(B))" Jo (« *B 5 | 


1 gbt+Dr gy 0 yB\" Z 
= ——_ log ;\—_—___, (1+2) eer dye. 
1-y ~((a0% + rp)" Jo B 


Using generalized binomial theorem, we have 


soa att Ree ae ae Gee 
TRY) = 60) (sr aa. yP e vy dy}. 


T= 





After solving the integral, we get 


2, yy ( cot) \" r(Br +1) 
NE ayy D3 (7) (a n re) apy Wort 
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4. Tsallis Entropy of Burhan Distribution 


Tsallis entropy of order y for Burhan distribution (1.1) is given as 


1 ioe) 
s=—1- | roddy} 
. al ‘A aie 
where y > 0 andy #1 


1 gbrDy pay BY 
Sy =——})1-—____, (2+2] e-Y9Y dy). 
(a6? +(B)) Jo 


After solving the integral, we get 
ee (B+1) \" 
eae a 1-Y) EO ES See es Sy 
es a! a: ab® +(B)) (ap) (ya)er*4 
T= 
5. Mean Deviation from Mean of Burhan Distribution 


The quantity of scattering in a population is evidently measured to some extent by the 
totality of the deviations. Let Y be a random variable from Burhan distribution with mean 


pL. Then the mean deviation from mean is defined as 


D(u) = EY — ul) 


s I IY - ul FO)ay 
Lu ee) 
=| w-nrordy+ | o-wfoday 
0 Lu 
in in 00 oo 
n[r0) y [ 970) y+ | sf0) y [F009 y 
L ee) 
=nFW)-[ yfo)dy-utt-Fwl+ | yforay 
0 Lu 


L 
= 2uF(u) 2 i yf) dy. (5.1) 


Now 


vi 7 u gBbt1 yh apy 
[ sroray= | orem z) dy 
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ae [ 6 Le B+1,-0 
= ————_ ja e "Yd +2] e "Yd \ 
wa Rea f Bia 4 


After solving the integral, we get 
B+1 


a 1 
rere Ou) + GopraVB + 2, n6)}. (5.2) 


Lu 
[ sro) y= BoB? 


Using equation (5.2) and (1.2) in equation (5.1), we obtain 


2u (apore-% + (B+ 1,6n)) 


ae B(a0® +1 (B)) 


298+1 


a 1 
* Pew TTC y(2, Ou) + aaa VB + 2, 16) 


BO 
6. Mean Deviation from Median of Burhan Distribution 


Let Y be a random variable from Burhan distribution with median M. Then the mean 
deviation from median is defined as 


D(M) = E(|Y — M|) 


= ly — MI fay 
0 


M ioe) 
= (M—y) foray + [ (y — M) f)dy 


M 


M foo) 
~ MF(M) - I yf(y) dy — M[1 — F(M)] + ip yf”) ay 


M 
=f I yf (”) dy. (6.1) 


Now 


M M gBt1 yh nis 
[ xroray=| yore t ee dy. 
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After solving the integral, we get 
fa) Bt+1 


a 1 
war Foy gz”? OM) + pea VB + 2,mo)}. (6.2) 


M 
i yf) dy = Bobe 


Using equation (6.2) and (1.2) in equation (6.1), we obtain 


208+1 


DM) eR aT) 


a 1 
{ozv(2,6M) + pore’ B + 2,m0)} 


7. Reliability Measures 


Suppose Y be a continuous random variable with cdf F(y), y = 0. Then its reliability 
function which is also called survival function is defined as 


S(y) = p,(¥ > y) = | f(y) dy =1- FO). 
y 


Therefore, the survival function for Burhan distribution is given as 
SO, a, B, 0) =1- F(y,a, 8,0) 


(apo’e-% ere ed, ay)) 


a 7.1 
B(a0* + T(B)) ae 
The hazard rate function of a random variable y is given as 
fa, B, 8) 
0) = ———_.. . 7.2 
h(y, @, B, ) S(y, a, B, 0) ( ) 


Using equation (1.1) and equation (7.1) in (7.2), we get 


Boat (a + )e% 
h(y, a, B,@) = (apofe-% +T(B +1,0y)) 


Figure (7.1) and (7.2) expound few layouts of hazard rate function of Burhan distribution 
for varying parameters. 
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Figure 7.1: HRF of Burhan distribution Figure 7.2: HRF of Burhan distribution 
under different values to parameters. under different values to parameters. 


The cumulative hazard rate function is given as 


H(y,@,B, 0) = —log|F(y, @, B, )]. 


Using equation (1.2), we get 


(apo? eV 4T(B + 1,6y)) 
B(ao8 + F(B)) 


Mean residual function of random variable y can be obtained as 





H(y,a, 8,0) = —log|1— 


m(y, a, 8,0) = ae ome tf (t,a,B,0) dt — 


SQy,a 


— B(ao8 +78) pe ght ie 
~ Gpore-97 +B +1,0y) J, a0P+T(B) (« oF 5) aa | 


= ees Ac [ ter at + aI tte tat - 
apore-Y +T(B +1,0y) Wo BJ, . 


After solving the integral, we get 
apo’ (@y + 1)e—°” + F'(B + 2,0y) 
6(apoFe-% +T(B + 1,6y)) 
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8. Order Statistics of Burhan Distribution 


Let Y(4), ¥(2), Y(3),-+-+» Ym) denotes the order statistics of a random sample drawn 
from a continuous distribution with cdf F(y) and pdf f(y), then the pdf of Y(,) is given 
by 


frFroy 411 — Fy) "* k = 1,2,3....,. (8.1) 





n!} 
fro. ®) = Gaara! 


Substitute the equation (1.1) and (2.2) in equation (8.1), we obtain the probability 


function of k“" order statistics of Burhan distribution is given by 








2) i f (apo%e + (8 +1, ay" 
e- _— 


7 n! gB+1 
fren 9) (« ar B(a6F +I(B)) 


(k-— 1)! (n—k)! a6 +1(f) 


(apo e-® +r(B+1,6y)) Bok as 
p(aoF+r(p)) ae 


Then, the pdf of first order statistics Y(3) of Burhan distribution given by 


noet ( *) ie (apobe-% +P (B+ 1,6y)) 
J a ee ee eh 


aoP+T@\* B B(aoP +1 (By) 


fray, 8) = B 


And the pdf of nth order statistics ¥(,) of Burhan distribution is given by 


fray, 8) = 


1-k 
nobti yh (apore-% +T(6 + 1,6y)) 
e~9y 1 8 ———— EEE EE ee 
B 


aoP+TBy\ * B B(ad® +T(B)) 


9. Maximum Likelihood Estimator of Burhan Distribution 


Let Y,,Y2,..,Y, be a random sample of size n from Burhan distribution. Then its 
likelihood function is given by 


l= | [roweb.0) 
i=1 


nn 


B+1 ; 2 
~ (arsrH) [teres 


i=1 





http://www.earthlinepublishers.com 


Burhan Distribution with Structural Properties and Applications ... 439 





The log likelihood function is given as 


logl=n(6 +1) log 6 —nlog (a0? + r(B)) + ih, log (a + 20) 0 1%; - (9-1) 


The partial derivatives of equation (9.1), with respect parameters are given as 


n 








dlogl _—noF B ne 
da ~ GOP +1) Leap +yP oe 
dlogl _ n (ao? log @ +I'(B)) — (B log yi — 1)y8 
dlogl_nB+1) napoht 
00. 6 abP+r(B) pa ae 


From equations (9.2), (9.3) and (9.4), we have obtained a system of non-linear equations 
which cannot be expressed in compact form and is difficult to solve explicitly for a, 6 
and @. Applying the iterative methods such as Newton-Raphson method, secant method, 
Regula-falsi method etc. The MLE of theparameters denoted as ¢ (a, B, 6) of ((a@, 8,8) 


can be obtained by using the above methods. 


For interval estimation and hypothesis tests on the model parameters, an information 
matrix is required. The 3 by 3 observed matrix is 


0? logl 07 logl a7 logl 

fe: af EE Vy geal es 
0a? dadBp dade 

1 0? logl a7 logl 0? logl 
i@2=-le ot \ # ot). F d 
n| \ dBaa ap? apae 

a7 logl d7 log l 0? logl 
E a E ae E = 
000a O00p 002 


The elements of above information matrix can obtain by differentiating equations (9.2), 




















(9.3) and (9.4) again partially. Under standard regularity conditions when n > oo the 
distribution of ¢ can be approximated by a multivariate normal N (0, I (¢ i) distribution 
to construct approximate confidence interval for the parameters. 


Hence the approximate 100(1—w)% confidence interval for a, and @ are 


respectively given by 
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at 2y |tad(Q). B+ Ze lig) and 0+ 2» |158(0). 


where Zy denotes the ¢ th percentile of the standard normal distribution. 
2 


10. Applications 


This section exhibits the adaptability of the formulated distribution by applying real- 
world data sets. The suggested distribution is compared to the Shanker distribution 
(SHD), the Akash distribution (AD), the Ishita distribution (ID), the Pranav distribution 
(PD), the Rani distribution (RD), the Prakaamy distribution (PKD), and the Lindley 
distribution. We were efficient in outranking the specified distribution. 


In order to compare the above distribution models, we consider the criteria like AIC 
(Akaike information criterion), CAIC (corrected Akaike information criterion), BIC 
(Bayesian information criterion). Among the above distributions, the better distribution is 
considered having lesser values of AIC, CAIC, HQIC and BIC. 


Data set 1: This data is obtained from Murthy et al. [13], which represents failure 
times of 84 Aircraft Windshield. The data follows: 


0.040, 1.866, 2.385, 3.443, 0.301, 1.876, 2.481, 3.467, 0.309, 1.899, 2.610, 3.478, 0.557, 
1.911, 2.625, 3.578, 0.943, 1.912, 2.632,3.595, 1.070, 1.914, 2.646, 3.699, 1.124, 1.981, 
2.661, 3.779, 1.248, 2.010, 2.688, 3.924, 1.281, 2.038, 2.82,3, 4.035, 1.281,2.085, 2.890, 
4.121, 1.303, 2.089, 2.902, 4.167, 1.432, 2.097, 2.934, 4.240, 1.480, 2.135, 2.962, 4.255, 
1.505, 2.154, 2.964, 4.278,1.506, 2.190, 3.000, 4.305, 1.568, 2.194, 3.103, 4.376, 1.615, 
2.223, 3.114, 4.449, 1.619, 2.224, 3.117, 4.485, 1.652, 2.229, 3.166,4.570, 1.652, 2.300, 
3.344, 4.602, 1.757, 2.324, 3.376, 4.663. 


Table 10.1: Descriptive statistics of data 1°. 





Min Q, | Median | Mean] Q3 | Skew. | Kurt. | Max 





0.040 | 1.866 | 2.385 | 2.563 | 3.376 | 0.086 | 2.365 | 4.663 
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Table 10.2: The ML estimates of the unknown parameters for data set 1™. 








































































































Model a B ) S.E 

@ p f 
BHD | 0.0206 | 5.6641 | 2.4893 | 0.0129 | 1.0566 | 0.4044 
SHD | sess . |esates 0.6440 | ...... | eee eee 0.0459 

AD ~ |.eutees  \|eSeites O.9337) |isessesr || ee we'ene 0.056 
ID — | |easasvet || regenes 0.9181 |... | ce ee ee 0.0495 
PD |uceseder |] evans 1.2104 | wc... | ceeeee 0.0530 
RD > | cegewe jl evouas LD1DDe «| |seadesey Ill evens 0.0563 
PKD | 2.0... | eee Le7287 ° \\seeeae. dill ncesss 0.0710 
ED. 5 nseeee. — |teweee’s 0.6296 | ...... | coeeee 0.0502 

Table 10.3: Performance of distributions for data set 1“. 

Model —logl AIC CAIC BIC HQIC 
BHD 127.63 261.26 261.55 268.59 264.21 
SHD 151.55 305.10 305.14 307.54 306.08 

AD 147.46 296.93 296.97 299.37 297.91 
ID 146.83 295.66 295.71 298.10 296.64 
PD 144.40 290.81 290.85 293.25 291.79 
RD 143.12 288.24 288.29 290.68 289.22 
PKD 188.70 379.40 379.44 381.84 380.38 
LD 153.96 309.93 309.98 312.37 310.91 








Data set 2: The data set represents the survival times (in days) of 72 guinea pigs 


infected with virulent tubercle bacilli, observed and reported by Bjerkedal [12]. The data 


are as follows: 
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0.1, 0.33, 0.44, 0.56, 0.59, 0.59, 0.72, 0.74, 0.92, 0.93,0.96, 1, 1, 1.02, 1.05, 1.07, 1.07, 
1.08, 1.08, 1.08, 1.09, 1.12, 1.13,1.15, 1.16, 1.2, 1.21, 1.22, 1.22, 1.24, 1.3, 1.34, 1.36, 
1.39, 1.44, 1.46,1.53, 1.59, 1.6, 1.63, 1.68, 1.71, 1.72, 1.76, 1.83, 1.95, 1.96, 1.97,2.02, 
2.13, 2.15, 2.16, 2.22, 2.3, 2.31, 2.4, 2.45, 2.51, 2.53, 2.54, 2.78,2.93, 3.27, 3.42, 3.47, 
3.61, 4.02, 4.32, 4.58, 5.55, 2.54, 0.77. 


Table 10.4: Descriptive Statistics of data 2" 





Min 


Q 


Median 


Mean 


Qs 


Skew. 


Kurt. 


Max 








0.100 





1.077 





1.450 


1.754 








2.240 





1.328 


4.913 








5.550 








Table 10.5: The ML estimates of the unknown parameters for data set 2". 





































































































Model a B 6 S.E 
@ B 6 
BHD 0.0067 2.4143 1.9124 0.0085 0.6517 0.3729 
SHD iste leases 0.8667 | .seeee | ese 0.0675 
AD? © | veesine "|, senses LQ2220 5 |i deesse. — ioetivews 0.0815 
IDs = |iviwese, — __ [asewe 1.1664." | vestes 2 lls eee 0.0680 
PID 2 iseusee: = Wseedies 1.4780 | .cseee | wee 0.068 
RD* + © | Seee. —— ili devas B7994s. — | veseee: < —_|heeese's 0.0699 
PKD | ...... | ew eee AAVTB. —vasvec, |lsdetes 0.0853 
ED? =I venveam > | esses 0.8744 | eee | eee 0.0771 
Table 10.6: Performance of distributions for data set 2™. 
Model —logl AIC CAIC BIC HQIC 

BHD 94.06 194.12 194.48 200.95 196.84 

SHD 105.98 213.97 214.02 214.87 216.24 

AD 107.06 216.13 216.19 218.41 217.04 

ID 108.03 218.06 218.12 220.34 218.97 

PD 112.53 227.07 227.12 229.34 227.97 

RD 118.26 238.53 238.58 240.80 239.43 

PKD 176.87 355.75 355.81 358.03 356.66 

LD 106.52 215.05 215.10 217.32 215.95 
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Figure 10.1: Fitted pdf’s for data set I. 
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Figure 10.2: Fitted pdf’s for data set II. 


It is clear from Table 10.3 and 10.6 that Burhan distribution (BHD) has least values 
of - log, AIC, CAIC, BIC and HQIC when compared with competitive distributions. We 
accomplish that Burhan distribution provides an adequate fit than compared distributions. 
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11. Conclusion 


In the present paper a novel distribution has been formulated by employing convex 
combination of exponential and gamma distribution stated as Burhan distribution (BHD). 
Its several properties including moments, moment generating function, survival function, 
hazard rate function, mean residual life function, mean deviations, order statistics, Renyi 
entropy and Tsallis entropy have been discussed. The parameters of the distribution have 
been estimated by known method of maximum likelihood estimator. Finally the 
performance of the model has been examined through two data sets and compared which 
shows that Burhan distribution gives an adequate fit for the data sets. 
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